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The goal of this thesis effort was to develop an 
interactive electromagnetic scattering algorithm for 
personal computer use in the analysis and computer aided 
design (CAD) of scattering structures. Extensive 
interactive computer aided design type features have been 
added to a previously developed finite element computer 
algorithm which considers scattering from inhomogeneous 
axisymmetric objects. These computer aided design features 
allow structures composed of multiple dielectric and/or 
magnetic materials to be graphically constructed or 
modified by the user by way of either a keyboard or mouse 
interface. Bistatic scattering patterns are graphically 
presented upon completion of the computation for a variety 
of user specified incident field orientations. 
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I. 



INTRODUCTION 



The study of resonance region electromagnetic 
interaction with arbitrarily shaped and often inhomogeneous 
penetrable objects has increased dramatically in the past 
few years. This new concentration of effort is driven by 
the likelihood of expanded implementation of high frequency 
over-the-horizon radars, very high frequency satellite-borne 
radars and medium frequency surface wave radars. Such 
radars have the potential to observe "stealth" type 
platforms whose low radar cross section (RCS) design was 
limited to more conventional radar frequencies, above 500 
MHz. In addition, these new studies of electromagnetic 
interaction have proven to be of benefit in other areas 
such as propagation through rain, power absorption by 
biological tissues, focusing by microwave lenses, radiation 
from reentry vehicles and electromagnetic pulse (EMP) 
coupling into vehicles, installations and weapon systems. 
[Refs. 1 and 2] 

As frequencies become sufficiently low, structures 
having a low radar cross section may exhibit strong 
resonance characteristics. Design of such structures is 
complicated by these low frequency resonance 
characteristics, due to the lack of accurate analytical 
methods. Even attempts at structure design utilizing 
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"shaping" and "edge treatments" are ineffective at providing 
a highly reduced radar cross section at resonant 
frequencies. At higher frequencies, "asymptotic" 

formulations are available to rectify these problems [Ref. 
3], A more comprehensive approach is needed, one which 
allows the designer to rigorously evaluate the interaction 
of the scatterer with the illuminating field. At the same 
time, such a tool should be highly interactive with the 
user. 

The subject of this thesis is the development of a 
generalized, interactive, finite element scattering 
algorithm to be used in the analysis and computer aided 
design of scattering structures. In this, the first stage 
of effort, the programs developed for use on a personal 
computer calculate electromagnetic scattering by 
axisymmetric inhomogeneous structures. These structures are 
composed of user specified lossy or lossless penetrable 
(i.e., dielectric or magnetic) materials. The resulting 
programs provide significant improvement in computational 
speed and versatility over the current integral equation 
based computer codes being used. The physical structure can 
be quite general, other than the requirement for rotational 
symmetry. This general construction avoids the 

computational restrictions found in surface integral 
equation formulations, where, as frequency and geometric 
complexity increase, the time required to obtain a solution 
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increases. The finite element algorithms, based upon a 
volume variational formulation, are relatively immune to 
structural complexity driving the computation time. Another 
advantage of these algorithms is their method of matrix 
manipulation. The finite element algorithms produce very 
"sparse" matrix structures, whereas the integral equation 
formulations produce "full" matrices [Ref. 4]. Finally, 
when considering the construction of realistic material 
configurations, the computational time required utilizing 
the finite element algorithms can be substantially less than 
the time required by the use of integral equation 
formulations [Ref. 5]. 

The development of the axisymmetric finite element 
program is based upon Morgan's previous effort [Ref. 1], 
which restricted structures to non-metallic configurations. 
Major programming modifications to improve computational 
efficiency and to incorporate computer aided design features 
such as interactive graphics are now integral parts of the 
finite element computer program. 

A brief overview of the finite element scattering 
algorithm, which constitutes the basis of this thesis, is 
presented in Chapter II. In addition to the appropriate 
background information, the concept of coupled azimuthal 
potentials, the finite element algorithm and the unimoment 
method are reviewed. 
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Having established, in Chapter II, the necessary theory 
utilized in the construction of the original version of the 
computer program EMSCAT, Chapter III presents methods to 
construct program input data. Included in this chapter are 
discussions of computer programs written specifically for a 
given geometric shape, data file input and keyboard input. 
Basic flow charts provide an overview of the input-output 
interactions of the all computer programs. 

The first half of Chapter IV discusses the extension of 
the methods used to construct program input data to those 
methods that utilize graphics. Of special interest is the 
incorporation of computer aided design features. The second 
half of this chapter addresses the user interactive output 
graphics features provided. These features allow the user 
to obtain one of four different graphs, either on the 
computer screen or printer, of the EMCAD output data of 
choice. 

In Chapter V, the validation of modifications and 
improvements to EMCAD is presented. Graphical comparisons 
are made between the output of the modified version of 
EMSCAT versus both theoretical calculations and actual 
experimental data. In addition, validation of new programs 
written specifically to be used in conjunction with EMCAD is 
presented. 

The brief conclusions in Chapter VI will summarize the 
work previously presented, as well as provide 
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recommendations for future modifications to the interactive 
package of computer programs being used. 

The Appendices provide copies of the source code for 
all the computer programs utilized in this interactive 
package. A sample listing of the output of EMCAD is 

attached. As an added feature, user manuals with hints, 

warnings and directions for all of the major computer 
programs are included. 
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II. FINITE ELEMENT SCATTERING ALGORITHM 



A. BACKGROUND 

In recent years, great attention and effort has been 
directed toward the attempt to obtain the numerical 
electromagnetic scattering solutions for inhomogeneous 
material structures. Early efforts dealt primarily with 
improving the efficiency of existing computational 
algorithms. These algorithms utilized surface integral 
equation formulations, which required the inversion of large 
’‘full" matrices. Accordingly, methods using these 
formulations were computationally slow while being highly 
restricted in the level of inhomogeneity that could be 
handled. [Refs. 1 and 2] 

A major advance in matrix manipulation occurred with 
the development by Mei [Ref. 6] of the unimoment method. 
This new method decoupled the interior and exterior regions 
associated with radiation or scattering problems, resulting 
in uniformly banded system matrices for the interior region 
solution. The primary advantage of uniformly banded 
matrices is the capability to invert them directly using 
fast banded matrix inversion techniques, such as the Riccati 
transform. Mei's breakthrough opened the door for the 
development of applications such as the finite-difference 
solution of radiation by a biconical antenna, and the 
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finite-element solution of scattering by two-dimensional 
inhomogeneous dielectric cylinders. [Refs. 1 and 2] 

Attempting to improve upon these applications, Morgan, 
Chang and Mei [Ref. 7] conducted research with the goal of 
extending the unimoment technique to problems involving 
three-dimensional, axisymmetric inhomogeneous structures. 
During the initial period of this research, another major 
breakthrough occurred. This breakthrough was the discovery 
of the coupled azimuthal potential (CAP) formulation 
[Ref. 7]. The CAP formulation became the basis for 
construction of several numerical algorithms. These 
algorithms incorporated the capability to perform scattering 
computations for inhomogenous axisymmetric metallic and 
penetrable structures [Ref. 2]. 

Utilizing a triregional unimoment method, these 
numerical algorithms were incorporated into a versatile 
finite-element computer program. The triregional unimoment 
method was chosen because of its ability to minimize the 
bandwidth of the global system matrix, and its ability to 
minimize the number of nodes in the finite-element mesh. 
The finite-element computer program generates scattering 
solutions for an arbitrarily shaped inhomogeneous lossy 
penetrable body of revolution being impinged upon by 
multiple incident fields. [Ref. 2] 
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B. COUPLED AZIMUTHAL POTENTIALS 

The solution of electromagnetic field problems 
involving generalized axisymmetric penetrable structures 
using the classical TE and TM potential formulations are 
restricted to a handful of configurations, and have at most 
inhomogeneities in one-dimension. By utilizing the coupled 
azimuthal potential (CAP) formulation, solutions are valid 
in lossy inhomogeneous isotropic rotationally symmetric 
media, as shown in Figure 2-1, where the complex material 
parameters are invariant to the azimuthal coordinate. 
[Ref. 1] 

Utilizing normalized circular cylindrical coordinates, 
as also shown in Figure 2-1, the first step in developing 
the necessary equations is to decompose the total 
electromagnetic field into azimuthal modes through the use 
of an exponential Fourier series, 

co 

E(R,Z S §) = X ^(RjZ) exp(jmf) (2.1) 

m = -co 



and 



© 

n„H(R,Z,§) = 7 h ( R ) Z ) exp ( jm§) (2.2) 

u ro 

m = -co 
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NORMALIZED CIRCULAR CYLINDRICAL COORDINATES 



where 



(R,Z,§)- = (kgP.kgZ.f) 



(2.3) 



with 



k 



0 




(2.4) 



and 



Tig = 1 2 0 TT J2 



(2.5) 



Two uniformly continuous coupled azimuthal potential 
functions 



t (RiZ)fn) and V ( R , Z > m ) 
1 2 
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can be used to generate the modal electromagnetic field 
vectors, 



e ( R * Z ) 
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where the two-dimensional gradient operator is 



? = 





( 2 . 8 ) 



while the multiplicative function is given by [Refs. 1, 2, 
and 7] 



f (R,Z) 
m 



[ji (R,Z)€ (R,Z)R‘ 
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( 2 . 9 ) 
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For a two-dimensional cut in the (R,Z) plane of Figure 
2-1, the coupled azimuthal potential functions satisfy a 
formally self-adjoint system of second-order partial 
differential equations. Utilizing either Euler-Lagrange or 
stationary-theorem approaches, it can be shown that there is 
an easily managed variational principle which can be applied 
in addition to the partial differential equation 
formulation. 

This variational formulation is based upon the 
existence of a stationary function which manifests itself as 
a surface integral over the planar cross section, S, 



F 



"i p 

dR dZ 

J J S 



( 2 . 10 ) 



where the Lagrangian L, which is unique to within an 
arbitrary constant multiplier and arbitrary independent 
additive function, is of the form [Refs. 1 and 2] 
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C. FINITE-ELEMENT ALGORITHM 

Having presented an overview of the coupled azimuthal 
potential formulation, it is now appropriate to explain its 
application to the electromagnetic scattering problem. 
Figure 2-2 illustrates the finite-element mesh geometry 
employed by the numerical algorithm. The (R,Z) planar 
configuration shown is named the semi-annular conformal 
(SAC) mesh (Ref. 1]. This mesh is comprised of numerous 
linear triangular elements, which conform to the surface of 
the axisymmetric scattering body. This surface conformation 
simplifies analysis of normal derivative discontinuities at 
the surface boundary. The semi-annular conformal mesh 
provides a distinct computational advantage to those 
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FIG. 2-2. SEMI-ANNULAR CONFORMAL 
FINITE-ELEMENT MESH FOR AN 
INHOMOGENEOUS BODY OF REVOLUTION 
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scattering structures which can support a fairly large 
homogeneous spherical core [Ref. 2]. 

Figure 2-2 graphically demonstrates the regions of 
interest when calculating the solution of the scattering 
problem. The homogeneous spherical core occupies Region I, 
the inhomogeneous cross sectional surface of the 
axisymmetric penetrable scattering body is in Region II, and 
an unbounded exterior region is in Region III. The surface 
of interest is shown with the capability of being both 
banded and layered. The layers are depicted as "films" that 
run the entire length of the surface, while bands are 
depicted as isolated "pockets". The possibility exists, as 
shown, that all the bands and all the layers have differing 
material parameters, thus establishing a structure that is 
truly inhomogeneous. Note that the surfaces of Region I and 
Region III bound the surface of interest located in Region 
II. These surfaces, r = a and r = b, can therefore be used 
to specify Dirichlet boundary conditions [Ref. 2]. The 
numerical solution for the electromagnetic fields in the 
finite-element domain (Region II) can be computed for each 
of a set of boundary conditions. By requiring continuity of 
the electromagnetic fields at the interfaces where r = 
and r = r 2 , the unknown expansion coefficients for the core 
and scattered fields can be obtained [Ref. 1]. This is the 
essential concept behind the unimoment method. 
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In Region II of Figure 2-2, the serai-annular conformal 
mesh is composed of nodes, with each interior node 
surrounded by six triangular regions, known as "elements". 
The azimuthal field components of the mesh are described by 

N 

§ • e < R , 2 ) = 7 e ( n ) u < R , Z ) (2.12) 

m u m n 

n= 1 

and 

N 

§ • h (R,Z) = 7 h ( n ) u (R,Z) (2. 13) 

m ^ m n 

n = 1 

where u n (R,Z) represents the area coordinate linear pyramid 
basis functions. Each of these basis functions has linear 
variation in the six triangular elements around each node, 
while maintaining a value of one at the central node and a 
value of zero at all bounding nodes. Each set of e m (n) and 
h m (n) represent one of the N complex nodal values found in 
the azimuthal field. To initiate the variational finite- 
element technique, the linear shape-function expansions in 
Equations 2.12 and 2.13 are substituted into the surface 
integral stationary function F described in Equation 2.10. 
Imposing the null derivative condition with respect to each 
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of the unknown coordinates representing the interior nodal 
fields, while constraining the coordinates of the boundary 
nodal fields, allows the computation of the stationary 
point of the function F. The resulting linear system for 
the nodal fields is described by 
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where n = each internal node. When these equations are 
combined with the semi-annular conformal finite element mesh 
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analysis, each equation, for a given interior node n, will 
relate e m (n) and h m (n) to 14 adjacent node azimuthal field 
values e m (i) and h TO (i). Some of these adjacent field values 
may be specified on a boundary. The result of conducting 
these calculations for each value of n is a sparse global 
system matrix, containing complex elements and displaying a 
diagonal tri-block submatrix structure. Applying the 
Dirichlet boundary conditions to the semi-annular conformal 
finite-element mesh shown in Figure 2-2, the numerical 
solution for the electromagnetic fields in Region II is 
obtained by utilizing a two-sweep block-by-block matrix 
inversion algorithm related to the Riccati transform for 
difference equations. The diagonal tri-block submatrix 
structure is given by 



fi. * 

k 



k-1 



Wk * c kVi 




( 2 . 16 ) 



A^, Bfc, and Cfc are sparse submatrices which constitute 
the tri-block diagonal structure of the global system 
matrix. is generated from Equations 2.14 and 2.15 

through the use of the Dirichlet boundary conditions. 
[Ref. 2] 
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The mesh construction and matrix generation are 
conducted on a local basis in a sequential manner. The 
analysis starts at the positive Z-axis, and works clockwise 
in equal ©-increments. Using a simple minimum distortion 
scheme at each step, the mesh is constructed to conform to 
multiple approximate piecewise smooth surfaces. The Riccati 
transform sparse matrix inversion technique conducts a 
"first sweep" and a "backsweep". On the "first sweep", the 
values of R^ and Sj^ are calculated, where 

C K (2.17) 
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k k k 
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(2.18) 



On the "backsweep", the values of the solution vector are 
obtained, where the solution vector is given by [Refs. 1 and 
2 ] 
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k - 1 






R, 

k k 










(2.19) 
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A principal advantage of this numerical method is that 
the Riccati transform technique provides a systematic 
procedure. This procedure divides a larger area into 
several smaller areas defined by the tri-block matrix 
substructure, existing between any three adjacent radial 
rows of nodal values. Analysis using this type of smaller 
area procedure allows all operations including mesh 
construction, matrix element generation and matrix 
inversion to be performed at the same time. [Ref. 2] 

D . UNIMOMENT METHOD 

The unimoment method is a computational procedure which 
allows numerical values to be calculated for the complex 
coefficients of the incident, scattered and spherical core 
electromagnetic field expansions when dealing with the semi- 
annular conformal finite-element mesh. The scattering 
solution is obtained by representing the electromagnetic 
fields of Region I and Region III, as shown in Figure 2-2, 
in terms of TE and TM spherical harmonic expansions, with 
unknown core and scattered field coefficients. Using 
superposition of boundary conditions and interior solutions 
for the azimuthal fields, the core expansion is established 
as a boundary condition along r = a while the scattered 
field expansion added to the incident field is established 
as a boundary condition along r = b. The two existing sets 
of boundary conditions are modal spherical harmonic core 
fields along r = a, with zero fields along r = b, and 
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spherical harmonic and incident fields at r = b, with null 



fields at r = a. The multiple boundary conditions will 
generate multiple solutions at the nodes where r = r^ or r = 
r 2 , (i.e., along the boundaries). The method used to obtain 

the expansion coefficients sets the weighted superposition 
of interior solutions to the appropriate original field 
expansions. Along the boundaries at r = and r = r 2 , this 
method enforces equality utilizing least squares analysis. 
[Refs. 1 and 2] 

This method of analysis is only valid if the scattered 
electromagnetic fields everywhere outside of the interior 
mesh region can be represented by field harmonics. These 
representations appear as infinite series, and must be 
mathematically truncated to facilitate numerical 
calculations. Decomposing the plane-wave incident fields 
into azimuthal modes yields 



js i n ( mf) 
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( 2 . 21 ) 
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1 if in = 



where e m = 



0, e m = 2 if m > 1; M is the truncation 



index, and the polarization convention (TE/TM) , depicted in 
Figure 2-3, is used in the bracketed sinusoids. [Ref. 2] 

During computations, the numerical scattering solution 
algorithm is run concurrently with the sequential Riccati 
transform finite-element solution in Region II. On the 
forward sweep, boundary conditions for both the incident and 
modal fields are generated at each node row being used in 
that particular sweep step. On the backsweep, boundary 
conditions along the r = r^ and r = r 2 contours are used to 
compute each nodal value. The incident field and spherical 
harmonic field for each node is formed, and used to generate 
the nodal residual fields. During the backsweep, inner 
product integrations are formed, and these inner product 
results are added to the field moment matrix, T, and the 
incident field driving vector, D corresponding to each 
incident field. The scattering solution is obtained by 
multiplying the inverse of the field-moment matrix with the 
driving vector to obtain the field expansion coefficients. 
The complex scattering amplitude f(9,$) is defined in terms 
of the far-zone E-field, 
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FIG. 2-3. POLARIZATION CONVENTION 



where d max is the maximum scattering body dimension and 
where the components of f(8,$) are given by 
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The scattering cross section is found to be [Refs. 1 and 2] 
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Utilizing the theoretical principles and mathematical 
calculations just discussed, Morgan [Ref. 1] wrote the 
computer program EMSCAT to compute electromagnetic field 
scattering from layered bodies of revolution. Extensive 
interactive computer aided design features have been added 
to EMSCAT by Morgan and Connolly. These computer-aided 
design features, driving the change of name from EMSCAT to 
EMCAD, allow structures composed of multiple dielectric 
materials to be graphically constructed or modified by the 
user, and then evaluated by EMCAD for scattering. The 
remainder of this thesis is concerned with the utilization 
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and application of EMCAD, as well as with modifications and 
improvements of this transformation. 
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III. EMCAD INPUT DATA CONSTRUCTION 



As has been previously explained, EMCAD computes 
electromagnetic field scattering from penetrable bodies of 
revolution using the coupled azimuthal potential formulation 
in conjunction with a variational finite-element technique 
and a tri-regional unimoment method. The numerical solution 
is performed using a two-sweep Riccati transform with 
circumferential marching. In order to provide relatively 
"fast”, consistently accurate answers, yet be flexible 
enough to allow the user to vary any of a number of input 
parameters as conditions warrant, EMCAD requires sufficient 
input data to adequately depict the . penetrable body of 
revolution, as well as the desired incident and scattered 
field orientations of the user. This accuracy and 
flexibility is accomplished through the utilization of 
keyboard input and/or mouse input to form both material and 
structure input data files. Mouse input deals with the 
graphical interface, and as such will be discussed in 
Chapter IV. 

A. KEYBOARD INPUT 

Information obtained by keyboard input is fairly 
diversified in nature. This data ranges from input and 
output file names, to desired mesh parameters, from 
selection of phi angles, to selection of incidence angles. 
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To facilitate ease in understanding and running the program, 
EMCAD allows a user to input the necessary information 
according to the level of expertise needed. For the 
"Novice" level, the program gives brief explanations of 
input values, describes the required format, and gives 
typical values where applicable. For the other option, the 
"Expert" level assumes the user is familiar with the input 
parameters and formats, and simply prompts for the required 
inputs. A more detailed explanation of each facet of the 
keyboard input is contained in Appendix A, EMCAD USER'S 
MANUAL. 

B. DATA FILE INPUT 

Once the information obtained by keyboard input is 
established, the remaining data becomes that which describes 
the penetrable body of revolution being considered. This 
descriptive data can be separated into two distinct areas of 
information. The first area of information is the data 
which conveys the material composition of each layer of the 
structure. Call this the "material file". The second area 
of information is the data which defines the structural 
shape of each layer of the penetrable body, referred to as 



the "structur 


file" . 


As shown in 


Figure 


3-1, 


there 


are 


three methods 


of data 


construction 


which 


can 


be used 


to 


generate the 


"Material" 


file and the 


"Structur" 


file . 


The 



quotation marks indicate that the data files can be named at 
the user's discretion. Method 2 and Method 3 utilize 
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SCATTERING CALCULATION 
INPUT-OUTPUT RELATIONSHIP 



computer aided design features, and will therefore be 
discussed in Chapter IV, Interactive Graphics Features. 

Data construction by way of Method 1 is depicted in 
Figure 3-2 and is composed of special case programs which 
generate data files for specific geometrically shaped 
objects. For example, program SPHERE allows the user to 
define a layered penetrable sphere, offset from center if 
desired, with as many penetrable layers as can be handled by 
EMCAD at this time. The user enters material dielectric 
parameters, layer radii, and required resolution as 
prompted, and the program writes the necessary "Material" 
and "Structur" data into files in the required EMCAD 
format. CONE and CYLINDER were written specifically for 
EMCAD. As with SPHERE, they allow the user to establish the 
body of revolution as a multilayered structure. Source code 
for SPHERE, CONE, and CYLINDER can be found in Appendices J, 
K, and L respectively. Whenever a computer program 
generates data and data files, the capability exists for the 
data to be generated manually, or by any other means 
available, and then entered into the appropriate data files 
manually through the use of the keyboard. Keyboard entry is 
shown as a viable method of data construction in Method One, 
Two and Three. 
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DATA SCATTERING. 
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METHOD NUMBER ONE 



Spheres, cones, and cylinders are simple geometric 
shapes, and are extremely useful for elementary 
applications; however, nature is complex, and not always 
made up of simple geometry and elementary shapes. To obtain 
realistic results, realistic penetrable bodies are needed. 
It would indeed be tedious to attempt to write a structure 
defining computer program for each new shape encountered. 
As such, the optimum solution indicates the need for the 
implementation of computer aided graphics. 
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IV. INTERACTIVE GRAPHICS FEATURES 



Interactive graphics simply refers to the ability to 
take data at hand, and depict that data in some visual form, 
where a graph is the most commonly utilized depiction. A 
special subset of interactive graphics, known as computer 
aided design, also allows the user to draw a picture on the 
screen, and convert that picture to data. 

EMCAD was written in the programming language of Ryan- 
McFarland FORTRAN (Appendix M) . Microcompatibles GRAFMATIC 
(Appendix M) , and later Jewell Technologies GRAFPLUS 
(Appendix M) were chosen as the commercially available 
screen graphics and hard copy support for EMCAD output 
graphics. Because the output format that was used readily 
allowed conversion to an acceptable format, West Coast 
Consultants CURVE DIGITIZER (Appendix M) was chosen as the 
computer aided design package. 

A. COMPUTER AIDED DESIGN 

There are three means of constructing the data for 
penetrable bodies of revolution, as previously shown in 
Figure 3-1. Method 1 has already been discussed, and is the 
predecessor to Method 2. In Method 2, the "Material” data 
can be generated by one of the three programs written 
specifically for EMCAD, or by direct keyboard input into a 
data file if desired. The change from Method 1 occurs in 
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the construction of the structure defining data. As shown 
in Figure 4-1, the "Curveout" data is obtained by using 
CURVE DIGITIZER. In simple terms, this means that the user 
draws the picture on the screen using a mouse with the 
computer aided design package. Once the user is satisfied 
with the single, or multilayered drawing on the screen, the 
data is saved in the format established by the computer 
aided design package. An example of the two-dimensional 
cross section of a three-dimensional body of revolution 
which was designed using CURVE DIGITIZER is shown in 
Figure 4-2. 

The data file that is saved is a valid description of 
the structure, but is not in a form which can be used by 
EMCAD. This format conversion is accomplished by a program 
written specifically for this purpose, EMCADIN. The 
computer aided design output data is written in X,Y format, 
with special codes interspersed throughout. EMCAD requires 
radius format and a normalizing factor of 2 n must be taken 
into account. The output data is not necessarily evenly 
spaced, whereas EMCAD requires evenly spaced input data. It 
is conceivable that the data is of insufficient resolution 
for use by EMCAD. Primitive data should be of high 
resolution (i.e., 361 points). The program EMCADIN 
rectifies all of these problems, as well as computes the 
minimum and maximum radius values required. The output of 



34 



DATA SCATTERING 
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FIG. 4-1. SCATTERING CALCULATION 

METHOD NUMBER TWO 



X 
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FIG. 4-2. EXAMPLE OF TWO-DIMENSIONAL CROSS 
SECTION OF BODY OF REVOLUTION DESIGNED 
USING CURVE DIGITIZER 



EMCADIN is the data file "Structur". 



A more detailed 



description can be found in Appendix D, EMCADIN USER'S 
MANUAL. 

Method 3, as shown in Figure 4-3, is actually a union 
of Method 1 and Method 2. "Material" data is generated as 
described in both previous methods. Structural data is 
generated as in Method 1, but is then imported to the 
computer aided design package for modifications. The data 
to be modified must be changed by CRVDXGIN, written 
especially for this purpose, into a format which can be 
accepted as input by CURVE DIGITIZER. CRVDIGIN is explained 
in greater detail in Appendix F, CRVDIGIN USER'S MANUAL. 
Once the data has been read into the computer aided design 
package, the remainder of the development is essentially 
Method 2 . 

B. INTERACTIVE OUTPUT GRAPHICS 

Having constructed the "Material" and the "Structur" 
input data files using one of the three methods discussed, 
it is now possible to conduct the scattering calculations. 
Figure 4-4 shows that the execution of EMCAD yields output 
data and output graphics. 

As also shown in Figure 4-4, EMCAD generates four data 
files which are used as input to the graphics program 
MCADGRAF. It is significant to note that the extensions of 
these filenames have special meaning. TEP refers to TE 
incidence (Figure 2-3), F-Phi; TET refers to TE incidence, 
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FIG. 4-3. SCATTERING CALCULATION 

METHOD NUMBER THREE 



SCATTERING GRAPHICAL 
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FIG. 4-4. SCATTERING CALCULATION 

OUTPUT AND GRAPHICS 



F-Theta; TMP refers to TM incidence (Figure 2-3), F-Phi ; TMT 
refers to TM incidence, F-Theta. 

MCADGRAF is a menu driven output graphics program. 
Appendix H, MCADGRAF USER'S MANUAL, gives examples of all 
the different input menus. The menu structure gives the 
user some general information concerning the use of the 
program, as well as user choice of aspect ratio, data file 
type selection, data file name input, incidence angle 
selection, and phi angle selection. Once these selections 
have been made, MCADGRAF produces four output graph types, 
giving the user the capability to make a printout of any or 
all of these. The first graph type is magnitude versus 
theta in a rectangular coordinate system, as shown in Figure 
4-5. The second graph type is phase versus theta in a 
rectangular coordinate system, as shown in Figure 4-6. The 
third graph type is magnitude versus theta in a polar 
coordinate system, shown in Figure 4-7. The final type is 
simply an axis rotated version of the third type, for a 
different perspective, as shown in Figure 4-8. 
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Offset Sphere Validation Run 
Scattered Field Conplex Anpl itudes 
IE Incidence = 90. Degrees Phi = 0. Degrees 

F-Phi Magnitude Versus Theta 

Magnitude 




Hardcopy — > P or p 



FIG. 4-5. MCADGRAF GRAPHICAL OUTPUT 
MAGNITUDE VERSUS THETA 
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Composite Body 

Scattered Field Conplex Amplitudes 
IH Incidence = 135, Degrees Phi = 180, Degrees 
Phase F- Theta Phase Versus Iheta 




Theta 



FIG. 4-6. MCADGRAF GRAPHICAL OUTPUT 
PHASE VERSUS THETA 
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Offset Sphere Ualidation Run 
Scattered Field Conplex Amplitudes 
IE Incidence = 90, Degrees Phi = 0. Degrees 

Polar Plot, F-Phi Magnitude Uersus Theta 
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FIG. 4-7. MCADGRAF GRAPHICAL GUTPUT 
MAGNITUDE VERSUS THETA, POLAR PLOT 
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Conposite Body 

Scattered Field Conplex Anplitudes 
IH Incidence = 135. Degrees Phi - 180. Degrees 
Polar Plot, F- Theta Magnitude Versus Theta 
Axis Rotated Minus Ninety Degrees 



S Theta=90 




Theta = 0 Theta=180 
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FIG. 4-8. MCADGRAF GRAPHICAL OUTPUT 
MAGNITUDE VERSUS THETA. POLAR PLOT 
AXIS ROTATED MINUS NINETY DEGREES 



44 



V. VALIDATION 



In order to transform the existing program EMSCAT to 
the modified version entitled EMCAD, it was imperative to 
ensure that both the accuracy and the reliability of EMCAD 
were at least as good as that of EMSCAT. It was also 
essential to establish flawless compatibility between EMCAD 
and any program written to support EMCAD. 

The initial step used to validate the operation of 
EMCAD consisted of a comparison between the results obtained 
from EMCAD and the exact theoretical values of a Mie series 
far-field complex amplitude E-plane pattern for an offset 
lossless (e r = 2.6) penetrable dielectric sphere. Run 
parameters chosen included broadside incidence (incidence 
angle = 90 degrees), sphere radius equal to 0.8 (as a factor 
of wavelength), and z-axis offset = 0.16 (as a factor of 
wavelength). The series truncations at M=7, N M I= 8 and 
Nj^ II =15 were chosen for the variable mesh scheme. These 
values were the same as those used by Morgan [Ref.l] in his 
initial validation of EMSCAT. The results of this 
comparison are shown in Figure 5-1 for a phi angle of 0 
degrees, and in Figure 5-2 for a phi angle of 90 degrees. 

Having demonstrated the close resemblance between the 
results obtained from EMCAD, and those predicted by theory, 
the next step used to validate EMCAD required a comparison 
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Offset Sphere Validation Run 

Incidence Angle — 90 Degrees 
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FIG. 5-1. OFFSET SPHERE SCATTERING (#1) 




spn;ju5o^ 
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FIG. 5-2. OFFSET SPHERE SCATTERING (#2) 



between the results obtained from EMCAD and actual measured 
data. Morgan [Ref. 1] validated the original version of 
EMSCAT with extensive experimental verification. Having the 
results of this experimental verification readily available 
presented the unique opportunity to validate the operation 
of EMCAD with the same experimental data used to validate 
EMSCAT. For each axisymmetric body constructed, parameters 
chosen included TE and TM incidence with incidence angles of 
0, 45, 90, 135 and 180 degrees. The series truncations at 
M=9, N m i= 15 and N M II= 20 were chosen for the variable mesh 
scheme. All experimental bodies were constructed of 
plexiglass (polymethylmethacrylate) with material parameters 
defined by e r = 2.59 - j0.017. [Refs. 1 and 2] 

The decision was made to demonstrate the versatility 
and reliability of the program by utilizing a number of 
significant shapes which could be readily fabricated for 
experimental measurements. Since a sphere had already been 
used, and shown to work, a bicone was chosen as an 
appropriate figure. A representation of the solid bicone 
utilized is depicted in Figure 5-3. The comparison for the 
solid bicone validation is shown in Figure 5-4. The hollow 
bicone of Figure 5-5 demonstrates increasing complexity. 
The sphere shown inside the bicone is a free space cavity in 
the object, although it could be made of any material. The 
comparison for the hollow bicone is presented in Figure 5-6. 
The next figure chosen was that of a hollow cylinder. The 
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SCATTERING CALCULATION VALIDATION 
(COORDINATES IN TERMS OF WAVELENGTH) 



So!id Bicone Validation Run 

Incidence Angle « 45 Degrees 
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SOLID BICONE SCATTERING 
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Hollow Bicone Validation Run 
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FIG. 5-6. HOLLOW BICONE SCATTERING 



internal spherical cavity had the same dimension as that 
found in the hollow bicone. The height of the cylinder was 
1.9 (in terms of wavelength) and the radius of the base of 
the cylinder was .945 (in terms of wavelength). The 
validation comparison for a phi angle of 0 degrees is 
depicted in Figure 5-7 and the validation comparison for a 
phi angle of 90 degrees is presented in Figure 5-8. 

To this point, all validation runs have correlated 
extremely well with the experimental data, and all 
experimental figures have consisted of combinations of 
simple geometric shapes such as spheres, cones, and 
cylinders. These are extremely useful for elementary 
applications, however, nature is complex, and not always 
made up of nice geometry or elementary shapes. To prove 
that EMCAD could handle more realistic bodies, the composite 
structure of Figure 5-9 was constructed. Since a program 
was not written to generate the "structur" data required for 
EMCAD input, the decision was made to utilize CURVE 
DIGITIZER (CAD package) to produce the required information. 
The representation of the composite body of Figure 5-9, 
required as input for EMCAD, is depicted in Figure 5-10. 
Note that this representation is composed of data defining a 
two dimensional axisymmetric body. The actual comparison 
for this structure is shown in Figure 5-11. As was the case 
before, the data generated by EMCAD correlates extremely 
well with the experimental data. 
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COMPOSITE SCATTERING BODY 
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FIG. 5-10. COMPOSITE BODY USED FOR 
SCATTERING CALCULATION VALIDATION 
(COORDINATES IN TERMS OF WAVELENGTH) 



Composite Body Validation Run 
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In addition to demonstrating the validity of the 
present version of EMCAD, the composite body of Figure 5-9 
helped demonstrate program compatibility. Since the output 
data of CURVE DIGITIZER was not in the format required for 
use by EMCAD, the computer program EMCADIN was developed and 
written to transform the output data of CURVE DIGITIZER into 
a format which could be used by EMCAD. Figures 5-12, 5-13 
and 5-14 present figures which were generated by the 
programs SPHERE, CONE, and CYLINDER, transformed by CRVDIGIN 
into proper format, and imported into the computer aided 
design package CURVE DIGITIZER. The pictures shown are the 
graphical representations of the imported data as seen by 
CURVE DIGITIZER. This generation and transformation not 
only shows the validation of the programs written to 
generate specific geometric shapes, but also demonstrates 
the compatibility of these programs with CRVDIGIN, and the 
compatibility of CRVDIGIN with CURVE DIGITIZER. 



59 



MULTI-LAYERED 



4 " 

N 



CO 


CM 










cr 


cr 


cr 


LU 


LU 


LU 


>- 


>- 




< 


< 


< 




! 







LU 

LO CT 
ZD LU 
— JZ 
□ Q_ 
< lo 
cr 

i i — ■ 

h LU 
I— LO 
_J Ll 
ZD Ll 
S CD 



o 

(\ J 

r—i 

I 

LO 



LD 

I } 

Ll 



60 



GRAPHICAL REPRESENTATION OF OUTPUT 
DATA OF PROGRAM SPHERE 



MULTI-LAYERED. MULT I -HEIGHT 
MULTI-APEX-ANGLE CONE 
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GRAPHICAL REPRESENTATION OF OUTPUT 
DATA OF PROGRAM CONE 



MULTI-LAYERED, MULT I -HEIGHT 
MULTI-RADIUS CYLINDER 
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GRAPHICAL REPRESENTATION OF OUTPUT 
DATA OF PROGRAM CYLINDER 



VI. CONCLUSIONS 



This thesis has described a very useful numerical 
algorithm for computing resonance region electromagnetic 
interaction with arbitrarily shaped and even inhomogeneous 
axisymmetric penetrable objects. The background material 
associated with the development of this interactive, finite 
element scattering algorithm was presented in Chapter II. 
Included in this background information was the concept of 
coupled azimuthal potentials, the finite element algorithm 
and the unimoment method. Chapter III then presented 
procedures for constructing input data for the scattering 
program. Computer programs written specifically to generate 
a given geometry, and to allow data file input and keyboard 
input were examined in detail. The first half of Chapter IV 
discussed the extension of the methods used to construct 
program input data to those methods utilizing graphics. Of 
special interest was the incorporation of computer aided 
design features. The second half of the chapter presented 
user interactive output graphics features. Validations of 
the modifications and improvements to EMCAD were presented 
in Chapter V. Extensive graphical comparisons were made 
between the output of the modified version of EMSCAT, and 
both theoretical calculations and actual measured 
experimental data. The impressive accuracy of the 
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comparisons validated the changes and modifications to 
EMSCAT used in developing EMCAD. Inter-program 
compatibility was also extensively demonstrated. The 
Appendices provide detailed users' manuals for four of the 
major computer programs developed. 

This was only the first stage of a multistage project. 
Future developments of the axisymmetric code will proceed 
with major programming modifications to allow metallic 
surfaces to be embedded in the inhomogeneous penetrable 
material of the scattering object. Also, the near fields of 
the object will be made available for display. This is 
useful for identifying "hot spots" and field caustics. 

A second effort has already begun to develop a fully 
three-dimensional scattering object capability. The basis 
of this development was the discovery of a generalized 
coupled azimuthal potential formulation which does not 
require an axisymmetric material configuration [Ref. 8]. 
Such a computer algorithm will ultimately employ three- 
dimensional graphics in a powerful computer aided design 
package for user design of scattering objects as well as the 
output of scattering patterns and near field quantities. 
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APPENDIX A 



EMCAD USER'S MANUAL 

Welcome to EMCAD. This program computes 
electromagnetic field scattering from layered bodies of 
revolution using the coupled azimuthal potential formulation 
in conjunction with a variational finite element technique 
and a tri-regional unimoment method. The numerical solution 
is performed using a two-sweep Riccati transform with 
circumferential marching. 

A . BACKGROUND 

EMCAD is based upon an earlier code named EMSCAT which 
was written during the period 1974-1976 by Dr. M. A. Morgan 
while at the University of California at Berkeley. 
Extensive interactive computer aided design (CAD) type 
features have been added in this effort. These CAD features 
allow structures composed of multiple dielectric materials 
to be graphically constructed or modified by the user, and 
then evaluated by EMCAD for scattering. EMCAD has been run 
in both Ryan-McFarland FORTRAN (Appendix M) and Microsoft 
FORTRAN (Appendix M) , with the current version using 
Microsoft. Microcompatibles GRAFMATIC (Appendix M) and 
later Jewell Technologies GRAFPLUS (Appendix M) were chosen 
as the commercially available screen graphics and hard copy 
support for EMCAD output graphics. Because its output 
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format allowed easy conversion to an acceptable format, West 
Coast Consultants CURVE DIGITIZER (Appendix M) was chosen as 
the computer aided design package. Any questions concerning 
colors, formats, changes to the program, etc... should be 
referred to the respective user's manuals, Dr. M. A. Morgan 
or LT. E. M. Connolly (Appendix M) . 

B. COMPUTER SETUP 

EMCAD was written for the following computer 
configuration: 

* IBM AT or compatible, including 80386 based 
compatibles 

* hard disk 

* EGA graphics card 

* EGA monitor 

* math coprocessor - not required, but 
highly recommended 

* memory - 64 OK RAM minimum 

EMCAD can be placed in any directory, at the user's 
discretion, however, if a RAM-drive is to be used, is highly 
recommended that during each run the user allocate a minimum 
of 4M of RAM for EMCAD and EMCAD input/ output files. EMCAD 
input is taken from a variety of sources. The data files 
from these sources are required to be in the same 
directory/sub-directory as EMCAD. 
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C . PROGRAM INPUT 



EMCAD requires both keyboard input and data file input. 
As shown in Figure A-l, two input data files, "MATERIAL" . DAT 
and "STRUCTURE" . DAT , are required as input to EMCAD. The 
quotation marks indicate that the data files can be named at 

the user's discretion. These two data files contain 

information which describes the penetrable body of 
revolution being considered. The "material file" contains 
the data which conveys the material composition of each 

layer of the structure. The "structure file" contains the 
data which defines the structural shape of each layer of the 
penetrable body. The three methods which are used to 

generate these data files are discussed in much greater 
detail in both Chapter II and Chapter III. 

D. PROGRAM EXECUTION 

EMCAD is a menu driven program. Once the correct 
directory has been selected, EMCAD can be invoked by typing 



EMCAD 

at the DOS prompt. This is shown in the top line of Figure 
A-2. The program will load itself, and the user will see 
the next two statements shown in Figure A-2 requesting the 
names of the "material file" and the "structure file". 
Sample keyboard inputs are shown for these two requests. It 
must be pointed out that the extensions are included (.dat) 
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SCATTERING CALCULATION 
INPUT-OUTPUT RELATIONSHIP 



CAIHESISSMSFQRDEMCAD 



ENTER MATERIAL PARAMETER FILE (D: FILENAME. EXTENSION) 
ndatl.dat 



ENTER STRUCTURE DATA FILE (D'.FILENAME. EXTENSIONS 
sdatl.dat 

Offset Sphere R=8.8 D=0. 16 
Er=(2.6,0.) Ur=(l.,0.) 



PLEASE PRESS AN¥ REV TO CONTINUE. 



FIG. 



A-2. 



INPUT DATA FILES 



as required by EMCAD. The two lines beginning with "Offset 
Sphere..." refer to two label lines which were previously 
inserted into the "material file" to allow the user to 
differentiate this set of data from every other set of data. 
By choosing labels carefully, the user can discern from the 
very first screen whether the proper set of data is being 
loaded. This becomes very important because as figures 
increase in complexity, run time for the program increases. 
It is not uncommon to have runs which take several hours on 
an AT class computer. It would be very disturbing to run 
the program for seven or eight hours only to find out that 
the incorrect set of data had been used. Upon completion of 
this screen, the user depresses ENTER to see the following 
screens: 

* Welcome To EMCAD (Figure A-3) 

This screen presents some brief background information. 

* EMCAD Level of Expertise (Figure A-4) 

To facilitate ease in understanding and running the program, 
EMCAD allows a user to input the necessary information 
according to the level of expertise needed. For the 
"Novice" level, the program gives brief explanations of 
input values, describes the required format, and gives 
typical values where applicable. The "Expert" level 
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****** WELCOME TO EMCAD ****** 

THIS PROGRAM COMPUTES EH FIELD SCATTERING FROM PENETRABLE 
BODIES OF REUOLUTION USING THE COUPLED AZIMUTHAL POTENTIAL 
(CAP) FORMULATION IN CONJUNCTION WITH A VARIATIONAL FINITE- 
ELEMENT TECHNIQUE AND A TRI-REG I ONAL UNIMOMENT METHOD. THE 
NUMERICAL SOLUTION IS PERFORMED USING A TWO-SWEEP RICCATI 
TRANSFORM WITH CIRCUMFERENTIAL MARCHING. 

UPDATED VERSION OF ORIGINAL U.C. BERKELEY CODE. 
MODIFICATIONS BY M.A. MORGAN MAR 1987 - FEB 1988 

E.M. CONNOLLY JUL 1987 - APR 1988 . 



PLEASE PRESS ANY KEY TO CONTINUE. 



FIG. A-3. WELCOME TO EMCAD 
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a**** ** ENCAD INPUT ******* 



ENCAD ALLOWS A USER TO INPUT THE NECESSARY INFORNAIION 
ACCORDING TO THE LEVEL OF HIS EXPERTISE. 

1 NOVICE LEVEL - GIVES BRIEF EXPLANATIONS/DESCRIPTIONS 

OF INPUT VALUES. DESCRIBES THE REQUIRED 
FORNAT FOR THE INPUT VALUE. GIVES 
TYPICAL VALUES WHERE APPLICABLE. 

2 EXPERT LEVEL - ASSUNES THE USER IS FAMILIAR WITH THE 

INPUT PARANETERS AND FORNATS. SINPLY 
PRONPTS FOR REQUIRED INPUTS. 

3 EXIT ENCAD 

PLEASE SELECT THE LEVEL OF EXPERTISE BY ENTERING 1, 2, OR 3. 



FIG. A-4. EMCAD LEVEL OF 
USER EXPERTISE 
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assumes the user is familiar with the input parameters and 
formats, and simply prompts for the required inputs. The 
remaining menus will use the "Novice" level to cover all 
ranges of user. 

* Input Name for Output Data File 

(Figure A-5) 

This menu prompts the user for the filename, WITHOUT 
EXTENSION, that will be used to hold the output data of 
EMCAD. The program will append the extension .OUT to the 
user's filename. 

* Input Name for Output Graphics Data Files 

(Figure A-6) 

EMCAD produces four output data files for input to the 
graphing routines of MCADGRAF . All four data files have the 
same name, but each has a distinct extension related to the 
type of data the file contains. This menu prompts the user 
for the name of the output graphics data files, WITHOUT 
EXTENSION. This filename may be the same as that used for 
the output data. 

* Graphics Caption (Figure A-7) 

This menu affords the user the opportunity to input a 
graphics caption to help identify this set of data from all 
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THE OUTPUT DATA FILE IS THE FILE CONTAINING THE OUTPUT 
RESULTS OF ALL NUNER1CAL CALCULATIONS CONDUCTED BY ENCAD. 
THE FORMAT FOR THIS INPUT IS FILENAME ONLY. NO EXTENSION IS 
REQUIRED OR DESIRED. EMCAD AUTOMATICALLY APPENDS AN 
EXTENSION OF .OUT TO YOUR FILENAME. 

PLEASE ENTER THE FILENAME OF THE OUTPUT DATA FILE. 



FIG. A-5. INPUT NAME FOR 
OUTPUT DATA FILE 
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THE OUTPUT GRAPHICS -DATA FILE IS THE FILE CONTAINING THE 
OUTPUT DATA FROM ENCAD TO BE USED AS INPUT TO GRAPHING 
ROUTINES. IHE FORMAT FOR THIS INPUT IS FILENAME ONLY. NO 
EXTENSION IS REQUIRED OR DESIRED. ENCAD AUTOMATICALLY 
APPENDS AN EXTENSION OF .TNT, .TUP, .IET, AND .TEP TO YOUR 
FILENANE AS IT PRODUCES FOUR OUTPUT FILES FOR GRAPHICS. 

.TNT ) TN INCIDENCE, F-THETA 

.IMP ) TN INCIDENCE, F-PHI 

.IET ) TE INCIDENCE, F-THETA 

.IET — > TE INCIDENCE, F-PHI 

PLEASE ENTER IHE FILE NANE OF IHE OUTPUT DATA FILE. 



FIG. A-6. INPUT NAME FOR 
OUTPUT GRAPHICS DATA FILES 
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THE GRAPHICS CAPTION IS A PERSONALIZED CAPTION ALLOWING THE 
USER TO IDENTIFY IHIS SET OF GRAPHS FROM ALL OTHER SETS. 

THE MAXIMUM LENGTH OF IHIS CAPTION IS 64 CHARACTERS. 

NOTE: WHEN USED WITH THE GRAPHICS PACKAGE, THE PROGRAN IS 
ABLE TO DIFFERENTIATE BETWEEN UPPER CASE AND LOWER CASE 
CHARACTERS. 

PLEASE ENTER ANY GRAPHICS CAPTION YOU DESIRE 



FIG. A-7. GRAPHICS CAPTION 
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other sets. The graphics caption appears on the first line 
of the header of each graph produced by MCADGRAF. 

* Mesh Density (Figure A-8) 

As discussed in Chapter II, EMCAD utilizes a semi-annular 
conformal finite-element mesh in order to conduct scattering 
calculations. This menu prompts the user for the minimum 
and maximum mesh density parameters. At azimuthal modes M = 
0 and 1, D max is used as the density parameter. The density 
decreases linearly to D m ^ n at the final mode step. Notice 
that these values are REAL. 

* Incidence Angles (Figure A-9) 

EMCAD allows up to five incident fields to impinge on the 
object of interest. This menu asks for the number of 

incidence angles, and then the value for each in degrees. 

* Number of Scattering Field Theta Points 

(Figure A-10) 

EMCAD conducts scattering calculations by starting at 9 = 0, 
and progressing clockwise to 9 = 180 in equal increments. 
This menu allows the user to determine the size of the 
increments, and indirectly the accuracy. 
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DNIN AND DNAX ARE PARAMETERS OF MESH DEMSITV IN TERNS OF 
ELEHENIS/ INTERIOR LAMBDA . INPUT UALUES ARE EXPECTED TO BE 
REAL, I.E. THE DECIMAL POINT MUST BE INCLUDED. 

TVPICAL UALUES ARE DNIN = 10. AND DNAX = 15. 

PLEASE INPUT DNIN 

10. 



PLEASE INPUT DNAX 
15. 



FIG. A-8. MESH DENSITY 



THE NUMBER OF INCIDENT FIELD ANGLES IS THE TOTAL NUMBER OF 
INCIDENT FIELDS THAT IMPINGE ON THE OBJECT OF INTEREST. 
THIS PROGRAM ALLOWS A MAXIMUM OF FIVE INCIDENT FIELD 
ANGLES. WHEN ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE A 
DECIMAL BECAUSE THE INPUT MUST BE IN INTEGER FORMAT, I.E. 



NA = 3 , NA = 1 



PLEASE INPUT THE NUMBER OF INCIDENT FIELD ANGLES. 
2 



ENTER INC FLD ANGLE (DEG) FOR 8 1 



0 . 



ENTER INC FLD ANGLE (DEG) FOR 8 2 
90. 



FIG. A-9. INCIDENCE ANGLES 
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THE NUMBER OF SCATTERING FIELD THETA POINTS DETERMINES THE 
SPACING BEIUEEN THETA POINTS DURING ENSCAI ITERATIONS AND 
CALCULATIONS. 

DELTA THETA = 180 / (NUMBER THETA POINTS - 1) SO... 

NUMBER THETA POINTS - 37 > DELTA THETA = 5 DEGREES 

NUMBER THETA POINTS = 19 > DELTA THETA = 10 DEGREES 

UHEN ENTERING TOUR ANSUER, PLEASE DO NOT INCLUDE A DECIMAL 
BECAUSE THE INPUT MUST BE IN INTEGER FORMAT. I.E. NT = 19 



PLEASE INPUT THE NUMBER OF SCATTERING FIED THETA POINTS 



FIG. A-10. NUMBER OF 
SCATTERING FIELD THETA POINTS 



so 



* Phi Angles 



(Figure A-ll) 



EMCAD allows up to eight phi angles which define the 
location from which the user views the scattering for all of 
the incident angles. This menu asks for the number of phi 
angles, and then the value for each in degrees. 

* Expansion Modes (Figure A-12) 

EMCAD partitions the incident and scattered fields into 
azimuthal modes, indexed by "M" . The solutions for these 
modes are accumulated during the "M-loop" . If 0 degrees and 
180 degrees are the only incidence angles of interest, then 
the M = 1 mode is the only mode that is needed. The values 
of N1 and N2 represent the number of spherical harmonic 
modes in the spherical core and the external region. 

E . PROGRAM OUTPUT 

As previously mentioned, EMCAD produces text output data 
and output data for graphics. The graphics data is 
discussed in Appendix H, MCADGRAF User's Manual, and the 
text output data is demonstrated in Appendix B, Portion of 
an EMCAD Sample Output. 
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THE NUMBER OF PHI ANGLES IS THE TOTAL NUMBER OF PHI ANGLES. 
THIS PROGRAM ALLOWS A MAXIMUM OF EIGHT PHI ANGLES. WHEN 
ENTERING HOUR ANSWER PLEASE DO NO! INCLUDE A DECIMAL 
BECAUSE THE INPUT MUST BE IN INTEGER FORMAT. I.E. 



HP - 3 , HP = 1 



PLEASE INPUT THE NUMBER OF PHI ANGLES. 
2 



ENTER PHI FLD ANGLE (DEG) FOR I 1 



90. 



ENTER PHI FLD ANGLE (DEG) FOR I 2 
180. 



FIG. A— 11. PHI ANGLES 
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ENTER MSTOP C.LE. 13) 



ESTIMATED "MINIMUM" VALUE IS: 1 

1 

ENTERING NO. OF INTERNAL AND EXTERNAL EXPANSION MODES 
ESTIMATED (KI*RMIN) "MINIMUM" INTERNAL N1 IS: 12 

ENTER N1 (.GE. MSTOP) 

12 

ESTIMATED (UNWIN) "MINIMUM" INTERNAL N2 IS: 21 

ENTER N2 (.GE. MSTOP) 

21 

FIG. A-12. EXPANSION MODES 
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APPENDIX B 



PORTION OF AN EMCAD SAMPLE OUTPUT 

This Appendix demonstrates a portion of an EMCAD 
sample output for the specific case chosen where 
M indexes from 0 to 7 inclusive. Only the M = 1 
case is fully shown, as well as the accumulated 
solution for all eight modes. 



xxxxxxxxxxxxxxxxxxxx EMCAD OUTPUT DATA xxxxxxxxxxxxxxxxxxxxxxxxxxx 
COMPLEX Er(n) AND UrCn) 

1 C 2.610E+00 -1.000E-06) C 1.000E+0Q -1.000E-06) 

INCIDENT FIELD ANGLES 
N ALPHACN) 

1 0 . DEG 

2 90. DEG 

3 180. DEG 

SERIES SOLUTION PARAMETERS 
MSTART = 0 

MST0P= 7 

Nl= 8 

N2= 15 
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1 ***x*x***x*x*x 



MESH PARAMETERS 
RADII (WAVELENGTHS) 



RADMIN S 


. 64000E+00 


RADMAX= 


. 96000E+00 


RADA = 


. 55241 E+QO 


RADB = 


. 10476E+01 


INTEGER 


CONSTANTS 


NDS = 


13 


IMX = 


59 


NR1 = 


1 


NR2 = 


1 



MESH DENSITY AND ELEMENT SIDE RATIOS 
DENSITY/LAMBDA-OUT = 24.233 

DENSITY/LAMBDA-IN = 15.000 



DSA/DR = 


= .725 


DSB/DR : 


s 1.375 
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RICCATI 


MARCHING 


1 


. 6 03E+0 1 


1 


. 6 03E+0 1 


2 


. 6 02E+01 


3 


. 6 02E+01 


4 


. 6 00E+01 


5 


. 599E+ 01 


6 


. 597E+01 


7 


. 595E+01 


8 


c 592E+01 


9 


. 539E+01 


10 


.5S6E+01 


11 


. 533E+01 


12 


.579E+01 


13 


.575E+01 


14 


.571E+01 


15 


. 566 E+ 0 1 


16 


.562E+01 


17 


. 557 E+01 


18 


. 552E+0 1 


19 


. 547 E+0 1 


20 


. 542E+ 01 


21 


. 537 E+0 1 


22 


. 531 E+0 1 


23 


. 526 E+0 1 


24 


. 52 0E+0 1 


25 


. 515E+01 


26 


. 509E+01 


27 


. 503E+01 


23 


. 498 E+0 1 


29 


. 492E+0 1 


30 


.487 E+ , 0 1 


31 


. 432E+01 


32 


. 477 E+0 1 


33 


. 47 1 E+0 1 


34 


. 466 E+0 1 


35 


. 46 1E+01 


36 


. 457 E+01 


37 


« 452E+0 1 


33 


. 443 E + 01 


39 


.443E+01 


40 


. 439E+01 


41 


. 435E+ 0 1 


42 


. 432E+0 1 


43 


. 428 E+0 1 


44 


. 425 E+01 


45 


. 422 E+0 1 


46 


.419 E+01 


47 


.416 E + 01 


43 


. 414E+01 


49 


. 412E+01 


50 


. 41 0 E+0 1 


51 


. 403 E+0 1 


52 


. 406 E+0 1 


53 


. 405E+01 


54 


. 404E+01 


55 


.403E+01 


56 


. 403E+01 


57 


. 402E+0 1 


53 


. 402E + 01 


59 


. 402E+01 
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INTERIOR AND SCATTERED FIELD 
AZIMUTHAL MODE NUMBER M= 1 

INCIDENCE ANGLE ALPHA = 0. 



EXPANSION COEFFICIENTS 
DEGREES 



xxxxxxx INTERIOR COEFFICIENTS xxxxxxx 



STAGGERED 


COEFFICIENTS: 


TM 


INCIDENCE 






N 


A 


CN) 






BCN) 


1 


. 500 E+00 




. 266 E+0 0 


- . 318E+00 


. 298 E+0 0 


2 


. 209E+00 


- 


. 773E-01 


. 233E+0 0 


. 192E+00 


3 


- . 1 03E+G0 


- 


. 1 77 E+0 0 


. 172E+00 


- . 66 0 E-0 1 


A 


- . 207 E+0 0 




. 582E-0 1 


.163E-01 


-. 1A7E+00 


5 


. 996 E-02 




. 169E-01 


- . 1 35E+0 0 


. 98 1 E-02 


6 


.886E-01 




. 390E-01 


-, 109E+00 


. 599E-01 


7 


.61AE-01 




. 1 38 E“0 1 


-.373E-01 


. A37E-01 


8 


.197E-01 




. 325E-Q2 


- . 6 25E-02 


. 136E-01 


STAGGERED 


COEFFICIENTS: 


TE 


INCIDENCE 






N 


A 


(N) 






BCN) 


1 


.266E+00 


- 


. 500 E+0 0 


. 298E+0 0 


. 318 E+0 0 


2 


- . 778E-01 


- 


. 209 E+0 0 


. 192E+00 


- . 233E+0 0 


3 


- . 177 E+0 0 




. 103E+00 


-.660E-01 


-. 172E+00 


A 


.582E-01 




. 207 E+0 0 


1A7E+00 


-.163E-01 


5 


. 169E-01 


- 


. 996 E-02 


. 98 1 E-02 


. 135E+00 


6 


. 39 OE-0 1 


- 


.886E-01 


.599E-01 


. 109E+00 


7 


. 138E-01 


- 


.61AE-01 


. A37E-01 


.373E-01 


8 


- . 325E-02 




. 197E-01 


.136E-01 


. 6 25E-02 


xxxxxxx SCATTERED 


COEFFICIENTS 


xxxxxxx 




STAGGERED 


COEFFICIENTS-. 


TM 


INCIDENCE 






N 


C(N) 






DC N ) 


1 


- . 9 02E-0 1 


- 


. 1AAE-01 


. 10AE+00 


. 22 A E+0 0 


2 


-.115E-01 


- 


. 195E+00 


. 955E-02 


-. 123E+00 


3 


- . 1 03E+0 0 




. 157 E+00 


99AE-01 


- . 850 E-0 1 


A 


. 138 E+0 0 




. 131E+00 


- . 928 E-0 1 


. 6 96 E-0 1 


5 


. A72E-01 


- 


.95AE-01 - 


.670E-01 


. 158E+00 


6 


.168E-01 


- 


.616E-01 


.510 E-0 1 


. 6 59 E-0 1 


7 


. 36 1 E-02 


- 


. 171E-01 


. 1 59 E-0 1 


.135E-01 


8 


. A66E-03 


- 


. 306 E-02 


. 298 E-02 


. 160E-02 


9 


. 372E-0 A 


- 


. A16E-03 


.398E-03 


. 103E-03 


10 


. 192E-05 


- 


. A97E-0A 


. A51E-0A 


- . 239 E-05 


11 


-.151E-06 


- 


.692E-05 


. 6 36 E-05 


-.951E-06 


12 


-.606 E-07 


- 


. 119E-05 


.111E-05 


-. A22E-06 


13 


-.3A9E-07 


- 


.207E-06 


. 255E-06 


-.117E-06 


1A 


- . 150E-07 




. A33E-07 


. 523 E-07 


-.3A2E-07 


15 


-.336E-08 


— 


.776E-08 


. 112E-07 


- .626E-08 


STAGGERED 


COEFFICIENTS: 


TE 


INCIDENCE 






N 


C(N) 






DC N ) 


1 


-.1AAE-01 




.902E-01 


. 22AE+00 


10AE+00 


2 


- . 1 9 5E+00 




. 115E-01 


-. 123E+00 


-.955E-02 


3 


. 157 E+0 0 




. 103E+00 


- . 850 E-0 1 


. 9 9 A E- 0 1 


A 


. 131E+00 


- 


. 138E+00 


. 6 96 E-0 1 


. 928 E-0 1 


5 


-.95AE-01 


- 


. A72E-01 


. 158 E+0 0 


-.670 E-0 1 


6 


-.616E-01 


- 


. 168E-01 


. 6 59 E-0 1 


- . 51 0 E-0 1 
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7 


171E-01 


-.361E-02 


. 135E-01 


159E-01 


3 


306E-02 


-.466E-03 


. 160E-02 


-.293E-02 


9 


-.416E-03 


- . 37 2E~ 04 


. 103E-03 


398E-03 


10 


-.497E-04 


192E-05 


239E-05 


- . 451E-04 


11 


- . 6 9 2E-0 5 


. 151E-Q6 


- . 951 E-Q6 


- .636E-05 


12 


119E-05 


.606E-07 


-.422E-06 


111E-G5 


13 


- . 207E-06 


. 349E-07 


- , 117E-06 


-.255E-06 


1$ 


933E-07 


. 150E-07 


-.342E-07 


- . 528 E- 07 


15 


“ . 776 E-08 


. 336E-08 


-.626E-03 


112E-07 
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INTERIOR AND SCATTERED FIELD EXPANSION COEFFICIENTS 
AZIMUTHAL MODE NUMBER M= 1 

INCIDENCE ANGLE ALPHA = 90. DEGREES 



xxxxxxx INTERIOR COEFFICIENTS xxxxxxx 



STAGGERED 


COEFFICIENTS: 


TM 


INCIDENCE 






N 


A 


CN) 






BCN) 


1 


. 457 E+0 0 




.212E+00 


. 555E-02 


. 8 3 0 E-0 1 


2 


.913E-02 


- 


. 474E-01 


- . 246 E+0 0 


-. 152E+00 


3 


.354E-01 




.375E-01 


-.262E-01 


. 429 E-0 1 


4 


.260E-03 


- 


. 239 E~0 2 


- . 588 E-0 1 


- . 7 03E-01 


5 


. 650E-02 




.110E-01 


-.653E-01 


-.473E-01 


6 


. 782E-02 




.751E-02 


. 204E-02 


-. 190E-01 


7 


.198E-02 




, 29 0 E-02 


. 651 E-0 2 


- . 6 26 E-0 2 


8 


.708E-03 




. 1 Q4E-G2 


.371E-02 


- .277E-02 


STAGGERED 


COEFFICIENTS: 


TE 


INCIDENCE 






N 


A 


CN) 






BCN) 


1 


.662E-01 


- 


.344E-01 


. 624E+00 


. 267 E-0 1 


2 


. 294E+0 0 




. 353 E-0 1 


. 397E-01 


-. 379E-01 


3 


. 833E-0 1 


- 


. 504E-01 


. 342E-0 1 


. 253 E-0 1 


A 


.445E-01 




. 346 E-0 1 


. 21 9 E-0 1 


-. 115E-01 


5 


.653E-01 




. 237 E-0 1 


-. 173E-02 


. 1 19 E-0 1 


6 


-.126E-01 




. 127E-01 


- « 56 1 E-03 


.916E-02 


7 


131E-01 




. 504E-02 


- . 8 96 E-03 


.3 38E-02 


8 


- . 574E-02 




. 199E-02 


-.384E-03 


. 1 0 0E-02 


xxxxxxx SCATTERED 


COEFFICIENTS 


xxxxxxx 




STAGGERED 


COEFFICIENTS : 


TM 


INCIDENCE 






N 


C(N) 






DC N ) 


1 


.901E-01 


- 


. 229 E+0 0 


. 536 E-0 1 


-.517E-01 


2 


--376E-01 




. 20 0 E-0 1 


105E+00 


. 169E+00 


3 


.411E-01 


- 


. 318 E-0 1 


.40 0 E-0 1 


-.698E-02 


A 


-.387E-03 


- 


.606E-03 


- . 847 E-0 1 


.700 E-0 1 


5 


.127E-01 


- 


.868E-02 


342E-01 


. 348 E-0 1 


6 


.497E-02 


- 


.408E-02 


- . 895E-02 


-.529E-02 


7 


.103E-02 


- 


. 8 06 E-0 3 


- . 1 95E-0 2 


- . 393E-02 


8 


. 16 4E-03 


- 


. 112E-03 


- . 297 E-03 


- . 1 05E-02 


9 


. 124E-04 


- 


. 197E-05 


- . 402E-0 4 


181E-03 


10 


. 339E-05 




.149E-06 


- . 346 E-0 5 


- . 26 5 E- 0 4 


11 


- . 352E-06 




.528E-06 


- . 403E-O6 


- . 292E-05 


12 


. 186 E-Q6 




. 523E-07 


- .977E-07 


-.256E-06 


13 


- . 314E-07 




.132E-07 


. 357 E‘ 08 


-. 189E-07 


14 


. 795E-08 




.301E-08 


- . 6 0 2E- 0 8 


-.652E-09 


15 


-.137E-08 


“ 


. 282E-09 


. 7 30 E-0 9 


-. 169E-09 


STAGGERED 


COEFFICIENTS: 


TE 


INCIDENCE 






N 


CCN) 






DC N ) 


1 


-.758E-02 




. 29 1 E-0 1 


. 722E-02 


. 198E-03 


2 


.260E-01 


- 


. 267 E-0 1 


- . 1 52E-0 1 


.294E-01 


3 


-.582E-01 




.473E-01 


.212E-01 


- . 237 E-0 1 


4 


.726E-01 


- 


.648E-01 


-. 108E-01 


-.879E-02 


5 


.301E-01 


- 


. 344 E-0 1 


. 157E-01 


. 513E-02 


6 


. 768E-02 




.348E-02 


.66 7 E-0 2 


.194E-02 
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7 


. 162E-02 


. 324E-02 


. 143E-02 


.572E-03 


8 


.236E-03 


.848E-03 


.222E-03 


. 114E-03 


9 


. 330E-04 


. 135E-03 


. 187E-04 


. 26 1 E- 0 4 


10 


. 241E-05 


. 186E-04 


. 394E-05 


. 294E-05 


11 


. 50 1 E- 06 


. 171E-05 


-.363E-06 


.332E-06 


12 


.533E-07 


. 955E-07 


. 176E-06 


. 7 03 E- 0 7 


13 


- . 925E-03 


- . 222E- 08 


374E-07 


. 113E-07 


14 


, 457E-0S 


206E-08 


.757E-08 


176E-03 


15 


- . 86 0 E-09 


-.157E-09 


154E-08 


528E-09 
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ACCUMULATED FORWARD AND BACK SCATTERED AMPLITUDES 
M = 1 



STAGGERED AMPS: INCIDENT TM 



ALPHA 

(DEG) 




FWD 




BK 




0 . 


- . 8 05E+0 1 


-.177E+02 


139E+01 




-.550E+01 


90 . 


- . 29 AE+01 


- . 359 E+ 0 1 


. 200E+01 




. 1A1E+01 


180 . 


- . 77 6 E+0 1 


- . 172E+02 


- . 371E+01 




. A74E+01 


STAGGERED 


AMPS: INCIDENT TE 








ALPHA 

(DEG) 




FWD 




BK 




0 . 


-.805E401 


- . 177 E+ 02 


139E+01 




-.550E+G1 


90 . 


- . 213E+01 


~ . 295E+0 1 


. I96E-01 




- . 2A7 E+ 0 0 


180. 


- , 776E+01 


172E+02 


371E+01 




. A7AE+01 



SCATTERING CROSS SECTION CONVERGENCE 



M= 1 



NINC = 1 



N 


NEW TERM 


PCT CHANGE 


1 


.6E-01 


100.0000 


2 


. 2E+Q 0 


80 .6270 


3 


. 7 E+ 0 0 


69 . 3592 


A 


. 1E+01 


58 . 7 A50 


5 


. 2E+0 1 


A 7 . 0760 


6 


. 1E+01 


17.3836 


7 


. 1E+00 


1.7682 


8 


. AE-02 


. 0730 


9 


. 9 E-0 A 


.0017 


10 


. 2E-05 


.0000 


11 


. AE“07 


.0000 


12 


. 2E-08 


.0000 


13 


. 1 E-09 


. 0000 


1 A 


. 6 E-l 1 


.0000 


15 


.3E-12 


.0000 



NINC 



N 


NEW TERM 


PCT CHANGE 


1 


. 6 E-0 1 


100 . 0000 


2 


. 2E+0 0 


77 .2227 


3 


. 6 E-0 1 


18 . 8 1 0 A 


A 


.3E+00 


52, 99 3 A 


5 


. 1 E+G 0 


17 . A57 1 


6 


. IE-01 


1 .6261 


7 


. 3E-02 


. 3A95 


8 


.2E-03 


. 0301 


9 


. 9E-Q5 


.0012 


10 


.3E-06 


. 0000 


11 


. AE-08 


.0000 


12 


.7E-10 


. 0000 


13 


.IE-11 


. 0000 


1A 


.IE-12 


.0000 


15 


.3E-1A 


. 0000 
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NINC = 


3 


4 




N 




NEW TERM 


PCT CHANGE 


1 




. 5E-01 


1 00 .0000 


o 

c 




. 3E+00 


85 .8396 


3 




. 1E+01 


76 . 4532 


4 




. 3E+01 


63.6074 


5 




. 1E+01 


22 9305 


6 • 




. 2E+00 


4 . 1769 


7 




. IE-01 


.1771 


8 




. 3E-03 


.0063 


9 




.IE-04 


, 00 02 


10 




.4E-06 


.0000 


11 




.9E-08 


. 0000 


12 




.4E-09 


.0000 


13 




.2E-10 


.0000 


14 




.IE-11 


. 0000 


15 




.5E-13 


. 0000 


NINC = 
N 


4 


NEW TERM 


PCT CHANGE 


1 




.6E-01 


100 .0000 


2 




. 2E+00 


80.6270 


3 




. 7 E+00 


69 . 3592 


4 




. 1E+01 


58.7450 


5 




. 2E+01 


47 . 0760 


6 




. 1E+01 


17.3836 


7 




. 1 E+00 


1.7682 


8 




.4E-02 


. 0730 


9 




. 9 E~04 


.0017 


10 




. 2E“05 


.0000 


11 




.4E-07 


. 0000 


12 




. 2E- 08 


.0000 


13 




. IE-09 


. 0000 


14 




.6E-11 


. 0000 


15 




. 3E-12 


.0000 



NINC = 5 



N 


NEW TERM 


PCT CHAN' 


1 


.8E-03 


100 . 0000 


2 


. IE-01 


93.3376 


3 


. 9E-01 


87 .6869 


4 


. 3E+00 


73.3983 


5 


. 1 E+00 


24.8194 


6 


.IE-01 


2.0363 


7 


.2E-02 


.4062 


8 


.2E-03 


. 0320 


9 


.6E-05 


.0011 


10 


. IE-06 


0000 


11 


. 2E-08 


.0000 


12 


. 3E~1 0 


. 0000 


13 


. IE-11 


.0000 


14 


.8E-13 


. 0000 


15 


. 4E-14 


. 0000 
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MODAL CROSS SECTION EVOLUTION 
NORMALIZED W.R.T. WAVELENGTH SQUARED 



INCIDENCE ANGLE = 0. DEGREES 

INCIDENT TM POLARIZATION 



INCIDENT TE POLARIZATION 



M 


SCT 


ABS 


TOT 


SCT 


ABS 


TOT 


0 


. 000E+00 


. 000E+00 


. 0 0 0 E+ 0 0 


. 000E+00 


. 000 E+0 0 


. 00 0 E+0 0 


1 


. 558E+01 


.705E-01 


. 565E+01 


. 558 E+0 1 


.705E-01 


. 565E+01 


2 


. 553E+Q1 


. 7 05E-0 1 


.565E+01 


. 558E+01 


.705E-01 


.565E+01 


3 


. 558E+01 


. 7 05E-01 


.565E+01 


. 5 58 E+0 1 


.705E-01 


.565E+01 


4 


. 558E+01 


.705E-01 


. 565E+01 


. 558E+0 1 


. 7 05E-01 


.565E+0I 


5 


.55SE+01 


. 7 05E-01 


.565E+01 


. 558E+0 1 


. 7 05E-01 


. 56 5E+0 1 


6 


. 558E+01 


. 7 05E-01 


. 565E+0 1 


. 558E+ 0 1 


. 7 05E-01 


.565E+0I 


7 


. 558E+01 


. 705E-01 


.565E+01 


. 558 E+0 1 


. 7 05E-01 


. 56 5 E+0 1 



INCIDENCE ANGLE = 90. DEGREES 

INCIDENT TM POLARIZATION 



INCIDENT TE POLARIZATION 



M 


SCT 


ABS 


TOT 


SCT 


ABS 


TOT 


0 


.247E+00 


. 1 00E+Q0 


. 347 E+0 0 


. 326 E x Q0 


. 105E+00 


. 431 E+0 0 


1 


. 104 E+0 1 


.998E-01 


. 1 14E+01 


. 834E+00 


. 1 05 E+0 0 


. 940E+00 


2 


.197E+01 


.110 E+0 0 


. 208E+01 


. 202E+01 


. 1 1 0E+00 


.213E+01 


3 


. 317 E+0 1 


.106E+00 


. 328 E+0 1 


. 285E+0 1 


. 106E+00 


. 296 E+ 0 ] 


4 


. 440 E+0 1 


. 113E+00 


. 451 E+0 1 


.396E+G1 


. 1 18 E+0 0 


. 407 E+0 1 


5 


. 521 E+0 1 


. 1 44E+0 0 


. 536 E+ 0 1 


. 519E+01 


. 199 E+0 0 


. 539E+01 


6 


. 525E+0 1 


. 166 E+0 0 


. 542E+0 1 


. 523E+0 1 


. 216 E+0 0 


. 545E+01 


7 


. 525E+01 


. 166E+00 


. 542E+0 1 


. 523E+0 1 


. 216 E+0 0 


. 54 5 E+0 1 
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INCIDENCE ANGLE 



18 0 . DEGREES 





INCIDENT 


TM POLARIZATION 


INCIDENT 


TE POLARIZATION 


M 


SCT 


AB$ 


TOT 


SCT 


ABS 


TOT 


0 


.220E-12 


- . 220E-1 2 


. 000E+00 


.209E-12 


-.209E-12 


. 000E+C 


1 


. 546 E+0 1 


. 255E-0 1 


. 548E + 0 1 


. 546 E+01 


. 255E-01 


. 548 E+C 


2 


. 546 E+01 


.255E-01 


. 548 E+01 


. 546 E+01 


. 255E-01 


. 548E+C 


3 


. 546 E+0 1 


. 255E-0 1 


. 548E+01 


. 546 E+01 


. 255E-01 


. 548E+C 


4 


. 546 E+01 


. 255E-0 1 


. 548 E+01 


. 546 E+01 


. 255E-0 1 


. 548 E+G 


5 


. 546 E+01 


.255E-01 


. 548 E+01 


. 546 E+01 


. 255E-0 1 


. 548E+G 


6 


. 546 E+01 


. 255E-01 


. 548E+01 


. 546 E+01 


.255E-01 


. 548 E+0 


7 


. 546 E+01 


. 255E-01 


. 548 E+01 


. 546 E+01 


. 255E-0 1 


. 548 E+0 
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***** SCATTERED FIELD COMPLEX AMPLITUDES ***** 



INCIDENCE ANGLE = 0. DEGREES 

PHI = 90. DEGREES 

STAGGERED RAD FLDS : TM INCIDENCE 



THETA F-THETA F-PHI 





MAG 


DEG 


MAG 


DEG 


0 . 


.851E-06 


.00 


. 195E+02 


65.60 


5. 


. 821 E- 06 


. 00 


. 137E+02 


65.19 


10. 


.735E-06 


.00 


. 166E+02 


63.93 


15 . 


.6 03E-06 


.00 


. 133E+Q2 


61 71 


20 


. 96 1 E” 06 


.00 


. 925E+0 1 


58 . 17 


25. 


. 322E-06 


. 00 


. 999E+01 


51.69 


30 . 


. 229 E- 06 


.00 


. 119E+01 


19 . 66 


35. 


.207E-06 


. 00 


. 253E+01 


-112.31 


90 . 


. 225 E- 06 


. 00 


.9S9E+01 


-129.99 


95 . 


. 239 E- 06 


. 00 


. 6 03E+0 1 


-131 .99 


50 . 


.235E-06 


.00 


.608E+01 


-138 . 30 


55 . 


.217E-06 


.00 


. 515E+01 


-199.73 


60 . 


. 1 97 E" 06 


.00 


.353E+01 


-151 .82 


65. 


. 1 8 3 E-06 


.00 


.159E+01 


-161 .86 


70. 


. 173E-06 


.00 


. 382E+00 


99.96 


7 5 . 


. 1 6 0 E-06 


.00 


.137E+01 


18.57 


80 . 


. 138E-06 


. 00 


.281E+01 


10.27 


85 . 


.111E-06 


. 00 


.307E+01 


3.63 


90 . 


. 933E-07 


.00 


. 270E+01 


-2.13 


95. 


.993E-07 


.00 


. 186E+01 


-6 . 50 


100 . 


. 12GE-06 


.00 


.812E+00 


-5.21 


105. 


. 137E-06 


.00 


. 290E+00 


118.65 


110. 


. 138E-06 


. 00 


. 996E+00 


190.36 


115. 


. 121 E- 06 


. 00 


. 133E+01 


138 99 


120. 


.3S9E-07 


.00 


. 121E+01 


133.29 


125. 


.557E-07 


. 00 


.716E+00 


123.28 


130 . 


.533E-07 


.00 


. 189E+00 


39.52 


135. 


.939E-07 


. 00 


. 7 59E + 0 0 


-30.93 


190. 


.129E-06 


. 00 


. 125E+01 


-35.32 


195. 


. 1 53E-Q6 


. 00 


.192E+01 


-29.65 


150. 


.165E-06 


.00 


.139E+01 


-9 .95 


155. 


.173E-06 


.00 


. 199E+01 


26.39 


16 0 


.189E-06 


.00 


.231E+01 


55.01 


165. 


. 202E-06 


.00 


.398E+01 


67.53 


170. 


. 229E-06 


. 00 


. 96 0 E+ 0 1 


72.91 


175. 


.291E-06 


.00 


. 539E+01 


75.20 


130. 


.298E-06 


.00 


. 568E+01 


75.85 
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***** SCATTERED FIELD COMPLEX AMPLITUDES ***** 



INCIDENCE ANGLE = 0. DEGREES 



PHI = 90. DEGREES 



STAGGERED RAD FLDS : TE INCIDENCE 



THETA F-THETA 



F-PHI 





MAG 


DEG 


MAG 


0. 


. 195E+02 


-11A. AO 


.851E-06 


5. 


. 183E+02 


-115.57 


.819E-06 


10 . 


. 163E+02 


-119.26 


.724E-06 


15. 


. 139E+02 


-126 .13 


.58QE-06 


20 . 


. 105E+02 


-137.69 


. 4 04E-06 


25. 


. 7 38 E+ 0 1 


-157 .34 


. 218E-06 


30 . 


. 52AE+0 1 


169.80 


.497E-G7 


35. 


. A73E+01 


128.60 


. 1 1 0E-06 


AO . 


.515E+01 


95.19 


. 21 2E-06 


A 5 . 


.547E+01 


69.76 


.264E-06 


50 . 


.537E+01 


A6 .52 


. 266 E- 06 


55. 


. A97E+01 


21 .84 


. 225E-06 


60 . 


. A51E+01 


-5.65 


. 1 54 E- 06 


65, 


. A19E+01 


-34.83 


. 6 93E-07 


70 . 


.396E+01 


-63.23 


.167E-07 


75. 


.366E+01 


-90.06 


, 8 18E-07 


80. 


.316E+01 


-117.25 


. 1 23E-06 


85. 


. 255E+0 1 


-149.31 


.134E-06 


90 . 


.21AE+01 


168.63 


. 118E-06 


95. 


.227E+01 


123,71 


.814E-07 


100 . 


.275E+01 


89 .95 


. 355E-07 


105. 


.31AE+01 


65.78 


. 127E-07 


110. 


.317E+01 


45.98 


.435E-07 


115. 


. 277 E+0 1 


26.74 


.579E-07 


120 . 


. 203E+01 


3 . 03 


.527E-07 


125. 


. 127E+01 


-39.54 


. 3 1 3E-07 


130 . 


. 133E+01 


-107.69 


.828E-08 


135. 


.215E+01 


-147 . 00 


.329E-07 


1 AO , 


. 29 5E+ 0 1 


-169.25 


. 545E-07 


1A5. 


.3A9E+01 


172.65 


. 6 19E-07 


150 . 


. 378 E+ 0 1 


154.61 


. 585E-07 


155. 


. 396E+01 


135.25 


.653E-07 


160 . 


. A21E+01 


115.51 


.101E-06 


165. 


. A62E+01 


98 .19 


. 152E-06 


170 . 


.512E+01 


85.58 


. 201 E“ 06 


175. 


. 552E+0 1 


78.23 


.236E-06 


180 . 


. 568E+01 


75.85 


.248E-06 



DEG 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
,00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
.00 
. 00 
.00 
.00 
.00 
.00 
.00 
. 00 
.00 
.00 
. 00 
. 00 
00 
. 00 
. 00 
. 0 0 
00 
.00 
.00 
.00 
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ooononooooonooonoooooooooooooonoooonoooonooo 



APPENDIX C 



EMCAD SOURCE CODE 



PROGRAM EMCAD 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THIS PROGRAM COMPUTES EM FIELD SCATTERING FROM LAYERED 
BODIES OF REVOLUTION USING THE COUPLED AZIMUTHAL POTENTIAL 
(CAP) FORMULATION IN CONJUNCTION WITH A VARIATIONAL FINITE- 
ELEMENT TECHNIQUE AND A TRI -REGIONAL UNIMOMENT METHOD. THE 
NUMERICAL SOLUTION IS PERFORMED USING A TWO-SWEEP RICCATI 
TRANSFORM WITH CIRCUMFERENTIAL MARCHING. UNIT 1 IS USED FOR 
SEQUENTIAL UNFORMATTED STORAGE OF GENERATED RICCATI ARRAYS TO 
DISK o UPDATED VERSION OF CDC-7600 U.C. BERKELEY CODE ,T EMSCAT ,f 
WITH LAYERED MEDIA ADDITIONS. 

MODS BY M . A . MORGAN MAR 87-FEB 88 
E.M. CONNOLLY MAR 87-FEB 88. 

xxxxxxxxxxxxxxxxxxxx EMCAD PARAMETERS xxxxxxxxxxxxxxxxxxxxxxxxx 

xx NA - NO. OF INCIDENT ANGLES/ DALPHACJ) (DEGREES) 

XX DMIN/DMAX - MIN AND MAX MESH DENSITY ( EL EM/ LAMB DA ) 

XX NR1,NR2 - RADIAL NODE SPACING BETWEEN RA-R1 8 R2-RB 

xx NLAY - NO. OF LAYERED SURFACES, NLAY+1 IS FREE SPACE 

xx ER(NLAY) ,UR(NLAY) - COMPLEX ER AND UR OF PENETRABLE LAYERS 

XX MSTART/MSTOP - AZIMUTHAL MODE RANGE ON ff M ,f 

xx N1,N2 - UPPER LIMITS ON INTERIOR & EXTERIOR MODES "N" 

xx NT , NP - NO. OF FIELD PRINTOUT THETA 8 PHI POINTS 

xx NB, NBMX - NO. OF INTERIOR NODAL POTENTIALS/THETA STEP 

xx NCMAX - NO OF INTERIOR + EXTERIOR EXPANSION COEF'S 

xx NMMAX - NCMAX + 2XNA (NO. UNIMOMENT BC f S) 

xx N I NC , NINMX - 2XNA (2 POLARIZATIONS/INC ANGLE) 

XX BC - 0 OR 1 : Z-AXIS HOMO DIRICHLET OR CAUCHY BC f S 
xx I MX, I BIG - NO. OF THETA STEPS IN MARCHING 8 DATA FILE 

xxxxxxxxxxxxxxxxxxxx ARRAY DIMENSIONS xxxxxxxxxxxxxxxxxxxxxxxxxx 

XX R , RI NV - (NBMX, NBMX) RICCATI MATICES 

xx S, SI NV - (NBMX, NMMAX ) S-VECTOR ARRAYS 

xx FMOM/DUM - ( NCMAX, NCMAX ) COEFFICIENT MOMENT ARRAYS 

xx DMOM , COEF - (NCMAX, NINMX) DRIVING 8 COEF ARRAYS (SAME) 

xx FT1 TO FP2 - (NT/NP,NA) SCATTERED E-FIELDS 

xx RAD - (2,5) RADII(THETA) OF LAYERS AT I 8 1+1 STEPS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CHARACTERXl BEL L , FF , L EVEL 
CHARACTERXS GRFDAT , PRTDAT 

CHARACTERXl 2 GRFDT 1 , GRFDT2 , GRFDT3 , GRFDT4 , PRTDT1 
CHARACTERX64 DUMMY , GRA FLAB , STRDAT , HDR, MATPAR 
COMPLEX A,B,G,P, E1,H1, E2,H2, EB1,HB1, EB2,HB2 
COMPLEX F30, F21 , F12 , F20, Fll , F10, FI, F2, FT1, FP1, FT 2, FP2 
COMPLEX R(58,58),S(58,80),RINV(58,58),SINV(58,80),ER(6),UR(6) 
COMPLEX FMOM(70,70),COEF(70, 10) 

INTEGER BC, ELND, IDUMMY 
REAL RAD ( 2 , 5 ) 

COMMON/ ONE/X( 62) , Y( 62 ) 

COMMON/TWO/ ELND( 6 0, 3), NDEL( 62,3), NDCT( 62) 

COMMON/THREE/ EB1 ( 3, 8 0),HB1( 3, 80), EB2 ( 3 , 80 ) , HB2( 3 , 8 0 ) 
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COMMON/ FOUR/ AC 53, 4 ),B( 53, 6 ),G( 53,4), PC 4, 8 0),KA( 58), KB (53), KG (53) 
COMMON/FIVE/E1(60,3,3),H1(60,3,3),E2(60,3,3),H2(60,3,3) 

COMMON/ S IX/ FT1( 37, 8, 5) , FP1 ( 37 , 3 , 5 ), FT2( 37 , 8 , 5 ) , FP2( 37 , 3 , 5 ) 
COMMON/ S EVEN/ SCAT (13,10) , A BSB ( 1 3 , 10) , TOTL ( 1 3 , 10) 

COMMON/ EIGHT/ ALPHA (5) ,PHI(3),DALPHA(5), DPHIC3) 

COMMON/TEN/ F3 0, F21 , F12 , F20 , FI I , FI 0 , P30 , P21 , P 12 , P20 , P 1 1 , P 1 0 
EQUIVALENCE (RINV,SINV) 

BELL =CHAR ( 7 ) 

FF = CHARC12) 

C SETTING MAXIMUM MATRIX DIMENSIONS 

NCMAX=70 
NMMAX=80 
NBMX=58 
N I NMX= 1 0 
PI = 3 . i A 159 27 
DTR = P 1/18 0 . 

C READING INC. FIELD, MESH, AND SOLUTION PARAMETERS 

WRITEC*, 137 ) 

READC*, 101) MATPAR 
0PEN(2,FILE=MATPAR) 

WRITEC*, 112) 

READC *, 101) STRDAT 
OPENC 10, FILE=STRDAT) 

READC 2,101) HDR 
WRITEC*, 101) 

WRITEC*, 101) HDR 
READC 2,101) HDR 
WRITEC *, 101 ) HDR 
WRITEC*, 25) 

READC*, 101) DUMMY 
WRITEC*, 101) FF 
WRITEC*, 10) 

READC*, 101) DUMMY 
WRITEC*, 101) FF 
WRITEC*, 12) 

READC*, 101) LEVEL 
WRITEC*, 101) FF 
IFC LEVEL . EQ . ’ 3® ) GO TO 99 
IFC LEVEL . EQ . ’2® ) THEN 
WRITEC*, 113) 

READC*, 101) PRTDAT 
ELSE 

WRITEC*, 13) 

READC*, 101) PRTDAT 
ENDIF 

PRTDT1 = PRTDAT// s .OUT 1 
0PENC3, FILE=PRTDT1) 

I F( LEVEL , EQ . ’ 2 ’ ) THEN 
WRITEC*, 114) 

READC*, 101) GRFDAT 
ELSE 

WRITEC * , 14 ) 

READC*, 101) GRFDAT 
ENDIF 

GRFDT1 = GRFDAT//’ .TMT’ 

GRFDT2 = GRFDAT//’ . TMP 1 
GRFDT3 = GRFDAT//’ .TET 1 
GRFDT4 = GRFDAT//’ .TEP’ 

OPENC 4, FILE = GRFDT1) 

0PENC7 , FILE=GRFDT2) 

OPENC 8, FILE=GRFDT3 ) 

OPENC 9 , FILE=GRFDT4 ) 

I FC LEVEL . EQ . ’2’) THEN 
WRITEC*, 115) 

READC*, 101) GRAFLAB 
ELSE 
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WRITECX, 15) 

READ(x,101) GRAFLAB 
ENDIF 

W R I T E C 4 , 101) GRAFLAB 
W R I T E C 7 f 101) GRAFLAB 
H R I T E C 3 , 101) GRAFLAB 
WRITEC9,101) GRAFLAB 
IFCLEVEL .EQ. '2' ) THEN 
WRITECX, 116 ) 

READC x , X) DMIN 
WRITECX, 117) 

READC x, x) DMAX 
ELSE 

WRITECX, 16) 

READC X , X ) DMIN 
WRITECX, 17) 

READC X , x ) DMAX 
ENDIF 

WRITEC3,101) FF 
WRITEC3, 100) 

WRITEC 3, x ) 

WRITEC3, 111) 

READC 2 f 102) NL AY 
DO 22 L=1,NLAY 
READC 2, 103) ERCL),URCL) 
WRITEC 3,107) L,ERCL),UR(L) 
22 CONTINUE 

READC 10,136) IBIG 
READC1Q,103) RMIN,RMAX 
IFCLEVEL. EQ. '2') THEN 
WRITECX, 118) 

READC x , x ) NA 
ELSE 

WRITECX, 18) 

READC x, x ) NA 
ENDIF 
NINC=2XNA 
SAMX=0 .0 
WRITEC 3 / 109) 

DO 33 1=1, NA 
WRITECX, 23) I 
READC X, x ) DALPHACI) 
WRITEC3, 110) I, DALPHACI) 
ALPHACI)=DTRX DALPHACI) 

SA = SINCALPHACI) ) 
IFCSA.GE.SAMX) SAMX=SA 
33 CONTINUE 

IFCLEVEL . EQ. '2' ) THEN 
WRITECX, 119) 

READC X, X) NT 
ELSE 

WRITECX, 20) 

READC X, X) NT 
ENDIF 

IFCLEVEL . EQ. • 2 • ) THEN 
WRITECX, 120) 

READC x , X ) NP 
ELSE 

WRITECX, 21) 

READC X, * ) NP 
ENDIF 

DO 44 1=1, NP 
WRITECX, 24) I 
READC X, X) DPHICI) 

44 PHIC I ) = DTRXDPHIC I) 

MSTART = 0 
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MST0P=INTCRMAX*SAMX)+1 
NR I TEC * , 121 ) MSTOP 
READC *, * ) MSTOP 
IFCSAMX. LT. 1 . E-6 ) MSTART = 1 

C SETTING RADIAL STEPS BETWEEN RA to R1 AND R2 to RB 

NR1 = 1 
NR2 = 1 

C DEFINING OUTSIDE PERMITTIVITY WITH SLIGHT LOSS 

URCNLAY+1)=(1 . 0,-1 . OE-6) 

ERCNLAY+l)=C1.0,-1.0E-6) 

C ESTIMATING NUMBER OF CORE AND EXTERIOR MODES 

N1=INT(RMIN*CABSCC$QRTC2 . 0*ERC 1 )*URC 1 ) ) ) )+l 
N2=INTCRMAX*CABSCCSQRTC2. G*ERCNLAY)*URCNLAY) )) )+l 
NRITEC * , 122) N1 
READC*,*) N1 
NRITEC*, 123) N2 
READC*,*) N2 
NRITEC*,*) BELL 

PAUSE 'Check that N1+N2 is ,LE. 35 .... Otherwise Abort 1 

IFC CN1+N2) . GT . 35 ) STOP 

WRITEC 3,105) MSTART, MSTOP, N1,N2 

NRITEC *,124) 

M1=MSTART+1 
M2=MST0P+1 
MM=MST OP-MSTART-1 
IFCMM.LE.O) MM = 1 
DO 77 MD=M1 , M2 

OPENC 1 , FI L E= ' RS . DAT 1 , STATUS 3 1 UNKNOWN ' , ACCESS = 1 SEQUENTIAL ' , 
1 FORM= ' UNFORMATTED ' ) 

DD=DMAX-DMIN 

M=MD-1 

WRITEC3, 101) FF 
WRITE(3,204) M 
NRITEC*,*) BELL 
NRITEC*, 125) M 
BC = 0 

IF CM.EQ.l) BC=1 
MF = 0 

IFCM.EQ.O) MF=1 

C COMPUTING MODAL EXPANSION TRUNCATION LIMITS 
N M 1 = N 1 + 1 - M - M F 
IFCNM1.LT. 1) NM1 = 1 
NM2 -N2+1 -M-MF 
IFCNM2.LT. 1) NM2=1 
NC0EF=2*C NM1 + NM2 ) 

NMODES=NCOEF+NI NC 

C COMPUTING MESH PARAMETERS AND RADCI) FOR GIVEN MODE = M 
DENS=DMAX-DD*CM+MF-1 . 0 )/MM 
NRITEC*, 126) 

CALL MESH C DENS, ER, UR, NR1 , NR2, IMX, NDS, LMX, NB, DR, DT, 

1 RA,R1,R2,RB,RMIN,RMAX,NLAY) 

C CHECKING MAXIMUM MESH PARAMETERS 

IFCNB.LE.NBMX) GO TO 55 
NRITEC 3, 10A) 

STOP 

55 CONTINUE 

C LOADING AND MARCHING CCR)) AND CCS)) 

NRI TEC 3 , 101) FF 
NRITEC 3,201) 

NRITEC*, 127) IMX 
DO 66 1=1, IMX 
I R= IMX-I 

NRITEC*, 128) I , I R 
NRITEC*, 129) 

CALL LODER C I ,M, BC, ER, UR, NMODES, NM1 , NM2, NINC, NB, IMX, IBIG, 

1 NDS,RA,RB,R1,R2,DR,DT,NR1,NR2,RAD,NLAY) 
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IFCI.GT.l) GO TO 57 

WRIT EC 3 , 30 0 ) I, (RAD(1,L),L=1,NLAY) 

57 CONTINUE 

WRITE(3,300) I, (RAD(2,L),L=1,NLAY) 

WRIT EC x, 130 ) 

CALL MARCH C I , IMX, BC ,NB, NBMX, NMODES, NMMAX, R, S, RINV) 

66 CONTINUE 

WRITECX, 131 ) 

CALL SWEEP CM, IMX, BC, NM1 , NM2, NINC, NMGDES, NMMAX, NCOEF, NCMAX, 

I NR1,NR2,NB,NBMX,R1,R2,DT, ER , UR, R, S, SI NV, FMOM, COEF, NINMX) 
WRITECX, 132 ) 

CALL MODECO C NM1 , NM2, N INC, NCOEF, NCMAX, FMOM, COEF, M, NINMX) 

WRITECX, 133) 

CALL RADFLD CM, NMI , NM2, NT, NP, NA, NCMAX, COEF, NINMX) 

WRIT EC X , 1 34 ) 

CALL CRSSEC C M, NM1 , NM2, NA, NCMAX , COEF, NI NMX) 

REWIND 10 

READC 10,133) IBIG, IDUMMY 
READC 10,103) RMIN, RMAX 
77 CONTINUE 

CL OSEC 1 , STATUS^ * DEL ETE 1 ) 

CL0SEC2) 

CLOSECIO) 

WRITECX , 135 ) 

CALL DATOUT C MSTART , MST OP , NT , NP , NA ) 

WRITECX, X) BELL 
WRITECX, 136) 

10 F0RMATC/////////////////7X, *xxxxxx WELCOME TO EMCAD xxxxxx * , //IX, 

I I THIS PROGRAM COMPUTES EM FIELD SCATTERING FROM PENET RABL E 1 , /7 X, 

2 * BODIES OF REVOLUTION USING THE COUPLED AZIMUTHAL POTENT I AL 1 ,/7X , 

3 * C CAP ) FORMULATION IN CONJUNCTION WITH A VARIATIONAL F I NITE- 1 , /7 X , 
4 f EL EMENT TECHNIQUE AND A TRI-REGIONAL UNIMOMENT METHOD. THE*,/7X, 
5 1 NUMERICAL SOLUTION IS PERFORMED USING A TWO-SWEEP RICCATI * ,/7X, 

6 1 TRANSFORM WITH CIRCUMFERENTIAL MARCHING . 1 ,/7X, 

7 1 UPDATED VERSION OF ORIGINAL U.C. BERKELEY CODE.',/7X, 

8* MODIFICATIONS BY M.A. MORGAN MAR 1987 - FEB 1983 f ,/7X, 

9 * E.M. CONNOLLY JUL 1987 - APR 1988 . *,//////7X, 

1 1 PL EASE PRESS ANY KEY TO CONTINUE.*,/) 

12 FORMAT C //////////////// 7X, 1 xxxxxxx EMCAD INPUT xxxxxxx » , //7X , 

1 * EMCAD ALLOWS A USER TO INPUT THE NECESSARY INFORMATION 1 , /7X, 

2 1 ACCORDING TO THE LEVEL OF HIS EXPERTISE . 1 ,//7X, 

3*1 NOVICE LEVEL - GIVES BRIEF EXPLANATIONS/DESCRIPTIONS 1 ,/7X, 

4* OF INPUT VALUES. DESCRIBES THE REQUIRED » , /7X , 

5* FORMAT FOR THE INPUT VALUE. GIVES f ,/7X, 

6* TYPICAL VALUES WHERE APPLICABLE.* , //7X, 

7*2 EXPERT LEVEL - ASSUMES THE USER IS FAMILIAR WITH THE*,/7X, 

8* INPUT PARAMETERS AND FORMATS. SIMPLY*, /7X, 

9* PROMPTS FOR REQUIRED INPUTS .', //7X, 

1*3 EXIT EMCAD*, //7X, 

2 * PLEASE SELECT THE LEVEL OF EXPERTISE BY ENTERING 1, 2, OR 3.', 
3//7X) 

13 FORMAT C//////////////////////7X, 

1 * THE OUTPUT DATA FILE IS THE FILE CONTAINING THE OUTPUT*, /7X, 

2 * RESULTS OF ALL NUMERICAL CALCULATIONS CONDUCTED BY EMCAD. *,/7X, 

3 * THE FORMAT FOR THIS INPUT IS FILENAME ONLY. NO EXTENSION IS*,/7X, 
4 » REQUIRED OR DESIRED. EMCAD AUTOMATICALLY APPENDS AN *,/7X, 

5 * EXTENSION OF .OUT TO YOUR FILENAME .*, //7X, 

6 * P L EAS E ENTER THE FILENAME OF THE OUTPUT DATA FILE.*,//7X) 

1 *THE OUTPUT GRAPHICS DATA FILE IS THE FILE CONTAINING THE *,/7X, 

2 ' OUTPUT DATA FROM EMCAD TO BE USED AS INPUT TO GRAPHING »,/7X, 

3 * ROUTINES . THE FORMAT FOR THIS INPUT IS FILENAME ONLY. NO*,/7X, 

4 * EXTENSION IS REQUIRED OR DESIRED. EMCAD AUTOMATICALLY *,/7X, 

5 * APPENDS AN EXTENSION OF . TMT , . TMP , . T ET , AND .TEP TO YOUR*,/7X, 

6 * FI L ENAME AS IT PRODUCES FOUR OUTPUT FILES FOR GRAPHICS .*, //7X, 

7*. TMT > TM INCIDENCE, F-THETA » , /7X, 
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8 ' . TMP > TM INCIDENCE, F-PHI',/7X, 

9 ' . TET -> TE INCIDENCE, F-THETA ' , /7X, 

I'.TET > TE INCIDENCE, F-PHI',//7X, 

2 ' PLEASE ENTER THE FILE NAME OF THE OUTPUT DATA FILE.',//7X) 

15 FORMAT (///////////////////////7X, 

1 ' THE GRAPHICS CAPTION IS A PERSONALIZED CAPTION ALLOWING THE',/7X, 
2 ' USER TO IDENTIFY THIS SET OF GRAPHS FROM ALL OTHER SETS.',/7X, 

3 ’ THE MAXIMUM LENGTH OF THIS CAPTION IS 69 CHARACTERS . ’ ,/7X, 

9’NOTE: WHEN USED WITH THE GRAPHICS PACKAGE, THE PROGRAM IS',/7X, 
5'ABLE TO DIFFERENTIATE BETWEEN UPPER CASE AND LOWER CASE',/7X, 
k * * //7v 

7 1 PI EASE ENTEr’aNY GRAPHICS CAPTION YOU DESIRE' ,//7X) 

16 FORMAT (//////////////////////7X, 

l’DMIN AND DMAX ARE PARAMETERS OF MESH DENSITY IN TERMS 0F’,/7X, 

2 ’ ELEMENTS/ INTERIOR LAMBDA. INPUT VALUES ARE EXPECTED TO BE’,/7X, 

3 ’ REAL , I.E. THE DECIMAL POINT MUST BE INCLUDED. ' ,/7X, 

9 'TYPICAL VALUES ARE DMIN = 10. AND DMAX = 15.',//7X, 

9 ' PLEASE INPUT DMIN’,//7X) 

17 F0RMAT(//7X, 'PLEASE INPUT DMAX',//7X) 

18 FORMAT C//////////////////////7X, 

1 ' THE NUMBER OF INCIDENT FIELD ANGLES IS THE TOTAL NUMBER OF',/7X, 
2' INCIDENT FIELDS THAT IMPINGE ON THE OBJECT OF INTEREST. ' ,/7X, 
3'THIS PROGRAM ALLOWS A MAXIMUM OF FIVE INCIDENT FIELD ’,/7X, 

9 ’ ANGLES . WHEN ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE A',/7X, 

5 ' DECIMAL BECAUSE THE INPUT MUST BE IN INTEGER FORMAT, I.E.',//7X, 
5 » j^A = 3 , NA = 1 ' / /7X 

6 ' PL EASE INPUT THE NUMBER OF INCIDENT FIELD ANGLES. ' ,//7X) 

19 FORMAT (//////////////////////7X , 

1 ' THE INCIDENT ANGLES DEFINE THE DIRECTION FROM WHICH THE ',/7X, 

2 ' I NCI DENT FIELDS ORIGINATE .’ ,//7X, 

3*0 DEGREES > + Z DI RECTI ON ’, /7X, 

9*90 DEGREES > + X DIRECTION ' ,/7X, 

5*180 DEGREES > - Z DIRECTION' ,//7X) 

20 FORMAT C//////////////////////7X, 

1 'THE NUMBER OF SCATTERING FIELD THETA POINTS DETERMINES THE’,/7X, 

2 * SPACING BETWEEN THETA POINTS DURING EMCAD ITERATIONS AND',/7X, 

3 ' CALCULATIONS . ' ,//7X, 

9 * DELTA THETA = 180 / (NUMBER THETA POINTS - 1) S0...',//7X, 

5 ' NUMBER THETA POINTS = 37 ----- > DELTA THETA = 5 DEGREES ' ,//7X, 

6 ' NUMBER THETA POINTS = 19 ----- > DELTA THETA = 10 DEGREES ’ ,//7X, 

7 'WHEN ENTERING YOUR ANSWER, PLEASE DO NOT INCLUDE A DECIMAL ' ,/7X, 

8 ' BECAUSE THE INPUT MUST BE IN INTEGER FORMAT. I.E. NT = 19', ///IX, 
9 * PL EASE INPUT THE NUMBER OF SCATTERING FIELD THETA POINTS' ,//7X) 

21 FORMAT (//////////////////////7X , 

1 ' THE NUMBER OF PHI ANGLES IS THE TOTAL NUMBER OF PHI ANGLES .' ,/7X, 
2'THIS PROGRAM ALLOWS A MAXIMUM OF EIGHT PHI ANGLES. WHEN ',/7X, 

3 ' ENTERING YOUR ANSWER PLEASE DO NOT INCLUDE A DECIMAL ',/7X, 

9 ' BECAUSE THE INPUT MUST BE IN INTEGER FORMAT. I.E.',//7X, 

5' NP = 3 , NP = 1 ' ,//7X, 

6 ' PLEASE INPUT THE NUMBER OF PHI ANGL ES . ' , //7X ) 

23 F0RMATC/7X, 'ENTER INC FLD ANGLE (DEG) FOR f ’,I3,//7X) 

29 F0RMAT(/7X, 'ENTER PHI FLD ANGLE (DEG) FOR # *,I3,//7X) 

25 F0RMAT(//////7X, 'PLEASE PRESS ANY KEY TO CONTINUE.',////) 

100 FORMAT(/, ' xxxxxxxxxxxxxxxxxxxx EMCAD OUTPUT DATA xxxxxxxxxxxxxxxx 
ixxxxxxxxxxx ' ) 

101 FORMAT (A) 

102 FORMAT ( 15 ) 

103 FORMAT ( 9( E19 . 6 ) ) 

109 FORMAT ( ' *** PROGRAM ABORTED BECAUSE NB > NBMX xxx’ ) 

105 F0RMAT(/7X, 'SERIES SOLUTION PARAMETERS '/7X, 'MSTART= ', 15, /7X, 

1 ' MSTOP= ' ,I5,/7X, 'Nl= ', I8,/7X, 'N2= M 8 ) 

107 FORMAT (//7X,I5,2(2X, ' ( ' , 1PE11 . 3,2X, 1PE11 .3, ' ) ' )) 

109 F0RMAT(//7X, 'INCIDENT FIELD ANGL ES ', //9X, ’ N 1 OX, ’ ALPHA ( N) ' ) 

110 FORMAT (//7X,I3,7X,F5.0, ' DEG') 

111 F0RMAT(/7X, 'COMPLEX Er(n) AND Ur(n)') 

112 F0RMAT(//7X, 'ENTER STRUCTURE DATA FILE ( D : FI L ENAME . EXTENSI ON ) : ') 
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113 F0RMATC//7X, ’ENTER OUTPUT DATA FILENAME WITHOUT EXTENSION ’ , //7X , 

1 ’ FOR PRINTER OUTPUT ENTER "LPT1”: f ) 

114 F0RMATC//7X, ’ENTER GRAPHICS DATA FILENAME WITHOUT EXTENSION ’ ) 

115 F0RMATC//7X, ’ENTER CAPTI Ofi OR LABEL f ) 

116 F0RMATC//7X, ’ENTER DMIN ( EL EMENTS/ INTERIOR LAMBDA): r ) 

117 FORMATC//7X, ’ENTER DMAX (ELEMENTS/INTERIOR LAMBDA): f ) 

113 FORMATC//7X, 1 ENTER NUMBER OF INCIDENT FIELD ANGLES ( . L E . 5) ’) 

119 F0RMATC//7X, 'ENTER NUMBER OF THETA POINTS FOR SCAT FIELD C.LE 37.) 
l f ,4X) 

120 F0RMATC//7X, ’ENTER NUMBER OF PHI ANGLES C.LE. 3): f ) 

121 F0RMATC//7X, 1 ENTER MSTOP C.LE. 13)’, //7X, 

1 1 ESTIMATED "MINIMUM” VALUE IS : 1 , 4X, I 4, //7X) 

122 F0RMATC//7X, ’ENTERING NO. OF INTERNAL AND EXTERNAL EXPANSION MODES 
1 1 ,//7x, 

2 1 ESTIMATED CKlxRMIN) "MINIMUM” INTERNAL N1 = 1 , 2X, 14, //7X, 

3 T ENTER N1 C.GE. MSTOP): ’,//7X) 

123 F0RMATC//7X, ’ESTIMATED CKO*RMIN) ’’MINIMUM” EXTERNAL N2 =’,2X,I4, 

1//7X, ’ENTER N2 C.GE. MSTOP): ’,//7X) 

124 FORMATC//7X, ’INDEXING PROGRAM THROUGH VALUES OF M’) 

125 F0RMATC7X, ’M-LOOP . ... M = M5) 

126 FORMAT C7X, ’CALL MESH’) 

127 FORMAT C7X, ’ ENTER I-LGOP; NO STEPS TO COMPLETE: M6) 

123 F0RMATC7X, ’ I = ’ , 2X , 1 5 , 7X, ’NO STEPS TO G0:’,2X,I5) 

129 FORMAT C 7X , ’CALL LODER’) 

130 FORMAT C 7 X , ’CALL MARCH’) 

131 FORMAT C 7X, ’EX I-LOOP, SWEEP’) 

132 FORMAT C7X, ’CALL MODECO’) 

133 FORMAT C 7X, ’CALL RADFL D ’ ) 

134 FORMAT ( 7X , ’ CALL CRSSEC’) 

135 FORMAT C 7 X , ’EX M-LOOP, DATOUT’) 

136 F0RMATC//7X, ’ xxxxxxxxxxxxxxx EMCAD COMPLETED xxxxxxxxxxxxxx ’ ) 

137 FORMATC//7X, ’ENTER MATERIAL PARAMETER FILE C D : FI L ENAME . EXTENSI ON ) 

1 : ’ ) 

138 FORMAT C 2 C I 5 ) ) 

201 FORMAT C//7X, ’ R1CCATI MARCHING STEPS’) 

204 FORMAT C//7X, ’ xxxxxxxxxxxxxx M = ’,15,’ xxxxxxxxxxxxxx ’ // ) 

300 FORMATC7X, I5,5E12.3) 

99 STOP 
END 



The above constitutes only the main program of EMCAD. The 
remaining subroutines would take approximately 26 more pages to 
print. The entire program code can be obtained by sending an 
IBM PC compatible formatted 5 1/4 ” floppy disk to 

Prof. M. A. Morgan 
Code 62 MW 

Naval Postgraduate School 
Monterey, Ca . 93943 
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APPENDIX D 



EMCADIN USER'S MANUAL 

Welcome to EMCADIN. This program is designed to accept 
properly formatted input data from the West Coast 
Consultants commercial shareware product CURVE DIGITIZER 
(Appendix M) and convert it to a form which can be used by 
the scattering calculation program EMCAD. CURVE DIGITIZER 
is a computer aided design (CAD) package which will allow 
user manipulation of the data. 

A. BACKGROUND 

EMCADIN was written in the programming language of Ryan- 
McFarland FORTRAN (Appendix M) . Any questions concerning 
formats, changes to the program, etc... should be referred 
to the respective user's manuals. Dr. M. A. Morgan or LT. E. 
M. Connolly (Appendix M) . 

B. COMPUTER SETUP 

EMCADIN was not written for any specific computer 
configuration, and as such, should run on most IBM 
compatible machines. 

C. PROGRAM INPUT AND OUTPUT 

EMCADIN obtains its input from a previously generated 
data file. The input data file for EMCADIN must be located 
in the same directory/subdirectory as EMCADIN. This data 
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file is produced through the utilization of the computer 
aided design package, CURVE DIGITIZER. As shown in Figure 
D-l , the output of CURVE DIGITIZER is the input to EMCADIN. 
For EMCADIN and CURVE DIGITIZER to be compatible, it was 
necessary to establish and adhere to the predefined data 
format of CURVE DIGITIZER. In order for EMCADIN to be 
utilized successfully, it is imperative that the user 
understand the need for the particular format being used, as 
well as the format itself. The transition from CURVE 
DIGITIZER to EMCADIN provides more of an opportunity for 
catastrophic error than at any other point in this entire 
package of interactive computer programs. 

Understanding the format of the input of the data for 
EMCADIN necessitates a thorough understanding of some of 
the idiosyncracies of CURVE DIGITIZER, referred to 
henceforth for brevity as CD. Computer aided graphic design 
deals in the basic business of handling points. Any graphic 
figure can be represented as some combination of points. 
Lines can be formed by connecting two points. Arcs can be 
represented by placing more points closer together and by 
approximating the curvature with sufficiently small line 
segments between points. The required input-output format 
for CD is utilized to insure that adequate information is 
present to direct CD how to connect the points. A similar 
process may be that of the simple children's game connect- 
the-dots. A series of dots exists which define a picture, 
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CRVDIGIN AND EMCAOIN 
FLOW DIAGRAM 



if connected properly. As more and more dots are used, 
more complex pictures with better resolution can be 
constructed. The fundamental assumption employed by CD is 
that each point should be connected to the next point listed 
in the data file unless some code is present to indicate 
otherwise. For simplicity and unrestricted accessibility, 
CD input and output data are represented in terms of XY 
rectangular coordinates. CD flags the break between two 
data points by what shall be referred to as a separation 
code. Since two coordinates, X and Y, are required to 
establish a data point, two separation codes, one for the X 
data and one for the Y data are required to separate two 
distinct data points. The separation code employed is of 
the form 999990. The last two digits, 90, may change to 
reflect scaling, color, or any of a number of other things. 
It was not possible to determine the specifics from the user 
manual. In the specific case of EMCADIN input data, 999990 
999990 codes are required to distinguish the structural 
data describing one layer from that describing another 

r 

layer. 

CD stores the data in the exact order in which the data 
points were drawn. There are other methods of generating 
the simple geometric shapes demonstrated in the next two 
examples, however, these examples were specifically chosen 
to illustrate potential pitfalls that could occur if extreme 
care is not taken. For the first example, the simple square 
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of Figure D-2 was constructed utilizing CD. The base of the 
square was drawn by defining point number 1, and then by 
drawing a line to point number 2. The top was likewise 
established by defining point number 4, and then by drawing 
a line to point number 3. The two sides were then formed by 
connecting point number 4 to point number 1 with a line, and 
by connecting point number 2 to point number 3 with a line. 
This is one of many methods of connecting the four data 
points that define this square. Figure D-3 lists the output 
data file which contains the CD defining data for this 
square. Notice there are significantly more than four data 
points describing this object which requires only four data 
points. Recall from earlier discussion that EMCADIN 
utilizes the 999990 codes to distinguish structural data 
describing one layer from that describing another layer. 
Figure D-2 depicts only one layer, yet EMCADIN using the 
data file of Figure D-3 would interpret the data as four 
distinct layers. The 999990 codes of Figure D-3 are present 
due to the order in which the data points were entered. Had 
the square been generated by defining point number 3, 
connecting it to point number 2, then to 1, then to 4, 
there would have been no 999990 codes separating the data, 
and the data file of Figure D-3 would consist of only four 
pairs (XY coordinates) of data. 

The second and third major errors that might occur can 
also be demonstrated using Figure D-2. The sole 
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FIG. D-3. EXAMPLE OF AN OUTPUT 
DATA FILE FOR A CD GENERATED 
SQUARE 
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justification for utilizing CD was its ability to generate 
data to be used ultimately by the scattering program EMCAD. 
The main body of the thesis discusses the requirement of 
EMCAD for axisymmetric bodies of revolution. Any three 
dimensional axisymmetric body of revolution can be defined 
by knowing the axis of revolution, and then by defining the 
body in planar terms. As an illustration, a given sphere is 
a three dimensional body of revolution. The sphere has 
axial symmetry and can be described by a semi-circle in any 
meridian plane intersecting the axis. Rotating this planar 
semicircle about the axis would produce the original 
sphere. EMCAD requires this meridian representation of the 
three dimensional body as its input. For simplicity's sake, 
the CD representation of the meridian plane will be referred 
to as the XY plane, with X being the horizontal axis and Y 
being the vertical axis, as well as the axis of revolution. 
This coordinate system is for discussion purposes only, and 
is not to be confused with the three dimensional coordinate 
system of EMCAD theory. Viewing the. body of Figure D-2 as a 
two dimensional body, the Y axis is the required axis of 
revolution. Therefore all CD planar representations for 
EMCAD should be constructed with values on the X axis 
restricted to positive values. 

The third problem occurs when removal of data points is 
attempted in order to obtain a proper planar representation. 
EMCADIN performs a valuable service to the user of improving 
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the resolution of input data. It suffices to say that due 
to the inherent mathematics involved, EMCAD requires that 
two data points be defined, both on the Y axis. The first 
data point defines the beginning of the object, at X = 0, 
and the last data point defines the end of the object, also 
at X = 0. In the example of the square, the beginning data 
point would be (0,2), and the ending data point would be 
(0,-2). Without these beginning and ending data points for 
each layer of data, the EMCADIN interpolation will "bomb", 
and will produce erroneous results, if any results at all. 

For the second example, the simple circle of Figure D-4 
was chosen. CD allows the user to construct a circle by 
simply defining the coordinates of the center of the circle, 
and by specifying the desired radius. In every instance, CD 
generates the data points by starting at point 3 on Figure 
D-4, progressing to point 2, to point 5, to point 4, and 
back to point 3 in a counter clockwise pattern. A sample of 
the output data file for this circle is shown in Figure D-5. 
This data, for a chosen radius of 2, reflects the same 
counter clockwise pattern, which is directly opposite to the 
pattern required by EMCADIN input data. In EMCADIN, theta 
is defined in a clockwise pattern, starting on the positive 
y axis and progressing by equal step sizes clockwise to the 
negative y axis. By using this method of circle 
construction, the user must find someway to edit and 
rearrange the data before it can be used by CD. The other 
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FIG. 0-4. EXAMPLE OF CD GENERATION 
OF A SIMPLE CIRCLE 
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FIG. D-5. EXAMPLE OF A PORTION 
OF AN OUTPUT DATA FILE FOR A CD 
GENERATED CIRCLE 
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problem that arises, particularly with circles, is that as 
the radius becomes smaller and smaller, fewer and fewer 
points are used to construct the circle, causing the data 
points to become less and less accurate. 

For the case of multiple layered bodies, EMCADIN 
expects that the data point values for each layer have been 
entered in the input data file starting with the layer 
closest to the center of the object, labeled as Layer 1 in 
Figure D-6, then advancing outward from the center until all 
layers have been accounted for. It has already been stated 
that data points are listed in the order in which they are 
entered, which means that not only might data points in a 
given layer be out of proper order, but data points in 
differing layers, or even the layers themselves might be 
interspersed. At the present time, there appears to be only 
two methods to preclude the abundance of possible errors. 
The first of these dictates the use of extreme care when 
entering the data points. The values must be entered in 
order, starting at e = 0 degrees and progressing clockwise 
to 0 = 180 degrees for the layer closest to the center 
point, and then entered for the next closest layer, and then 
the next closest layer . . . until all layers have been 
described. Even a small band defined by a different 
material parameter must be defined for all 0 from 0 degrees 
to 180 degrees. If this process is followed exactly, there 
is no need for the second method. The second method 
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EMCADIN VIEW OF LAYERED DATA 



involves the use of some text editor to rearrange the data 
as necessary, remove unnecessary 999990 codes, insure data 
for layers as well as layers are in the proper order, etc. . . 
The prudent method would include extreme caution coupled 
with a data review utilizing a text editor. 

The EMCADIN program output data consists of the number 
of layers of the body being considered, the minimum and 
maximum values of radius, and the listing of radius values 
in a clockwise manner for 9 = 0 to e = 180 degrees for each 
layer. 

D. PROGRAM EXECUTION 

Once the correct disk or directory has been selected, 
EMCADIN can be invoked by typing 

EMCADIN 

at the DOS prompt. The program will load itself and the 
user will see the following screens: 

* Welcome to EMCADIN (Figure D-7) 

This screen presents some brief background information. 

* Input Filename (Figure D-8) 

This screen prompts the user for the name of the file 
containing the input data for EMCADIN. If the file is 
obtained from CD, the file must have an extension of .dat, 
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WELCOME 10 EMCADIN 

THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA 
FROM CURVE DIGITIZER AMD CONVERIS IT TO A FORM 
WHICH CAN BE USED BY ENCAD, 



PLEASE PRESS ENTER TO CONTINUE. 



FIG. 0-7. WELCOME TO 
EMCADIN MENU 
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IHPUI FILENAME 



1NPUTFN IS THE INPUT FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED ON AND 
CONUERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BV 
ENCAD. 



PLEASE ENTER THE NAME OF THE INPUT FILE. THE EXTENSION 
OF THE FILENAME NUSI BE .DAI, I.E. F1LENANE.DAT 



FIG. D-8. EMCADIN PROMPT 
FOR INPUT FILENAME 
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and the extension must be entered by the user as part of the 
filename. The input data file must be in the same 
directory/subdirectory as EMCADIN. 

* Output Filename (Figure D-9) 

This screen prompts the user for the name of the file that 
will hold the output data from EMCADIN. The extension of 
this file is left to user discretion, but must be included 
as part of the output filename. 

* Resolution (Figure D-10) 

EMCADIN uses linear interpolation to provide the capability 
to generate output data with better resolution than the 
input data. This screen prompts the user for the desired 
output data resolution, in the units of degrees. 
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OUTPUT FILENAME 



OUTPUTFN IS IHE OUTPUT FILENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND CONVERTED TO A 
FORM WHICH WILL BE USED BY EMCAD. 



PLEASE ENTER THE NAME OF THE OUIPUT FILE. THE EXTENSION 
OF THE FILENAME MUST BE INCLUDED I.E. FILENAME.DAT 



FIG. D-9. EMCADIN PROMPT 
FOR OUTPUT FILENAME 
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RESOLUTION 



DELTHE IS THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES. 

PLEASE ENTER THE DESIRED DELTA THETA VALUE IN DEGREES . 



FIG. D-10. EMCAOIN PROMPT 
FOR OUTPUT DATA RESOLUTION 
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APPENDIX E 



EMCADIN SOURCE CODE 



PROGRAM EMCADIN 

********************* ******************* ******** ******** ********** 
WELCOME TO EMCADIN WRITTEN BY LT E. M. CONNOLLY. 

EMCADIN IS A PROGRAM DESIGNED TO ACCEPT PROPERLY FORMATTED 
INPUT DATA FROM CURVE DIGITIZER AND TO CONVERT IT TO A 
FORM WHICH CAN BE USED BY THE PROGRAM EMCAD, WRITTEN BY 
DR. M. A. MORGAN . 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE DEFINITIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



INPUTFN 

OUTPUTFN 

1/ J, K, 
L, N 
RADMIN 
RADMAX 
PI 

DELTHE 

DELTHETA 

THETACN) 

DEGTORAD 

NUMPTSIN 

NUMSETS 

DUMMY1 

DUMMY2 

XARRAYCI) 

YARRAYCI) 

A(I, J) 

BCI, J) 

NUMPTOUT 

THETAI 

THETAK 

XFOUND 

YFOUND 



THE INPUT FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED 
ON AND CONVERTED TO A FORM WHICH CAN 
SUBSEQUENTLY BE USED BY EMCAD 
THE OUTPUT FILENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND 
CONVERTED TO A FORM WHICH WILL BE USED BY 
EMCAD 

INTEGER COUNTERS 

THE MINIMUM RADIUS VALUE 
THE MAXIMUM RADIUS VALUE 
A CONSTANT 

THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES 

THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN RADIANS 

AN ARRAY HOLDING THE DESIRED VALUES OF THETA 
AS DETERMINED BY DELTHE 

A CONVERSION FACTOR FROM DEGREES TO RADIANS 
THE NUMBER OF POINTS HELD IN A GIVEN INPUT 
DATA SET 

THE NUMBER OF INPUT DATA SETS 

AN INPUT VARIABLE ALLOWING SCREEN ADVANCING 

DURING MENU SELECTION 

A CHARACTER STRING USED TO LOCATE THE END OF 
EACH INPUT DATA SET 

AN ARRAY CONTAINING THE VALUES OF X FOR ALL 
THE INPUT DATA SETS 

AN ARRAY CONTAINING THE VALUES OF Y FOR ALL 
THE INPUT DATA SETS 

AN ARRAY CONTAINING THE VALUES OF X FOR ALL 
THE J INPUT DATA SETS 

AN ARRAY CONTAINING THE VALUES OF Y FOR ALL 
THE J INPUT DATA SETS 

THE NUMBER OF POINTS OUT, DETERMINED BY THE 
DESIRED RESOLUTION 

THE VALUE OF THETA DETERMINED BY THE I-TH INPUT 
X AND Y VALUES * 

THE VALUE OF THETA DETERMINED BY THE K-TH INPUT 
X AND Y VALUES 

THE DESIRED VALUE OF X CORRESPONDING TO A GIVEN 
VALUE OF THETACN) 

THE DESIRED VALUE OF Y CORRESPONDING TO A GIVEN 
VALUE OF THETACN) 
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RADIUS C N, J ) AN ARRAY CONTAINING THE OUTPUT CALCULATED 

VALUES OF RADIUS FOR ALL THE INPUT DATA SETS 
M SLOPE OF THE LINE CONTAINING THE I-TH AND K-TH 

DATA POINTS 

BB Y-IN7ERCEPT OF THE LINE CONTAINING THE I-TH AND 

K-TH DATA POINTS 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX**** 

VARIABLE DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DOUBLE PRECISION XINPUTC 2000 ) , YINPUTC 2000 ) , RADIUS (2000,5) 

DOUBLE PRECISION DELTHE, DELTHETA, THETAC2000), .THETAI, THETAK 
DOUBLE PRECISION A(2000,5), BC2000,5) 

DOUBLE PRECISION DEGTORAD, INTER1, INTER2, RADMIN, RADMAX 
DOUBLE PRECISION PI, DUMMY, XFOUND, YFOUND 
DOUBLE PRECISION BB, M 

INTEGER I, J, K, L, N, NUMPTSIN, NUMPTOUT, NUMSETS 
CHARACTERS DUMMY1 , DUMMY2 
CHARACTERX12 INPUTFN , OUTPUTFN 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PI = 3.14159 
J = 0 
K = 0 
L = 0 

NUMPTSIN = 0 
NUMSETS = 1 
DEGTORAD = PI / 180. 

RADMIN = 10000 
RADMAX = -10000 

xxxxxx*x*xxx*****xxxxx****xxxxx*xxxx*x*x*xxxxxxxxxxxxxxxx*x*xx*x* 
INPUT DATA OF INTEREST 

x*xxxx*x*xxxxxxxx*xx**x**xx*xx*xx*x*******xxxxx*****x*xxxxxxxxxx* 

WRITE(*,110) 

READC * , 10 0) DUMMY1 
HR I TEC x, 120) 

READC *,100) INPUTFN 
WRITE(X,130) 

READC ^,100) OUTPUTFN 
HRITEC * , 140) 

READC*,*) DELTHE 

xxx**x****x****x*xx*xx*x**x*x*x****xxx*x*xxx*xx**xxxxxxxx*x*xxxxxx 

*x*x*x*xx**x*x*xx*xxxxxxxx***xxx*xx*xx***x*xxxxxxxxxxxxxxxxxxxxxxx 
DETERMINE LOCATION OF END OF DATA 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
OPENC7, FILE = INPUTFN) 

OPENC 8 , FI L E = OUTPUTFN) 

10 READ C 7 , 1 0 0 ) DUMMY2 

NUMPTSIN = NUMPTSIN + 1 
IFCDUMMY2 ,NE. 'END') GOTO 10 
CL0SEC7) 

OPENC 7 , FILE = INPUTFN) 

N = 1 
J = 1 

DO 20 1= 1, NUMPTSIN - 1 

READ C 7 , * ) XINPUTCI), YINPUTCI) 

IFCABSCXINPUTCI) - 999990) .GT. 20.) THEN 
A C N, J ) = XINPUTCI) 

BCN,J) = YINPUTCI) 

NUMSETS = J 
ELSE 
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J = J + 1 
N = 0 
ENDIF 
N = N + 1 
20 CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTERMEDIATE CALCULATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

NUMPTOUT = INTC180. / DELTHE) + 1 
DELTHETA = DELTHE X DEGTORAD 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTERPOLATION FOR THETA RESOLUTION 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DO 30 J = 1, NUMSETS 
DO 40 N = 1, NUMPTOUT 

DETERMINE THETA WHERE DATA POINTS X AND Y ARE TO BE FOUND 
THETA(N) = ( N - 1 ) x DELTHETA 
DO 50 I " 1, NUMPTSIN - 1 
K = 1-1 

I F ( I .EQ. 1) B(K,J) = BCI,J) 

IF( ABS( BC I , J ) ) .LT. .001) THEN 
THETAI = PI/2 

ELSE IFCBCI,J) .LT. 0) THEN 
THETAI = PI - DATANCDABSCACI,J)/BCI,J))) 

ELSE 

THETAI = DATANCABSCACI, J)/BCI, J))) 

ENDIF 

DETERMINE IF K-TH DATA PAIR FORMS ANGLE GREATER THAN 90 DEGREES 
ASSUMES CENTER OF COORDINATE SYSTEM AT X = 0 Y = 0 
I F ( ABS( B C K , J ) ) .LT. .001) THEN 
THETAK = PI/2 

ELSE IFCBCK,J) .LT. 0) THEN 
THETAK = PI - DATANCDABSCACK,J)/BCK,J))) 

ELSE 

THETAK = DATANCABSCACK, J)/BCK, J))) 

ENDIF 

CHECK TO SEE IF K-TH POINT IS THE DESIRED DATA POINT 
IFCDABSCTHETACN) - THETAK) .LT. .001) THEN 
XFOUND = ACK,J) 

YFOUND = B(K,J) 

RADI US C N , J ) = 2*PI*DSQRTCXF0UND**2 + YF0UND**2) 

I FC RADI USC N , J ) .LT. RADMIN) RADMIN = RADIUSCN,J) 

I F( RADI US ( N, J ) .GT. RADMAX) RADMAX = RADIUSCN,J) 

GOTO AO 

CHECK TO SEE IF I-TH POINT IS THE DESIRED DATA POINT 
ELSE IFCDABSCTHETACN) - THETAI) .LT. .001) THEN 
XFO M ND = A C I , J ) 

YFOUND = B C I , J ) 

RADIUSCN,J) = 2*PI*DSQRTCXF0UND**2 + YF0UND**2) 

IFCRADIUSCN, J) .LT. RADMIN) RADMIN = RADIU$CN,J) 

IFCRADIUSCN, J ) .GT. RADMAX) RADMAX = RADIUSCN,J) 

GOTO AO 

CHECK TO SEE IF THE DESIRED DATA POINT IS BETWEEN THE K-TH 
DATA POINT AND THE I-TH DATA POINT 
ELSE IF(THETACN) .GE. THETAK .AND. 

1 THETACN) ,LE. THETAI) THEN 

C CHECK FOR VERTICAL LINE 

IFCDABSCACI,J)-A(K,J)) .LT. .0001) THEN 
XFOUND = ACI,J) 

YFOUND = XFOUND / DTANC THETACN) ) 

RADIUS C N, J ) = 2*PI*DSQRTCXF0UND**2 + YFQUND**2) 

IFCRADIUSCN, J ) .LT. RADMIN) RADMIN = RADIUSCN,J) 

IFCRADIUSCN, J) .GT. RADMAX) RADMAX = RADIUSCN,J) 

GOTO AO 
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C CHECK FOR HORIZONTAL LINE 

ELSE IFCDABSCBCI, J)) .LT. .0001) THEN 
YFOUND = BCI,J) 

XFOUND = YFOUND X DTANC THETACN) ) 

RADIUSCN,J) = 2*PI*DSQRTCXF0UND**2 + YF0UND*X2) 

IFC RADIUSC N, J ) .LT. RADMIN) RADMIN = RADIUSCN,J) 

IFCRADIUSCN, J) . GT . RADMAX) RADMAX = RADIUSCN,J) 

GOTO AO 

C SLOPE OF LINE CONTAINING I-TH AND K-TH DATA POINTS 

ELSE 

M = C BCI,J)-BCK,J) ) / ( ACI, J)-ACK, J) ) 

C Y-XNTERCEPT OF LINE CONTAINING I-TH AND K-TH DATA POINTS 
BB = BCI,J) *= C M x A C I , J ) ) 

RECALL EQUATION OF ANY LINE IS Y = MX + B 

THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS INTERSECTS 

THE LINE DRAWN FROM THE ORIGIN THROUGH THE DESIRED DATA POINT 

AT THE DESIRED DATA POINT XFOUND AND YFOUND. THE EQUATION FOR 
THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS MUST THEN 
ALSO WORK AT THE DESIRED DATA POINT. WRITING XFOUND AND 
YFOUND IN TERMS OF THEIR ANGLE THETACN) , AND SUBSTITUTING 
YFOUND = BB / C 1 - C M * DTANC THETACN) ) ) ) 

XFOUND = YFOUND x DTANC THETACN) ) 

RADI US C N , J ) = 2XPJXDSQRTCXF0UNDXX2 + YF0UNDXX2 ) 

I FC RADIUSC N , J ) .LT. RADMIN) RADMIN = RADIUSCN,J) 

IFC RADIUSC N , J ) .GT. RADMAX) RADMAX = RADIUSCN,J) 

GOTO AO 
ENDIF 

NEED TO TRY A NEW I-TH AND K-TH DATA PAIR 
ELSE 

DUMMY = 0 
ENDIF 

50 CONTINUE 
40 CONTINUE 
30 CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
OUTPUT TO DATA FILE 

****************************************************************** 
WRIT EC 8 , 103) NUMPTOUT, NUMSETS + 1 
WRITEC 8 ,102) RADMIN, RADMAX 
DO 400 N = 1, NUMPTOUT 

WRITEC 8,101) N , C RADIUSC N , J ) , J = l, NUMSETS) 

400 CONTINUE 

****************************************************************** 
CLOSEC7) 

CL0SEC8) 

100 FORMAT C A ) 

101 FORMAT CI5,5(E14.6)) 

102 FORMAT C2CE14.6)) 

103 FORMAT C2C 15) ) 

110 FORMAT C////////////////////5X, 1 WELCOME TO EMCADI N 1 , //5X , 
l'THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA 1 ,/5X, 

2 1 FROM CURVE DIGITIZER AND CONVERTS IT TO A F0RM. f ,/5X, 

3 1 WHICH CAN BE USED BY EMCAD. »,/5X, 

4////5X, 1 PL EASE PRESS ENTER TO CONTI NUE . 1 , //5X ) 

120 FORMAT C////////////////////5X , 1 INPUT FI L ENAME 1 ,////// 5X, 

1 f I NPUTFN IS THE INPUT FILENAME OF THE FILE CONTAINING f , /5X, 

2 1 TH E DATA OF INTEREST WHICH WILL BE INTERPOLATED ON AND',/5X, 

3 1 CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BY f ,/5X, 

4 f EMCAD. ’ ,///5X, 

5 1 PL EASE ENTER THE NAME OF THE INPUT FILE. THE EXTENSION f ,/5X, 

6 1 OF THE FILENAME MUST BE .DAT, I.E. FI L ENAME . DAT 1 , //5X ) 

130 F0RMATC////////////////////5X, f OUTPUT FI L ENAME 1 , //////5X, 

1 1 OUTPUTFN IS THE OUTPUT FILENAME OF THE FILE CONTA I N I NG f , /5X , 

2 1 THE DATA WHICH HAS BEEN INTERPOLATED AND CONVERTED TO A T ,/5X, 
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3 f FORM WHICH WILL BE USED BY EMCA D . 1 , /// 5X , 

4 1 PL EASE ENTER THE NAME OF THE OUTPUT FILE. THE EXTENSION 1 ,/5X, 
6 1 OF THE FILENAME MUST BE INCLUDED I.E. FILENAME.DAT 1 ,//5X) 

140 F0RMATC////////////////////5X, RESOLUTION 1 ,////// 5X , 



1 ' DEL THE IS THE USER I T4PUT VALUE 
2 1 RESOLUT I ON IN DEGREES . 1 , ///5X , 
3 T PLEASE ENTER THE DESIRED DELTA 
END ‘ 



OF THE DESIRED THETA 1 ,/5X, 
THETA VALUE IN DEGREES . ' ,// 5X ) 
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APPENDIX F 



CRVDIGIN USER'S MANUAL 

Welcome to CRVDIGIN. This program is designed to 
accept properly formatted input data, and to convert it to a 
form which can be used by the West Coast Consultants 
commercial shareware product CURVE DIGITIZER (Appendix M) . 
CURVE DIGITIZER is a computer aided design (CAD) package 
which will allow user manipulation of the data. Upon 
completion of using CURVE DIGITIZER, the user may employ the 
program EMCADIN to transform the data into a form which can 
be used by the scattering calculation program EMCAD. 

A. BACKGROUND 

EMSCAT was written by Dr. M. A. Morgan while at the 
University of California at Berkeley. Extensive interactive 
computer aided design type feature have been added to 
EMSCAT, which consider scattering from inhomogeneous 
axisymmetric objects. These computer aided design features, 
driving the change of program name from EMSCAT to EMCAD, 
allow structures composed of multiple dielectric materials 
to be graphically constructed or modified by the user, and 
then evaluated by EMCAD for the scattering. CRVDIGIN was 
written in the programming language of Ryan-McFarland 
FORTRAN (Appendix M) . Any questions concerning formats, 
changes to the program, etc. should be referred to the 
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respective user's manuals, Dr. M. A. Morgan or LT. E. M. 
Connolly (Appendix M) . 

B. COMPUTER SETUP 

CRVDIGIN was not written for any specific computer 
configuration, and as such, should run on most IBM 
compatible machines. 

C. PROGRAM INPUT AND OUTPUT 

CRVDIGIN obtains its input from a previously generated 
data file. The input data file for CRVDIGIN must be located 
in the same directory/ subdirectory as CRVDIGIN. This data 
file may be produced by one of two methods. The first 
method consists of special case computer programs which 

generate data files describing the penetrable body of 
revolution for the specific geometrically shaped object 
being considered. A few examples of these special case 
programs are SPHERE, CONE and CYLINDER listed in Appendices 
J, K and L respectively. Each of these programs generates 
two descriptive data files which can be separated into 

distinct areas of information. The first area of 

information is the data which conveys the material 

composition of each layer of the structure. This data is of 
no concern at this particular time, however, it will be 
required before running EMCAD. The second area of 
information, which is the data required as input for 
CRVDIGIN, is the data which defines the structural shape of 
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each layer of the penetrable body. This data is referred to 
as "STRUCTUR" . DAT, and is shown in the flow diagram of 
Figure F-l. The second method of input data file 
construction consists of user generation through direct 
keyboard input. 

As shown in Figure F-l, the output of CRVDIGIN, shown 
as "CURVE IN ".DAT, is the input to the CAD package CURVE 
DIGITIZER. For CRVDIGIN and CURVE DIGITIZER to be able to 
accept input data from a variety of sources, it was 
necessary to establish and adhere to a predefined data 
format. The special case computer programs were written to 
incorporate this format. In order for the second method to 
be utilized successfully, as well as for other special case 
programs to be written, it is imperative that the user 
understand the need for the particular format being used, as 
well as the format itself. The input format for CRVDIGIN is 
relatively uncomplicated. An example of a portion of an 
input data file for CRVDIGIN is shown in Figure F-2. The 
first line of data contains two numbers in integer format. 
The first integer represents the number of data points used 
to describe any layer. The second integer lists the total 
number of layers + 1. The '+ l' was used to correct some 
minor discrepancies during the design of the code. The data 
defining the physical structure of the boundary for each 
layer is written in radius versus theta format. The value 
of the angle theta associated with each radius is determined 
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CRVDIGIN AND EMCADIN 
FLOW DIAGRAM 



361 4 

0 . 250300E+01 0 . 382760E+02 

1 0 . 188496E+02 0.251327E+02 



2 0.176975E+02 

3 0 . 166793E+02 

4 0 . 157731E+02 

5 0 . 149613E+02 

6 0 . 14230 1E+02 

7 0.135680E+02 

8 0 . 129657E+02 

9 0 . 1 24 155E+02 

10 0 . 119109E+02 

11 0 . 1 14466E+02 

12 0 . 110180E+02 

13 0 . 10621 1E+02 

14 0 . 102525E+02 

15 0 . 990939B+01 

16 0 „ 958921E+01 

17 0 . 928975E+01 

18 0 . 900910E+01 

19 0 . 874555E+01 

20 0 . 849761E+01 

21 0 . 826396E+01 

22 0 . 804341E+01 

23 0. 783490E+01 

24 0 . 763750E+01 

25 0 . 745035E+01 

26 0 . 727270E+01 

27 0 . 710385E+01 

28 0 . 694318E+01 

29 0 . 679012E+01 

30 0 . 664416E+01 



0 . 241 175E+02 
0. 231827E+02 
0 . 223194E+02 
0. 215196E+02 
0.207767E+02 
0 . 200848E+02 
0. 194390E+02 
0. 188348E+02 
0. 182684E+02 
0. 177363B+02 
0. 172357E+02 
0 . 167638E+02 
0. 163182E+02 
0. 158969E+02 
0. 154980E+02 
0. 151197E+02 
0. 147605E+02 
0. 144191E+02 
0. 140941E+02 
0. 137845E+02 
0. 134893E+02 
0. 132074E+02 
0. 129380E+02 
0 . 126803E+02 
0. 124336E+02 
0. 121972E+02 
0. 119706E+02 
0. 117531E+02 
0. 115442E+02 



0.314159E+02 
0 . 304952E+02 
0 . 296292E+02 
0 . 288130E+02 
0 . 280428E+02 
0 . 273146E+02 
0 . 266253E+02 
0. 259719E+02 
0 . 253516E+02 
0 . 247621E+02 
0 . 242013E+02 
0 . 236670E+02 
0 . 231575E+02 
0 . 226712E+02 
0 . 222066E+02 
0 . 217622E+02 
0 . 213369E+02 
0. 208294B+02 
0 . 205387E+02 
0 . 201639E+02 
0. 198040E+02 
0 . 194582E+02 
0. 1912S6E+02 
0. 188057E+02 
0. 184978E+02 
0. 182009E+02 
0. 179149E+02 
0. 176390E+02 
0. 173728E+02 
0. 171158E+02 



FIG. F-2. PORTION OF SAMPLE 
INPUT DATA FOR CRVDIGIN.EXE 
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by starting at theta equal to 0 degrees and progressing to 
theta equals 180 degrees in step sizes determined by the 
number of points available. The program expects that the 
values of the radius for each layer have been entered 
starting with the layer closest to the center of the 
object, then advancing outward from the center until all 
layers have been accounted for. The second line of data 
contains two real values. The first real value is the 
minimum radius value present, and the second value is the 
maximum radius value present. Following the description for 
layer entry above, this implies that the minimum radius 
value should come from the first layer, and the maximum 
radius value should come from the last layer. Figure F-3 
demonstrates the indexing of layers. The next A lines of 
data consists of B columns of data. The value of A is the 
number of points, and the value of B is the number of layers 
+ 1. The first column value is an index value starting at 1 
and incrementing by 1. The second, third, ... columns 
contain values of the radius associated with theta = (index 
- 1) * deltatheta for the first layer, the second layer . . . 

For each successive row the index changes, the value of 
theta changes and the radii associated with that given value 
of theta are listed. 

The output of CRVDIGIN , used as an input to CURVE 
DIGITIZER, is just as simple in appearance as the input to 
CRVDIGIN, yet much more difficult in theory. The 
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FIG. F -3. GRAPHICAL REPRESENTATION OF 
ASSIGNMENT OF LAYER BOUNDARIES 



explanation of this difficulty actually is embedded in the 
operation of CURVE DIGITIZER. Although these comments may 
apply to other CAD packages, they are specifically directed 
towards CURVE DIGITIZER, referred to henceforth for brevity 
as CD. Computer aided graphic design deals in the basic 
business of handling points. Any graphic figure can be 
represented as some combination of points. Lines can be 
formed by connecting two points. Arcs can be represented by 
placing more points closer together and by approximating the 
curvature with sufficiently small line segments between 
points. The required input-output format for CD is 
utilized to insure that adequate information is present to 
direct CD how to connect the points. A similar process may 
be that of the simple children's game connect-the-dots . A 
series of dots exists which define a picture, if connected 
properly. As more and more dots are used, more complex 
pictures with better resolution can be constructed. The 
fundamental assumption employed by CD is that each point 
should be connected to the next point listed in the data 
file unless some code is present to indicate otherwise. For 
simplicity and unrestricted accessibility, CD input and 
output data are represented in terms of XY rectangular 
coordinates. CD flags the break between two data points by 
what shall be referred to as a separation code. Since two 
coordinates, X and Y, are required to establish a data 
point, two separation codes, one for the X data and one for 



135 



the Y data are required to separate two distinct data 
points. The separation code employed is of the form 999990. 
The last two digits may change to reflect scaling, color, or 
any of a number of other things. It was not possible to 
determine the specifics from the user manual. In the 
specific case of CRVDIGIN output data, 999990 999990 codes 

are used to distinguish the structural data describing one 
layer from that describing another layer. Other than the 
insertion of separation codes at the appropriate locations, 
CRVDIGIN performs the simple mathematical task of converting 
the input data generated in radius-theta format to output 
data represented in XY format. 

D. PROGRAM EXECUTION 

Once the correct disk or directory has been selected, 
CRVDIGIN can be invoked by typing 

CRVDIGIN 

at the DOS prompt. The program will load itself and the 
user will see the following screens: 

* Welcome to CRVDIGIN (Figure F-4) 

This screen presents some brief background information. 
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WELCOME 10 CRUDIGIN 



THIS PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DAIA 
AND CONUERTS IT TO A FORM WHICH CAN BE USED BY THE 
COMMERCIAL SHAREWARE PRGDUCI CUM DIGITIZER. 

CURUE DIGITIZER IS A CAD PACKAGE WHICH WILL ALLOW USER 
MANIPULATION OF THE DATA. UPON COMPLETION OF 
CURUE DIGITIZER, THE USER MAY UIILIZE THE PROGRAM 
EMCADIN TO TRANSFORM THE DAIA INTO A FORM WHICH CAN 
BE USED BY THE PROGRAM EMCAD WRITTEN BY DR. M. A. MORGAN. 



PLEASE PRESS ENTER TO CONTINUE. 



FIG. F-4. WELCOME TO 
CRVDIGIN MENU 
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Input Filename 



(Figure F-5) 



* 



This screen prompts the user for the name of the file 
containing the input data for CRVDIGIN. The file may have 
any extension that the user desires, however, the extension 
must be entered by the user as part of the filename. The 
input data file must be in the same directory as CRVDIGIN. 

* Output Filename (Figure F-6) 

This screen prompts the user for the name of the file that 
will hold the output data from CRVDIGIN. Because of the 
input requirements of CURVE DIGITIZER, the extension of this 
file must be .dat, and must be included as part of the 
output filename. 

* Resolution (Figure F-7) 

CRVDIGIN uses linear interpolation to provide the capability 
to generate output data with better resolution than the 
input data. This screen prompts the user for the desired 
output data resolution, in the units of degrees. 

E . PROGRAM OUTPUT 

A portion of an example of CRVDIGIN output is shown in 
Figure F-8. Notice the 999990 999990 codes separating the 
different layers of data. 
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INPUT FILENAME 



1NPUTFN IS THE INPUI FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH MILL BE INTERPOLATED ON AND 
CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BY 
THE CURVE DIGITIZER PROGRAM, 



PLEASE ENTER THE NAME OF THE INPUT FILE. THE EXTENSION 
OF THE FILENAME MUSI BE INCLUDED, I.E. FILENAME.DAT 



FIG. F-5. CRVDIGIN PROMP 
FOR INPUT FILENAME 
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OUTPUT FILENAME 



OUTPUTFN IS THE OUTPUT FILENAME OF THE FILE CONTAINING 
IHE DATA WHICH HAS BEEN INTERPOLATED AND CONUERIED TO A 
FORM WHICH WILL BE USED BV THE CURVE DIGITIZER PROGRAM. 



PLEASE ENTER THE NAME OF THE OUTPUT FILE. THE EXTENSION 
OF THE FILENAME MUST BE .DAT, I.E. FILENAME.DAT 



FIG. F-6. CRVDIGIN PROMPT 
FOR OUTPUT FILENAME 
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RESOLUTION 



DELIHE IS THE USER INPUT OALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES. 

PLEASE ENTER THE DESIRED DELTA THETA OALUE IN DEGREES. 



FIG. F-7. CRVDIGIN PROMP 
FOR OUTPUT DATA RESOLUTION 



4.69975948 

4.35180950 

4.00664711 

3.66403151 

3.32372093 

2.98551798 

2.64917397 

2.31447673 

1.98121083 

1.64916909 

1.31813526 

0.98790431 

0.65828347 

0.32906568 

4 . 32930647E-05 
999990 

0 . 00000000E-01 
0.40459359 
0.75102264 
1.05115640 
1.31384826 
1.54581928 
1.75229287 
1.93734324 
2.10425401 
2.25564241 
2.39365172 



-18.84953690 

-18.84956930 

-18.84958080 

-18.84959030 

-18.84950830 

-18.84953120 

-18.84954830 

-18.84953880 

-18.84956740 

-18.84959410 

-18.84956740 

-18.84953120 

-18.84950830 

-18.84952740 

-18.84959980 

999990 

25.13270000 
23.17916870 
21.50649070 
20.05727580 
18.78891940 
17.66880800 
16 . 67196660 
15.77840710 
14.97255710 
14.24157710 
13.57508280 



FIG. 



F-8. PORTION OF SAMPLE 



OUTPUT DATA FROM CRVDIGIN 
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APPENDIX G 

CRVDIGI N SOURCE CODE 



PROGRAM CRVDIGIN 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WELCOME TO CRVDIGIN WRITTEN BY LT E. M. CONNOLLY. 

CRVDIGIN IS A PROGRAM DESIGNED TO ACCEPT PROPERLY FORMATTED 
INPUT DATA, AND CONVERT IT TO A FORM WHICH CAN BE USED BY 
THE COMMERCIAL SHAREWARE PRODUCT CURVE DIGITIZER. 

CURVE DIGITIZER IS A CAD PACKAGE WHICH WILL ALLOW USER 
MANIPULATION OF THE DATA. UPON COMPLETION OF CURVE DIGITIZER, 

THE USER MAY UTILIZE THE PROGRAM EMCAD1N, ALSO WRITTEN BY 
LT E. M. CONNOLLY TO TRANSFORM THE DATA INTO A FORM WHICH CAN 
BE USED BY THE PROGRAM EMCAD WRITTEN BY DR . M. A. MORGAN. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

VARIABLE DEFINITIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



INPUTFN 



OUTPUTFN 

X, J, K, 
L, N 

NUMBROWS 

NUMBCOLS 

SEPRATN 

RADARRAY 

MINRAD 

MAXRAD 

RADIUS 

PI 

DELTHE 

DELTHETA 

THETAOLD 

THETA 

DEGTORAD 

X 



THE INPUT FILENAME OF THE FILE CONTAINING 
THE DATA OF INTEREST WHICH WILL BE INTERPOLATED 
ON AND CONVERTED TO A FORM WHICH CAN 
SUBSEQUENTLY BE USED BY THE CURVE DIGITIZER 
PROGRAM 

THE OUTPUT FILENAME OF THE FILE CONTAINING 
THE DATA WHICH HAS BEEN INTERPOLATED AND 
CONVERTED TO A FORM WHICH WILL BE USED BY 
THE CURVE DIGITIZER PROGRAM 
INTEGER COUNTERS 

THE NUMBER OF ROWS CONTAINED IN THE DATA 
MATRIX IN THE FILE INPUTFN 

THE NUMBER OF COLUMNS CONTAINED IN THE DATA 
MATRIX IN THE FILE INPUTFN 

A FLAG USED BY CURVE DIGITIZER TO DISTINGUISH 
BETWEEN DATA POINTS OR SETS OF DATA 
A WORKING ARRAY USED TO HOLD ALL THE VALUES 
OF THE ARRAY SPECIFIED ABOVE BY NUMBROWS AND 
NUMBCOLS 

THE MINIMUM RADIUS VALUE, NOT USED BY THIS 
PROGRAM BUT REQUIRED FOR EMCAD 
THE MAXIMUM RADIUS VALUE, NOT USED BY THIS 
PROGRAM BUT REQUIRED FOR EMCAD 
A WORKING VARIABLE USED TO HOLD THE CURRENT 
VALUE OF THE RADIUS OBTAINED FROM RADARRAY 
A CONSTANT 

THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN DEGREES 

THE USER INPUT VALUE OF THE DESIRED THETA 
RESOLUTION IN RADIANS 

A WORKING VARIABLE USED TO HOLD THE CURRENT 
VALUE OF THETA OBTAINED FROM RADARRAY 
AN ARRAY HOLDING THE DESIRED VALUES OF THETA 
AS DETERMINED BY DELTHE 

A CONVERSION FACTOR FROM DEGREES TO RADIANS 
THE VALUE OF X OBTAINED BY CONVERTING THE POLAR 
COORDINATES RADIUS, THETA TO X,Y COORDINATES 
THE X , Y COORDINATE SYSTEM IS REQUIRED FOR THE 
DATA TO BE COMPATIBLE WITH THE CURVE DIGITIZER 
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PROGRAM 

Y THE VALUE OF Y OBTAINED BY CONVERTING THE POLAR 

COORDINATES RADIUS, THETA TO X,Y COORDINATES 
THE X , Y COORDINATE SYSTEM IS REQUIRED FOR THE 
DATA TO BE COMPATIBLE WITH THE CURVE DIGITIZER 
PROGRAM 

INTER1, INTER2 INTERMEDIATE CALCULATION VARIABLES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE DECLARATION 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
CHARACTER*12 INPUTFN, OUTPUTFN, DUMMY 
INTEGER I, J, K, L, N 

INTEGER NUMBROWS, NUMBCOLS, NUMPTS, SEPRATN 
REAL RADARRAY (1000,10) 

REAL MINRAD, MAXRAD, RADIUS, PI, DELTHE, DELTHETA, THETAOLD 
REAL THETA(IOOO), DEGTORAD, X, Y 
REAL INTER1 , INTER2 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PI = 3.14X59 
DEGTORAD = PI / ISO. 

SEPRATN = 999990 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT DATA OF INTEREST 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRXTE(X,110) 

READC * , 100) DUMMY 
WRITEC*,120) 

READC *,100) INPUTFN 
WRXTEC 130 ) 

READC *,100) OUTPUTFN 
WRITEC X, 140) 

READC * , x) DELTHE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
READ DATA FILE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
0PENC7,FILE = INPUTFN) 

0PENC8, FILE = OUTPUTFN) 

READC 7 , * ) NUMBROWS, NUMBCOLS 
READC 7 , X ) MINRAD, MAXRAD 

READC 7 , X) C C RADARRAY C X , J ) , J~1 , NUMBCOLS) , I = 1, NUMBROWS) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTERMEDIATE CALCULATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

NUMPTS = I NT C 180 . / DELTHE) + 1 

DELTAOLD = C130. / CNUMBROWS - D) * DEGTORAD 

DELTHETA = DELTHE * DEGTORAD 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTERPOLATION AND TRANSFORMATION 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
J = 1 

DO 10 K = 2, NUMBCOLS 
DO 20 N = 1, NUMPTS 
THETA C N ) = CN-1) X DELTHETA 
DO 30 I = 1, NUMBROWS 

THETAOLD = C RADARRAY C I , J ) - 1) X DELTAOLD 
IFC ABSCTHETAC N) - THETAOLD) ,LT. .00001 ) THEN 
RADIUS = RADARRAY C I , K ) 

X = RADIUS X SINC THETACN) ) 
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Y = RADIUS * C0S( THETACN) ) 

WRITE(8,X) X , Y 

GOTO 20 

ELSE IF (THETACN) .GE. (THETAOLD - DELTAOLD) .AND. 

1 THETACN) .LE. THETAOLD) THEN 
L = I - 1 

INTER1 = RADARRAY(L,K) - RADARRAY ( I , K ) 

INTER2 = (THETACN) - THETAOLD) / (-DELTAOLD) 

RADIUS = INTER1 X INTER2 4 RADARRAY ( I , K ) 

REMOVES NORMALIZING FACTOR OF 2XPI GENERATED BY SPHERE, CONE, 
ETC.., FACTOR IS RESTORED IN EMCADIN 
X = RADIUS X SINC THETACN) ) X 1 / C 2 X PI ) 

Y = RADIUS X COSC THETACN) ) x 1 / C 2 X PI ) 

WRITEC8 , X) X , Y 

GOTO 20 
ELSE 

DUMMY BRANCH, TRY NEXT VALUE OF THETAOLD 

RADIUS = 0 

ENDIF 

30 CONTINUE 
20 CONTINUE 

HR I TEC 8 , x ) SEPRATN, SEPRATN 
10 CONTINUE 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
100 FORMATC A ) 

110 FORMAT (//////////////////// 5X, 1 HE L COME TO CRVDIGIN * ,//5X, 

1 1 THI S PROGRAM ACCEPTS PROPERLY FORMATTED INPUT DATA * , /5X, 

2 * A N D CONVERTS IT TO A FORM WHICH CAN BE USED BY THE 1 ,/5X, 

3 'COMMERCIAL SHAREWARE PRODUCT CURVE DIGITIZER . 1 , /5X, 

A 1 CURV E DIGITIZER IS A CAD PACKAGE WHICH WILL ALLOW USER*,/5X, 

5 f MANIPULATION OF THE DATA. UPON COMPLETION 0F*,/5X, 

6 * CURVE DIGITIZER, THE USER MAY UTILIZE THE PROGRAM *,/5X, 

7 * EMCADIN TO TRANSFORM THE DATA INTO A FORM WHICH CAN f ,/5X, 

8 1 BE USED BY THE PROGRAM EMCAD WRITTEN BY DR. M. A. MORGAN. 1 , 
9////5X , 1 PLEASE PRESS ENTER TO CONTI NUE . 1 , //5X ) 

120 FORMAT C //////////////////// 5X, 'INPUT FILENAME* ,//////5X, 

1 * I NPUTFN IS THE INPUT FILENAME OF THE FILE CONTAINING 1 , /5X, 

2 1 THE DATA OF INTEREST WHICH WILL BE INTERPOLATED ON AND*,/5X, 

3 * CONVERTED TO A FORM WHICH CAN SUBSEQUENTLY BE USED BY',/5X, 

A 1 THE CURVE DIGITIZER PROGRAM .*, ///5X, 

5 * P L EASE ENTER THE NAME OF THE INPUT FILE. THE EXTENSION *, /5X, 

6 f OF THE FILENAME MUST BE INCLUDED, I.E FILENAME, DAT* ,//5X) 

130 F0RMATC////////////////////5X, * OUTPUT FI L ENAME 5 , //////5X, 

1 1 OUTPUT FN IS THE OUTPUT FILENAME OF THE FILE CONTAINING 1 ,/5X, 

2 * THE DATA WHICH HAS BEEN INTERPOLATED AND CONVERTED TO A»,/5X, 

3 * FORM WHICH WILL BE USED BY THE CURVE DIGITIZER PROGRAM .*, ///5X, 

A * PL EASE ENTER THE NAME OF THE OUTPUT FILE. THE EXTENSI ON » , /5X, 

6 1 OF THE FILENAME MUST BE .DAT, I.E. FI L ENAME . DAT *, //5X) 

1 AO FORMAT C////////////////////5X, * RESOLUT I ON f , //////5X, 

1 f DELTHE IS THE USER INPUT VALUE OF THE DESIRED THETA*,/5X, 

2 * RESOLUTION IN DEGREES .*, ///5X, 

3 1 PL EASE ENTER THE DESIRED DELTA THETA VALUE IN DEGREES .*, //5X) 
CLOSEC7) 

CL0SEC8) 

END 
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APPENDIX H 



MCADGRAF USER'S MANUAL 

Welcome to MCADGRAF. This program was written to be 
used in conjunction with the FORTRAN program entitled EMCAD. 
MCADGRAF adds graphics features to support the analysis of 
data generated and formatted by EMCAD. 

A . BACKGROUND 

MCADGRAF was developed and written by LT. E. M. 
Connolly using the programming language of Ryan-McFarland 
FORTRAN (Appendix M) . Microcompatibles GRAFMATIC (Appendix 
M) , and later Jewell Technologies GRAFPLUS (Appendix M) were 
chosen as the commercially available screen graphics and 
hard copy support for EMCAD output graphics. Any questions 
concerning colors, formats, changes to the program, etc. 
should be referred to the respective user's manuals, Dr. M. 
A. Morgan or LT. E. M. Connolly (Appendix M) . 

B. COMPUTER SETUP 

MCADGRAF was written for the following computer 
configuration : 

* IBM AT or compatible, including 80386 based 
compatibles 

* hard disk 

* EGA graphics card 
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EGA monitor 



* 

* math coprocessor - not required 

* memory - 512K RAM minimum 

MCADGRAF can be placed in any directory, at the user's 
discretion, however, it is highly recommended that MCADGRAF 
and EMCAD be placed in the same directory. MCADGRAF input 
is taken directly from EMCAD output. This output is 
required to be in the same directory/sub-directory as 
MCADGRAF, or on the same floppy disk if using a floppy 
system instead of a hard disk. The only other requirement 
for initial setup deals with printer initialization. In 
order to obtain a printer graphical output, GRAFPLUS must be 
installed on the computer, configured for the specific 
printer of use, prior to running MCADGRAF. An EPSON LX- 
800 series printer was used when writing this code. This 
printer, as with any other printer to be used, was 
configured by inserting two lines in the AUTOEXEC.BAT file. 
The first of these lines was a path command to tell the 
computer where to find the appropriate files for GRAFPLUS. 
Assuming all the files for GRAFPLUS are contained in a 
directory on the C disk entitled Grafplus, an appropriate 
path entry would be 

path c: \;c: \Grafplus 
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This setup will also work with other IBM graphics compatible 
printers such as the Okidata ML - 192/193. The second of 
these lines inserted actually established the initialization 
values for the printer to be used. The second line must be 
inserted anywhere in the AUTOEXEC.BAT after the path 
command. The printer initialization used for the LX-800 was 

GRAFFLUS = #115 

The equals sign (=) directs the printer to print the 
graphics horizontally. This means that the top of the image 
is at the top left corner of the page. For vertical 
printing, the object appears to be rotated with the top of 
the picture at the bottom left corner of the page. The 
pound sign (#) directs the printer to print the text screens 
as graphics screens. The I conveys the information that the 
screen mode is EGA - 640x350 resolution, with 16 colors 
available. The number 15 is the number assigned by GRAFPLUS 
to correspond with the EPSON LX-800. Other options 
available are the capabilities to print in color and to 
reverse printing. All of these options are discussed in 
much greater detail on pages 3 through 8 of the GRAFPLUS 
user's manual. As a note of caution, some GRAFPLUS printer 
definitions are defined such that "graying" is used on a 
black and white print-out to represent different colors. On 
the LX-800, the "graying" effect reduced printer clarity and 
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resolution dramatically. This problem was rectified by 
running the EDITGRAF.COM program contained in GRAFPLUS, 
answering all the questions without changes until the 
question "Do you wish to print in Black & White only? 
(disable gray scale) " appeared on the screen. Changing the 
answer from no to yes removed the "graying" and improved the 
resolution. A much more detailed explanation of EDITGRAF 
can be found on pages 9 through 21 of the GRAFPLUS user's 
manual. Once these changes to the printer definitions and 
AUTOEXEC.BAT have been made, it is necessary to re-boot the 
computer to initialize these values. This re-boot can be 
accomplished by either turning the machine off, waiting 15 
to 20 seconds, and turning it back on, or by hitting the 
Ctrl, Alt and Del keys all at the same time. 

C. PROGRAM INPUT 

As previously mentioned, MCADGRAF obtains its input from 
the output of EMCAD. This EMCAD output is four files in the 
form of "FILENAME" .TET, "FILENAME" .TEP, "FILENAME" . TMT and 
"FILENAME" . TMP. The quotation marks indicate that the 
FILENAME is a value that was previously chosen by the user 
while running EMCAD. The four extensions of FILENAME have 
significant meaning. The TET refers to TE incidence, 
F-Theta; the TEP refers to TE incidence, F-Phi ; the TMT 
refers to TM incidence, F-Theta; the TMP refers to TM 
incidence, F-Phi. The values for these four files were 
taken directly from a much larger, more comprehensive EMCAD 
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output data file. The scattering calculation output and 
graphics relationship is depicted in Figure H-l. 

D. PROGRAM EXECUTION 

MCADGRAF is a menu driven output graphics program. Once 
the correct disk or directory has been selected, MCADGRAF 
can be invoked by typing 



MCADGRAF 



at the DOS prompt. The program will load itself, and the 
user will see the following screens: 

* Welcome To MCADGRAF (Figure H-2) 

This screen presents some brief background information. 

* EMCAD Output Graphics Presentation, Part 1 

(Figure H-3) 

This menu affords the user the opportunity to exit the 
program in the event that the printer initialization has not 
been accomplished. 

* EMCAD Output Graphics Presentation, Part 2 

(Figure H-4) 

This screen explains the means of obtaining an output graph 
on the printer. The user will be presented with four graphs 
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SCATTERING CALCULATION 
OUTPUT AND GRAPHICS 



xxxxxxxxxx UELCOHE TO MCADGHAF xxxxxxxxxx 



THIS PROGRAM WAS URITTEH TO BE USED IN CONJUNCTION WITH 
THE FORTRAN PROGRAM ENTITLED ENCAD. MCADGRAF ADDS 
GRAPHING FEATURES TO SUPPORT THE ANALYSIS OF DATA 
GENERATED AND FORMATTED BY ENCAD. ENCAD UAS WRITTEN BY 
DR. H. A. MORGAN, WITH MODIFICATIONS BY HIMSELF AND 
LT E. H. CONNOLLY. MCADGRAF WAS DEVELOPED AND WRITTEN BY 
LT E. H. CONNOLLY TO HEP FULFILL THESIS REQUIREMENTS WHILE 
UNDER THE GUIDANCE AND THE SUPPORT OF DR. MORGAN. MCADGRAF 
WAS WRITTEN AND COMPILED USING RYAN MCFARLAND FORTRAN, 
GRAFMATIC, AND GRAFPLUS. ANY QUESTIONS CONCERNING COLORS, 
FORMAT, CHANGES TO THE PROGRAM, ETC... SHOULD BE REERRED 
TO THE RESPECTIVE USERS MANUALS, DR. M. A. MORGAN, OR 
LT E. M. CONNOEY. 

PLEASE PRESS ENTER TO CONTINUE. 



FIG. 



H-2. WELCOME TO MCADGRAF 



ENCAD OUTPUT GRAPHICS PRESENTATION 



NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, PART 1 

This progran allows the user to obtain printer hardcopies 
of the graphs being viewed. Hardcopy graphic outputs are 
obtained through the use of the progran entitled GRAFPLUS. 
Proper utilization of HCADGRAF requires prior setup of 
GRAFPLUS. Directions for proper setup are contained 
in the GRAFPLUS users nanual. 

1. CONTINUE WITH THE NEXT SCREEN 

2. EXIT THE PROGRAN 

Indicate your selection by the choice of 1 or 2. 



FIG, H-3. EMCAD OUTPUT 
GRAPHICS PRESENTATION PART 1 
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ENCAD OUTPUT GRAPHICS PRESENTATIOH 



NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, PART 2 

The user nakes the decision as to whether a printer 
graphical output is desired after viewing the graph. Once 
a graph is drawn on the screen, the user nay request a 
printer hardcopy by typing either a "P" or "p" and pressing 
the "ENTER" key. Selection of any other key will simply 
advance to the next graph. Spacing on 8 1/2 by 11 paper 
allows up to 2 graphs per page. 



PLEASE PRESS ENTER TO CONTINUE. 



FIG. H-4. EMCAD OUTPUT 
GRAPHICS PRESENTATION PART 2 
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for each data run. As each graph is presented on the 
screen, the user is given the opportunity to obtain a 
printer copy of the screen. The user decides to obtain a 
printer copy by depressing a P or a p, followed by ENTER. 
The reminder line telling the user what to do to produce a 
hardcopy is erased, and the graph is sent to the printer. 
Any other character input, followed by ENTER will simply 
advance to the next screen. Setting up the paper in the 
printer is the sole responsibility of the user for each and 
every graph. If one graph is to placed on every page, it is 
recommended that the paper be advanced a few inches so that 
the graph is printed closer to the center of the page. For 
the LX-800, it is possible to place 2 graphs on every page, 
if the paper is positioned properly in the printer. The 
best advice is to experiment. 

* Data File Name Input Menu (Figure H-5) 

If the data file to be used was generated by EMCAD, the user 
need only choose option 1. If the file to be used is a 
self-generated file, it must be in the exact format of an 
EMCAD generated file, and it must have an extension of .TET, 
.TEP, .TMT or .TMP . If the file is not in the correct 
format, or if the file does not contain the proper 
extension, the user must exit the program and correct the 
problem prior to continuing. 
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ENCAD OUTPUT GRAPHICS PRESENTATION 



DATA FILE NAME INPUT MENU 

If utilizing a file op series of files created by ENCAD, 
please select option 1. If utilizing any file other than 
one created by ENCAD, please ensure that the file is in 
proper format > and that the filename extension is one of 
the four extensions utilized by this progran. These 
extensions are: .tet .tep .tnt and .tnp 

If you are unsure as to whether your file is in the proper 
fornat with the proper extension, please exit the progran 
and verify before continuing. 

1. INPUT THE FILE OF USER CHOICE 

2. EXIT THE PROGRAN 

Indicate your selection by the choice of 1 or 2. 



FIG. H-5. DATA FILE NAME 
INPUT MENU 
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* Data File Name Input (Figure H-6) 

This menu prompts the user for the filename of the input 
data file. Since the extension is limited to four cases, 
and has previously been selected in the Data File Type 
Selection Menu, no extension should be added to the filename 
to be typed in. 

* Data File Type Selection Menu (Figure H-7) 

As previously discussed, MCADGRAF requires that the four 
EMCAD output data files be located in the same directory. 
This menu gives the user the chance to choose which of the 
four files the data will come from. 

This menu also affords the user the opportunity to change 
the default setting of the aspect ratio. The menu for this 
change is shown in Figure H-8 . Aspect ratio is a term that 
relates the number of available pixel columns on the screen 
to the number of available pixel rows on a screen. Since 
most computer screens are not square, this representation 
becomes very important. In simple terms, it deals with 
making a square look like a square, and not like a 
rectangle, or making a circle look like a circle, and not an 
oval . 

GRAFPLUS comes with a program to help determine the 
aspect ratio for any monitor. In the case of the 
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ENCAD OUTPUT GRAPHICS PRESENTATION 



DATA FILE NAME INPUT 

Please input the name of the file containing data that 
you wish to work with. Please input only the file nane 
i.e. do not include the file extension. The file 
extension will be selected in the DATA FILE TYPE 
SELECTION MENU. 



FIG. H-6. DATA FILE NAME INPUT 
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ENCAD OUIPUI GRAPHICS PRESENTATION 
DATA FILE TVPE SELECTION MENU 



Please select the type of data file containing the information of 
interest. ENCAD automatically created these four data file type s. 



1 . 


TE INCIDENCE, 


F-IHETA 


2. 


TE INCIDENCE, 


F-PHI 


3. 


TN INCIDENCE, 


F-THETA 


4. 


TN INCIDENCE, 


F-PHI 


5. 


CHANGE DATA FILE 




6 . 


CHANGE ASPECT RATIO 


7. 


EXIT THE PROGRAM 





Indicate your selection by the choice of 1, 2, 3, 4, 5, 6 or 7 



FIG. H-7. DATA FILE TYPE 
SELECTION MENU 
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XXXXXXXXXX ASPECT RATIO XXXXXXXXXX 



IN ORDER TO AUOID DISTORTION ON IHE SCREEN, IT IS NECESSARY 
10 DEFINE THE NUMBER OF ROWS PER INCH AND THE NUMBER OF 
COLUMNS PER INCH. IHESE DEFINITIONS CHANGE AS IHE 
RESOLUTION OF THE SCREEN, AS WELL AS THE TVPE OF SCREEN 
CHANGE. THIS PROGRAM WAS WRITTEN IN MODE 16 (EGA) USING 
A GB-1 UIDEO BOARD AND A NEC MULTISYNC MONITOR. THE 
OPTIMUM ASPECT RATIO FOR THIS CONFIGURATION WAS .65 



PLEASE INPUT ASPECT RATIO. 



FIG. H-8. ASPECT RATIO 
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development of this program, the aspect ratio given by 
GRAFPLUS delivered erroneous results. As such, this screen 
of MCADGRAF allows the user to input and change the aspect 
ratio of choice. For users unsure of their aspect ratio, 
start with an initial value of 1, and work up or down with 
increments of .1 or .2 until close to the desired results, 
then fine tune as much as needed. 

* Selection of Incidence Angle (Figure H-9) 

One of the input parameters for EMCAD was the number of 
incident fields, and their respective angles. MCADGRAF 
reads the input data file, and presents in this menu the 
number of incidence angles, and their respective angles in 
degrees. It is up to the user to select the incidence angle 
of interest. 

* Selection of Phi Angle (Figure H-10) 

As with incidence angles, phi angles were input parameters 
for EMCAD. MCADGRAF presents in this menu the number of phi 
angles, and their respective angles in degrees. It is up to 
the user to select the phi angle of interest. 
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ENCAD OUTPUT GRAPHICS PRESENTATION 



SELECTION OF INCIDENCE ANGLE 
The data file you have selected contains data for 5 
incidence angles. Please select your choice fron the 
following menu. Note, you nay only uorh with one 
incidence angle at a tine. 



1 INCIDENCE ANGLE = 

2 INCIDENCE ANGLE = 

3 INCIDENCE ANGLE = 

4 INCIDENCE ANGLE = 

5 INCIDENCE ANGLE = 

6 RETURN TO THE DATA 

7 EXIT THE PROGRAN 



0.00 

45.00 

90.00 

135.00 

180.00 

FILE TYPE SELECTION MENU 



INDICATE YOUR SELECTION BY THE CHOICE OF A NUMBER I.E. 1 2 3 



FIG. H-9. SELECTION OF 
INCIDENCE ANGLE 
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ENCAD OUTPUT GRAPHICS PRESENTATION 



SELECTION OF PHI ANGLE 

The data file you have selected contains data for 3 
phi angles. Please select your choice from the following 
menu. Note, you nay only work with one phi angle at a tine. 

1 PHI ANGLE = 0.00 

2 PHI ANGLE = 90,00 

3 PHI ANGLE = 180.00 

4 RETURN TO THE DATA FILE TYPE SELECTION NENU 

5 EXIT THE PROGRAM 

INDICATE YOUR SELECTION BY THE CHOICE OF A NUMBER I.E. 1 2 3 ,,.) 



FIG. H-1D. SELECTION OF 
PHI ANGLE 
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E . PROGRAM OUTPUT 



For each separate data set, MCADGRAF produces four 
different graphs. 

* Magnitude Versus Theta (Figure H-ll) 

As can be seen in the first line of the title for this 
graph, as well as the other three graphs, a label specific 
to this data set appears. This label was used as an input 
parameter to EMCAD to assist the user in differentiating 
this set of data from every other set. In the third line, 
the values of the incidence angle and phi angle, as well as 
the type of incidence are listed. Either F-Phi or F-Theta 
appears in the fourth line followed by a short description 
of the type of graph. The line at the bottom referring to 
hardcopy is simply a reminder to the user of what is needed 
to request a hardcopy output. This line is erased before 
the printer graph is produced. This first graph depicts the 
scattering magnitude plotted versus the 180 degrees of 
Theta . 

* Phase Versus Theta (Figure H-12) 

The only change in this graph from the first graph is that 
phase instead of magnitude is being graphed versus Theta. 
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Offset Sphere Validation Run 
Scattered Field Conplex Anpl itudes 
IE Incidence = 90. Degrees Phi = 0. Degrees 

F-Phi Magnitude Versus Theta 

Magnitude 




Hardcopy — > P or p 



FIG. H-ll. MAGNITUDE VERSUS 
THETA 
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Composite Body 

Scattered Field Conplex Amplitudes 
Ifl Incidence = 135. Degrees Phi = 180. Degrees 
Phase F-Theta Phase Versus Iheta 




glheta 



FIG. H-12. PHASE VERSUS 
THETA 
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* Magnitude Versus Theta, Polar Plot 

(Figure H-13) 

This graph represents the same data as the first graph, but 
from a different viewpoint. The first graph was drawn using 
rectangular coordinates. This graph utilizes polar 
coordinates. This graph type is the one predominately used 
in the industry to illustrate scattering. 

* Magnitude Versus Theta, Polar Plot 
Axis Rotated Minus Ninety Degrees 

(Figure H-14) 

The only change in this graph from the third graph is the 
rotation of the axis. The limitations of GRAFMATIC made it 
difficult to generate the data, maintain spacing 
restrictions and still display all the data necessary to 
adequately describe the situation. Rotating the axis ninety 
degrees displays the data slightly differently from the 
third graph, giving the user a bit more useful information. 
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Offset Sphere Validation Run 
Scattered Field Conplex Amplitudes 
IE Incidence = 90. Degrees Phi = 0. Degrees 

Polar Plot, F-Phi Magnitude Versus Theta 




FIG. H-13. MAGNITUDE VERSUS 
THETA, POLAR PLOT 
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Conposite Body 

Scattered Field Conplex Amplitudes 
III Incidence = 135. Degrees Phi = 180 . Degrees 
Polar Plot, F- Theta Magnitude Versus Theta 
Axis Rotated Minus Ninety Degrees 



S Theta-90 




Theta = 0 



Hardcopy — > P or p 



Theta=180 



FIG. H-14. MAGNITUDE VERSUS 
THETA, POLAR PLOT, AXIS 
ROTATED MINUS NINETY DEGREES 
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APPENDIX I 



MCADGRAF SOURCE CODE 



PROGRAM MCADGRAF. FOR 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THIS PROGRAM WAS WRITTEN TO BE USED IN CONJUNCTION WITH THE 
FORTRAN PROGRAM ENTITLED EMCAD. MCADGRAF ADDS GRAPHING FEATURES 
TO SUPPORT THE ANALYSIS OF DATA GENERATED AND FORMATTED BY EMCAD. 
EMCAD WAS WRITTEN BY DR. M. A. MORGAN, WITH SUBSEQUENT UPDATE BY 
HIMSELF AND LT E. M. CONNOLLY WHILE AT THE NAVAL POSTGRADUATE 
SCHOOL (MODS MAR 87 - MAR 88). MCADGRAF WAS DEVELOPED AND WRITTEN 
BY LT E. M. CONNOLLY TO HELP FULFILL THESIS REQUIREMENTS WHILE 
UNDER THE GUIDANCE AND THE SUPPORT OF DR. MORGAN. MCADGRAF WAS 
WRITTEN AND COMPILED USING RYAN MCFARLAND FORTRAN, GRAFMAT IC , AND 
GRAFPLUS. ANY QUESTIONS CONCERNING COLORS, FORMAT, CHANGES TO THE 
PROGRAM, ETC... SHOULD BE REFERRED TO THE RESPECTIVE USERS 
MANUALS, DR. M. A- MORGAN, OR LT E. M. CONNOLLY. 

A NOTE OF CAUTION, TO OBTAIN CORRECT PRINTOUTS FROM MCADGRAF, 
GRAFPLUS MUST BE PROPERLY CONFIGURED TO SUPPORT YOUR HARDWARE. 
QUESTIONS CONCERNING THE PROPER HARDWARE SETUP SHOULD BE REFERRED 
TO THE GRAFPLUS USER’S MANUAL. 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
PROGRAM STATEMENT DEFINITIONS 

XX xxxxx xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



QSMODE (mode) 



Clears 


the screen and sets the 


mode 


of the CRT 


display. May be used in 


either 


the 


text or the 


graphics mode. The desired value of 


mode is 


chosen 


from the following 


table 


: 




MODE 


RESOLUTION (CxR) 




COLOR T/G 


0 


40 X 25 




B/W 


T 


1 


40 X 25 




C 


T 


2 


80 X 25 




B/W 


T 


3 


80 X 25 




C 


T 


A 


320 X 200 




C 


G 


5 


320 X 200 




B/W 


G 


6 


640 X 200 




B/W 


G 


7 


80 X 25 




B/W 


T(M) 


13* 


320 X 200 




C 


G 


14* 


640 X 200 




c 


G 


15* 


640 X 350 




50HZ 


MONO 


16* 


640 X 350 




4/16 


COLORS 


x - — — 


> Requires EGA card 








jcoll, j c o 1 2 , 


jrowl, jrow2, xmin. 


xmax, 


ymi n 


, ymax. 


xorg, yorg, iopt, yoverx, aspect) 
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j co 1 1 
j co 1 2 
j rowl 
j r ow2 
xmin 

xmax 

ymin 

ymax 



xorg 

yorg 

i opt 
yoverx 

aspect 



Plot initialization routine; sets plot window, 
correspondi ng x and y ranges, and scaling 
factors . 

left hand pixel column defining plot region 
right hand pixel column defining plot region 
lower pixel row defining plot region 
upper pixel row defining plot region 
value of x variable assigned to correspond to 
column jcoll 

value of x variable assigned to correspond to 
column jcol2 

value of y variable assigned to correspond to 
row jrowl 

If iopt - 0, then ymax is the value of the y 
variable assigned to row = jrow2. If iopt = 1, 
then on input ymax is the assumed maximum value 
of the y variable required by the plot, 
vertical (y) axis intersects the horizontal Cx) 
axis at x = xorg 

horizontal Cx) axis intersects the wertical (y) 
axis at y = yorg 

sets scaling method option, see ymax 
number of y unit changes for a unit change in 
value of x, ie.. the change in the y coordinate 
per inch divided by the change in the x 
coordinate per inch 

aspect ratio for the user's CRT, for NEC 
Multisync, aspect = 0.6 for mode = 16 



QXAXIS ( xst, xfin, xmajor, minor, label. 



ndec ) 



Draws the x axis along the horizontal direction 
including tic marks and tic mark labels 

xst value of x at the start of the x axis; 

If a logarithmic scale is being used, then 
xst and xfin must be one of the following 
with xfin > xst : 

0 . 001 , 0 . 01 , 1 . 0 , 10 . 0 , 100 . 0 , 1000 . 0 , 10000.0 

xfin value of x at the end of the x axis 

xmajor distance between the major (large) tic marks; 

If xmajor is positive/negative, then the tic 
marks will be placed above/below the x axis 
respecti vely . If xmajor is zero, then no tic 
marks are inserted 

minor integer flag for use of minor tic marks or 

logarithmic scale; If minor = 0, then no tic 
marks are inserted. If minor = +1, then minor 
tic marks are inserted halfway between the 
major tic marks, but are not labeled. If 
minor = -1, then a logarithmic set of tic marks 
are inserted. 

label integer flag for auto labeling of tic marks; 

If label = 0, then there is no auto labeling. 

If label = 1, then the major tic marks are 
automatically supplied with a label. If label 
= “1, then the origin is labeled as well, 
ndec number of decimal places for auto labeling; 

The format of the label is FlO.n with the 
decimal point lined up with the major tic marks 



QYAXIS ( yst, yfin, ymajor, minor, label, ndec ) 
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Creates the y axis for a plot. 



yst 
yf i n 
yma j or 

minor 

label 

ndec 

QPREG ( 

ireg 
i va 1 

QPTXT ( 

ncha r 
char 

kolor 
i col 
i row 

QPTXTB 
QPTXT C 



QCMOV ( 

QSETUP 

ndots 

lkolor 

isymbl 



value of y at the start of the y axis 
value of y at the end of the y axis 
distance between major (large) tic marks on 
the y axis 

integer flag for use of minor tic marks or 
logarithmic scale; see QXAXIS 

integer flag for auto labeling of tic marks; 
see QXAXIS 

number of decimal places for auto labeling; 
see QXAXIS; 0 <= ndec <= 3 

ireg, ival) 

Changes the color table by modifying the 
Pallette Registers; EGA card only 

pallette register index (0~X5) 
value to be placed in register (0-63) 



nchar, char; kolor; icol; irow ) 

Writes plot text. The character variable char 
of length nchar is written at ( column; row ) = 

( icol; irow ) in specified color. 

length of character variable 

character variable corresponding to desired 

text 

specifies color of the text 
starting text column of text 
text row containing text 

( nchar; char; kolor ) 

Writes plot text at position B. A character 
variable char of length nchar is written at 
position ("B" ) just beyond the positive X axis 

( nchar, char, kolor ) 

Writes plot text at position C. A character 
variable char of length nchar is written at 
position ( ,f C ,f ) just above and centered on the 
Y axis 

icol, irow ) 

Moves the cursor to a specified column and 
row 

( ndots, lkolor, isymbl, klrsym ) 

Sets up line type (solid or dashed) and symbols 
to be used by QTABLE when plotting data 

Number of pixels (dots) specifying length of the 
dash. If ndots = 0, then a solid line is drawn. 
The input ndots = 1 is not allowed. By varying 
ndots > = 2, you may create a line composed of 
variable length dashes 
Specifies color of the line 

Specifies symbol to be inserted at data points 
QTABL. If isymbl = -1, then a cross is drawn at 
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k 1 rsym 


the data points. If isymbl = -2; then no symbol 
is drawn. Otherwise, the character with the 
ascii number "isymbl” is drawn 
Specifies color of the symbol 


QTABL C itype. 


npt, xarray, yarray ) 

Plots a tabulated function 


1 type 


Specifies point or line plot. When itype = 0, 
then the plot is a series of symbols at the 
tabulated points. When itype = 1, the symbols 
are joined by straight lines 


npt 

xarray 
ya r ray 
selectl 


The number of points in the function table 
X values in the function table 
Y values in the function table 
A character string used to hold the input 
values of incidence angle and phi angle 
selection 


seal exy 
xymax 


A scale factor used in three of the four graphs 
The maximum scaled value out of all the scaled 
X values and all of the scaled Y values 


f namel 


A character string used to hold the value of 
the input data file name without extension 


f nameZ 


A character string used to hold the value of 
the input data file name with extension 


xt ensi on 


A character string used to hold the value of 
the input data filename extension, i.e. 

.TET, .TEP, . TMT , or . TMP 


gra f print 


A character string used to allow the user to 
print out a graph from the screen 


dummy 

titlel 

title2 

ti tle3 

titled 

xarray 


An input variable used to advance the screen 
The first line of the title for each graph 
The second line of the title for each graph 
The third line of the title for each graph 
The fourth line of the title for each graph 
Data files contain numbers in THET, MAG, PHASE 
format. This polar coordinate format can be 
transformed to rectangular coordinate format. 
XARRAY is the array established to hold all of 
the transformed X values. 


ya r ray 


Data files contain numbers in THET, MAG, PHASE 
format. This polar coordinate format can be 
transformed to rectangular coordinate format. 
YARRAY is the array established to hold all of 
the transformed Y values. 


xamax 

yamax 

gamax 


The maximum value of the X array (XARRAY) 
The maximum value of the Y array (YARRAY) 
The maximum value of the X array (XARRAY) 
and the Y array (YARRAY) 


pi 

f name 


A constant approximately equal to 22/7 
A character string which equals either F-Phi or 
F-Theta 


numsets 


A number used for indexing purposes, equal to 
number of incidence angles times the number of 
phi angles 


numpts 


The number of theta angles per given incidence 
angle and phi angle 


a , b , c 


Dummy variables used for incrementation, 

a Index for incidence angles 

b Index for phi angles 

c Index for theta angles 


i ncangi 


A dummy variable used to read and hold the 
value of the incidence angle presently being 
uti 1 ized 
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incangl (a ) 
i 

phii 

incangld 

thetai 

ma gi 
phase! 

title6 
title? 
ti tle8 
tarrayC c ) 

rmarrayC c ) 

parrayC c) 

calcvar 

seal efac 
dummyvar 

theta (a, b, c) 
mag(a,b,c) 
phase( a , b , c ) 
radC c ) 



A matrix containing the values of all the 
incidence angles 

A dummy integer used as a counter 
A dummy variable used to read and hold the 
value of the phi angle presently being u ti ized 
A dummy variable containing an incidence a.gle 
value used for comparison purposes 
A dummy variable used to read and hold the 
value of the theta angle presently being 
uti 1 ized 

The magnitude of the output for the existing 
incidence angle, phi angle, and theta angle 
A dummy variable used to read and hold the 
value of the phase angle presently being 

utilized 

A dummy character string used to hold an input 
title to be used in an output graph 
A dummy character string used to hold an input 
title to be used in an output graph 
A dummy character string used to hold an input 
title to be used in an output graph 
A matrix containing the values of all the theta 
angles in degrees for a given incidence angle 
and phi angle 

A matrix containing the values of all the 
magnitudes for a given incidence angle 
and phi angle 

A matrix containing the values of all the 
phase angles for a given incidence angle 
and phi angle 

A variable used in intermediate ca lculati oins 
to determine a 
A graphical sec 
A dummy variabJ 



A matrix c 
angles for 
phi angles 
A matrix c 
magnitudes 
all the ph 
A matrix c 
angles for 
phi angles 
A matrix c 
angles in 
and phi an 



graphica 


1 scale 


factor 


(sea 


1 efac ) 


le facto 


r 








e used a 


s a limi 


t in scale 


factor 
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VARIABLE DECLARATIONS 



INTEGER A, B, C, I, IREG, IVAL, SCAL EFAC 

INTEGERS ICOL, IOPT, IRON, ISYMBL, JC0L1, JC0L2, JR0W1 

INTEGERS JR0W2, KLRSYM, KOLOR, LABEL, LKOLOR, MINOR 

INTEGERX2 MODE, NCHAR , NDEC, SELECT1 

CHARACTERS SELECT, XTENSION 

CHARACTERS GRAFPRNT 

CHARACTERS DUMMY, FNAME1 , TITLE7 

CHARACTERS 2 FNAME2 

CHARACTERX20 CHAR 

CHARACTERS 4 TITLE1, TITLE2, TITLE3, TITLE4, TITLE5 
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CHARACTERX64 TITLE8, TITLE6 ( 5, 19 , 1 ) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCANGI 
REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS, NUMPTS 
REAL NUMSETS, RAD, THETA I 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX , YMIN, YORG, YOVERX, YST 

REAL AARRAY ( 37 ) , BARRAY C 37 ) , I NCANGL C 5 ) , PHI C 19 ) 

REAL MAGC 5,19,37), PHASEC 5 , 1 9 , 37 ) , TH ET A C 5 , 1 9 , 37 ) 

COMMON/ONE/ RMARRAY C 37 ) , SCALEFAC, TARRAY C 37 ) 

C0MM0N/TH0/ PARRAY ( 37 ) 

COMMON/THREE/ XARRAY ( 37 ) , YARRAY C 37 ) , SCALEXY, XYMAX 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

**x*xxxx**x* ******* ************* ******************** ******** ****** 
READ DATA FILES, LOAD ARRAYS, CALCULATE MAXIMUM VALUES 
x***** ******* ************ ************ ****** ****** ****** *********** 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
ESTABLISH INITIAL VALUES 
pi = 3.1416 
aspect = .65 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
ESTABLISH INITIAL SCREEN MODE AND COLORS 
mode = 3 

CALL qsmode(mode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WELCOME TO MCADGRAF 

writeC X, 103) 

read(*,100) dummy 

CALL qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
NOTES CONCERNING GRAPHICAL OUTPUT PART 1 
6 0 writeC X , 1 01 ) 

read(X,100) select 
CALL qsmodeCmode) 

ifCselect ,ne. f l f .and. select .ne. * 2 * ^ go to 60 
ifCselect .eq. f 2 f ) go to 1 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

NOTES CONCERNING GRAPHICAL OUTPUT PART 2 

writeC X , 1 02 ) 

readCX,]GQ) dummy 

CALL qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DATA FILE NAME INPUT MENU 
4 writeC x , 9 1 ) 

readCX,100) select 

ifCselect .ne. 'l'.and. select .ne. f 2 f ) go to 4 
ifCselect .eq. '2 1 ) go to 1 
CALL qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT DATA FILE NAME 
45 writeC X, 92) 

readCX,10G) fnamel 
CALL qsmodeCmode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DATA FILE TYPE SELECTION MENU 
44 writeC X , 9 0 ) 

readCX,100) select 
xtension = select 
CALL qsmodeCmode) 

ifCselect .ne. 'I 1 .and. select .ne. f 2 f .and. select .ne. 1 3 1 

1. and. select .ne. *4* .and. select .ne. f 5' .and. select .ne. ’6’ 

2. and. select .ne. 1 7 •) go to 44 
ifCselect .eq. *7*) go to 1 
ifCselect .eq. *6’) go to 52 
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ifCselect . eq . *5 f ) go to 45 

OPEN DATA FILE 
ifCxtension . eq . *1*) then 

fname2 = f namel// *. tet T 
else ifCxtension .eq. *2’) then 
fname2 = f namel// *. tep ' 
else ifCxtension . eq . T 3 f ) then 
friame2 = f namel// 1 . tmt 1 
else 

fname2 = f namel// *. tmp 1 
endi f 

CALL qsmodeCmode) 
openC4, file=fname2) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

READ DATA FILE 

readC4,100) titleS 

readC 4,100) title7 

readC4,86) numsets 

readC4, 86) numpts 

a = 1 

incangld = -1000. 
b = 1 

do 16 i = 1 f intC numsets) 
c" 1 

readC4,86) incangi 
ifCincangi ,ne. incangld) b = 1 
ifCincangi ,ne, incangld) a = a + 1 
ifCi .eq. 1) a = 1 
incanglCa) = incangi 
incangld = incangi 
readC4,86) phii 
phiCb) = phii 
readC 4, 100) title6Ca,b,c) 
do 17 c 5 1 , i ntC numpts ) 
readC4,88) thetai, magi, phasei 
thetaCa,b,c) = thetai 
magCa,b,c) = magi 
phaseCa,b,c) = phasei 
17 continue 
b = b + 1 
16 continue 
cl oseC 4 ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SELECT DATA FILE SUBSET 

C I . E . SELECT INCIDENCE ANGLE AND PHI ANGLE) 

5 writeC*,93) a 
do 25 i=l,a 

wri teC *, 94) i,incanglCi) 

25 continue 
writeC*,95) a + 1, a + 2 
wri teC x, 96 ) 
readC*,*) selectl 
CALL qsmodeCmode) 

ifCselectl .ne. 1 .and. selectl ,ne. 2 .and. selectl .ne. 3 
l.and. selectl .ne. 4 .and selectl .ne. 5 .and. selectl .ne. 6 
2 .and. selectl .ne. 7) go to 5 
ifCselectl .eq. a + 1) go to 44 
ifCselectl .eq. a+2) go to 1 
a = selectl 

6 writeC*,97) b~l 
do 26 i=l,b-l 

writeC*,98) i,phiCi) 

26 continue 
writeC*,99) b, b+1 
wri teC 96 ) 
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readC*,*) selectl 
CALL qsmodeCmode) 

CALL qpreg( 0,32) 
ifCselectl .eq. b) go to AA 
ifCselectl .eq. b+1) go to 1 
b = selectl 

****************************************************************** 
LOAD ARRAYS 
c = 1 

xamax = 0.0 
yamax = 0.0 
gamax = 0.0 

title2 = ti 1 1 e6 C a , b, c ) 
do 18 c = 1 , intC numpts ) 
tarray(c) = thetaCa,b,c) 
rmarrayCc) = magCa,b,c) 
parray(c) = phaseCa,b,c) 
rad = tarrayCc) * pi / 180 
xarray(c) “ rmarrayCc) * sin(rad) 
yarray(c) = rmarrayCc) * cosCrad) 

****************************************************************** 
CALCULATE MAXIMUM VALUES 

if CabsCxarrayC c) ) . gt . xamax) xamax = abs CxarrayC c ) ) 
if CabsCyarrayC c) ) . gt . yamax) yamax - absCyarrayC c) ) 

18 continue 
75 continue 

****************************************************************** 
****************************************************************** 
CALCULATE SCALE FACTORS 
gamax = absCxamax) 

if Cyamax . gt . xamax) gamax = absCyamax) 
if CabsCgamax - 0) .It. .000001) gamax = .000001 
ealevar = al ogl 0 C gamax) 
scalexy = intCcalcvar) 

if Cgamax .It. 1.0) scalexy = intCcalcvar) - 1 
scalefac = intCscalexy) 
dummyvar = gamax/ C 1 0**scal exy) 
xymax = intC dummyvar +1) 

****************************************************************** 

CLEAR SCREEN 

mode = 16 

CALL qsmodeCmode) 

****************************************************************** 
CALL rectmvtCtitle2,title7,title8,numpts,aspect) 
****************************************************************** 
CLEAR SCREEN 
mode - 16 
CALL qsmodeCmode) 

****************************************************************** 
CALL rectpvtCtitle2,title7, titleS, numpts ; aspect ) 
****************************************************************** 
CLEAR SCREEN 
mode - 16 
CALL qsmodeCmode) 

****************************************************************** 
CALL polarmvt C ti tle2, title7 , title8, numpts , aspect ) 
****************************************************************** 
CLEAR SCREEN 
mode = 16 
CALL qsmodeCmode) 

****************************************************************** 
CALL polarrotCtitle2,title7,title8, numpts , aspect ) 
****************************************************************** 
CLEAR SCREEN 
mode - 16 
CALL qsmodeCmode) 
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C 



C 

C 

C 



****************************************************************** 
goto 44 

****************************************************************** 



ASPECT RATIO 
52 wr i te(*, 104) 

read(*,*) aspect 
mode - 3 

CALL qsmode(mode) 
goto 44 

****************************************************************** 
1 continue 
CLEAR SCREEN 
mode - 3 

CALL qsmodeCmode) 

****************************************************************** 



80 forma t(6X, F5 . 0, 3X, Ell . 3,6X, F7 .2) 

86 forma t ( F5 , 0 ) 

88 format(F5.0,E11.3,F7.2) 

90 format(/5X, 1 EMCAD OUTPUT GRAPHICS PRESENTATION 1 , /5X, 



5s. f ,//8X, 



DATA FILE TYPE SELECTION MENU f ,//5X, 

Please select the type of data file containing the information of 
,/5X, 

interest. EMCAD automatically created these four data file type 



1. TE INCIDENCE, F-THETA 1 ,//8X, 

2. TE INCIDENCE, F-PHI f ,//8X, 

3. TM INCIDENCE, F-THETA 1 , //8X, 

4. TM INCIDENCE, F-PHI ',//8X, 

5. CHANGE DATA FILE f ,//8X, 

6. CHANGE ASPECT RATIO 1 , //8X, 

7. EXIT THE PROGRAM 1 , //5X, 

Indicate your selection by the choice of 1, 2, 3, 4, 5, 6 or 7’, 
5//5X) 

91 format(//5X, 'EMCAD OUTPUT GRAPHICS PRESENTAT ION 1 , //5X, 

'DATA FILE NAME INPUT MENU',//5X, 

If utilizing a file or series of files created by EMCAD, 1 , /5X, 
please select option 1. If utilizing any file other than',/5X, 
one created by EMCAD, please ensure that the file is in',/5X, 
proper format, and that the filename extension is one of',/5X, 
the four extensions utilized by this program. These 1 , /5X, 
extensions are: . tet . tep . tmt and .tmp',/5X, 

If you are unsure as to whether your file is in the proper',/5X, 
format with the proper extension, please exit the program 1 ,/5X, 
and verify before continuing . 1 , //8X, 

1. ‘INPUT THE FILE OF USER CHOICE 1 , //8X , 

2. EXIT THE PROGRAM 1 ,//5X, 

Indicate your selection by the choice of 1 or 2.'//5X) 



92 format(//5X, 1 EMCAD OUTPUT GRAPHICS PRESENTATION 1 ,//5X , 



DA^A FILE NAME INPUT ',//5X, 

Please input the name of the file containing data that',//5X, 
you wish to work with. Please input only the file name f ,//5X, 
i.e. do not include the file extension. The file',//5X, 
extension will be selected in the DATA FILE TYPE',//5X, 
SELECTION MENU. f ,//5X) 



93 format(//5X, 1 EMCAD OUTPUT GRAPHICS PRESENTATION 1 , //5X , 
^SELECTION OF INCIDENCE ANGLE' ,/5X, 

2 f The data file you have selected contains data for ',I2,/5X, 
3'incidence angles. Please select your choice from the',/5X, 
4'following menu. Note, you may only work with one',/5X, 
5'incidence angle at a time. 1 ,/) 

94 format(5X, II, 5X, 'INCIDENCE ANGLE = f ,F6.2) 

95 formate 5X , 1 1 , 5X , 

1 ' RETURN TO THE DATA FILE TYPE SELECTION MENU ' ,/5X, II , 5X, 
2'EXIT THE PROGRAM') 

96 format(/5X, 'INDICATE YOUR SELECTION BY THE CHOICE OF A NUMBER 

1.1 2 3 . . . ) ' , //5X) 
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97 format(/5X, *EMCAD OUTPUT GRAPHICS PRESENTATION* ^//5X, 

1 * SEL ECTI ON OF PHI ANGLE*,/5X, 

2*The data file you have selected contains data for 1 ,I2,/5X, 

3 * phi angles. Please select your choice from the f ol 1 owi ng 1 , /5X , 
4*menu. Note, you may only work with one phi angle at a time.*,/) 

98 format(5X, 12, 5X, f PHI ANGLE = *,F6.2) 

99 format(5X, 12, 5X, 

1 1 RETURN TO THE DATA FILE TYPE SELECTION MENU * , /5X , 1 2 , 5X , 

2 * EXIT THE PROGRAM*) 

100 format(A) 

101 format(/5X, *EMCAD OUTPUT GRAPHICS PRESENTATI ON ' , //5X, 



102 format(/5X, *EMCAD OUTPUT GRAPHICS PRESENTATION * ,//5X, 



103 format(/5X, *xxxxxxxxxx WELCOME TO MCADGRAF XXXXXXXXXX * , //5X , 



NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, PART 1*,//5X, 

This program allows the user to obtain printer ha r dcopi es ' , / 5X , 
of the graphs being viewed. Hardcopy graphic outputs are*,/5X, 
obtained through the use of the program entitled GRAFPLUS . * ,/5X, 
Proper utilization of MCADGRAF requires prior setup of*,/5X, 
GRAFPLUS. Directions for proper setup are contained * ,/5X, 
in the GRAFPLUS users ma nua 1 , * , //5X , 

1, CONTINUE WITH THE NEXT SCREEN *, //5X, 

2 . EXIT THE PROGRAM* ,//5X, 

Indicate your selection by the choice of 1 or 2,*//5X) 



NOTES CONCERNING GRAPHICAL HARDCOPY OUTPUT, PART 2*,//5X, 

The user makes the decision as to whether a printer *,/5X, 
graphical output is desired after viewing the graph, 0nce',/5X, 
a graph is drawn on the screen, the user may request a *,/5X, 
printer hardcopy by typing either a **P n or "p** and pr essi ng * , /5X , 
the **ENTER** key. Selection of any other key will simply*, /5X, 
advance to the next graph. Spacing on 8 1/2 by 11 paper *,/5X, 
allows up to 2 graphs per page . * , ///5X, 

PLEASE PRESS ENTER TO CONTI NUE .*, 5X ) 



THIS PROGRAM WAS WRITTEN TO BE USED IN CONJUNCTION WITH * ,/5X, 

THE FORTRAN PROGRAM ENTITLED EMCAD. MCADGRAF ADDS *,/5X, 

GRAPHING FEATURES TO SUPPORT THE ANALYSIS OF DATA *,/5X, 

GENERATED AND FORMATTED BY EMCAD. EMCAD WAS WRITTEN BY*,/5X, 

DR . M. A. MORGAN, WITH MODIFICATIONS BY HIMSELF AND f ,/5X, 

LT E. M. CONNOLLY. MCADGRAF WAS DEVELOPED AND WRITTEN BY*,/5X, 

LT E. M. CONNOLLY TO HELP FULFILL THESIS REQUIREMENTS WHILE*, /5X, 
UNDER THE GUIDANCE AND THE SUPPORT OF DR. MORGAN. MCADGRAF *, /5X , 
WAS WRITTEN AND COMPILED USING RYAN MCFARLAND FORTRAN ,*, /5X , 
GRAFMAT IC , AND GRAFPLUS. ANY QUESTIONS CONCERNING COLORS ,*, /5X, 
FORMAT, CHANGES TO THE PROGRAM, ETC... SHOULD BE REFERRED * ,/5X, 
TO THE RESPECTIVE USERS MANUALS, DR. M. A . MORGAN, 0R*,/5X, 

LT E. M. CONNOLLY. * ,//5X, 

PLEASE PRESS ENTER TO CONTINUE .*, 5X) 



104 format(/5X, *XXXXXXXXXX ASPECT RATIO XXXXXXXXXX » , //5X, 

1 * I N ORDER TO AVOID DISTORTION ON THE SCREEN, IT IS NECESSARY *, /5X, 

2 * TO DEFINE THE NUMBER OF ROWS PER INCH AND THE NUMBER 0F*,/5X, 

3 * COLUMNS PER INCH. THESE DEFINITIONS CHANGE AS THE*,/5X, 

4 * RESOLUTION OF THE SCREEN, AS WELL AS THE TYPE OF SCREEN*, /5X, 

5 * CHANGE . THIS PROGRAM WAS WRITTEN IN MODE 16 (EGA) USING*, /5X, 

6 * A GB-1 VIDEO BOARD AND A NEC MULTISYNC MONITOR. THE*,/5X, 

7 * OPTIMUM ASPECT RATIO FOR THIS CONFIGURATION WAS .65*,///5X, 

8 * PL EASE INPUT ASPECT RATIO ,//5X) 

105 format(F5.3) 
end 



SUBROUTINE RECTMVT ( ti tl e2 , ti tl e7 , ti tl e8 , numpts, aspect ) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx xxxxx 

X XXXXXXXXXX X XX XXX XX XXXXX XX XXX XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XX 

MAGNITUDE VERSUS THETA 
RECTANGULAR COORDINATE GRAPH 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXX XX XXXXX XXXXXXXXXX XXXXX XXXXXXXXXX XXXXXXXXXX XXXXX XXXXXXXXXX XX XX 
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********************** ************ ********** ************ ********** 
VARIABLE DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER A, B , C, I, IREG, IVAL, SCALEFAC 

INTEGERS ICOL, IOPT, IRON, ISYMBL, JC0L1, JC0L2, JR0N1 

I NTEGER*2 JR0N2, KLRSYM, KOLOR, LABEL , LKOLOR, MINOR 

I NTEGER'*2 MODE, N C H A R , NDEC , SELECT1 

CHARACTERS SELECT, XTENSION 

CHARACT ER*3 GRAFPRNT 

CHARACTERS DUMMY, FNAME1 , TITLE7 

CHARACTERS 2 FNAME2 

CHARACTERSO CHAR 

CHARACTERS A TITLE1 , TITLE2, TITLE3, TITLED, TITLE5 
CHARACTERX6A TITLE3, TIT LE6 (5,19,1) 

REAL ASPECT, CALCVAR , DUMMYVAR, GAMAX, GMAX, INCANGI 
REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS, NUMPTS 
REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN , YORG, YOVERX, YST 

REAL AARRAY ( 37 ) , BARRAY ( 37 ) , INCANGL ( 5 ) , PHIC19) 

REAL MAGC5, 19,37), PHASEC 5 , 1 9 , 37 ) , THETAC 5 , 1 9 , 37 ) 

COMMON/ONE/ RMARRAY ( 37 ) , SCALEFAC, TARRAYC 37 ) 

COMMON/TMO/ PARRAY ( 37 ) 

COMMON/THREE/ XARRAYC37), YARRAYC37), SCALEXY, XYMAX 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SET DISPLAY MODE 

mode = 16 

CALL qsmode(mode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET AXIS COLOR 
CALL qpreg(7,6) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DEFINE GRAPH REGION 

jcoll = 030 

j col 2 = 520 

jrowl = AO 

jrow2 = 320 

xmin = -0.5 

xmax = 180 

ymin = -0.5 

ymax = xymax 

xorg = 00 

yorg = 0 . 0 

iopt - 1 

yoverx = .7X(ymax-ymin)/(xmax-xmin) 

CALL qplot(jcoll,jcol2,j rowl , j row2 , xmi n , xmax, ymin , ymax, xorg, 

1 yorg, iopt, yoverx, aspect ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SET UP LINE TYPE 

ndots = 0 

lkolor = 3 

isymbl - -2 

klrsym = A 

CALL qsetupC ndots, lkolor, isymbl , kl rsym) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DRAN X AXIS 

xst =0.0 

xfin = 130 

xmajor = A5 

minor = 0 

label = -1 

ndec = 1 

CALL qxaxisCxst , xf in , xmajor, mi nor, label, ndec) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
DRAN Y AXIS 
yst =0.0 
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yfin = xymax 
ynajor = 2 

if (xymax . eq . 1 .or, xymax . eq. 2 .or. xymax . eq . 3 .or. xymax 

l.eq. 4 .or. xymax . eq . 5) ymajor = 1 
minor = 1 
label = -1 

CALL qyaxisC yst , yf i n , ymajor, minor, label, ndec) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 
nchar - 20 
cha r = T Theta f 
kolor = 14 

CALL qptxtbC ncha r , char, kolor) 

XXXXXXXX XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 
nchar = 20 
char = Magnitude 1 
kolor = 14 

CALL qptxtcC nchar , char, kolor) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 
nchar = 64 
kolor = 10 
i col = 15 
irow = 24 

CALL qptxt ( ncha r , ti tle8, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE TEXT ON SCREEN 
nchar = 64 

titlel = ’Scattered Field Complex Amplitudes’ 
kolor = 10 
icol = 15 
irow = 23 

CALL qptxt (nchar, titlel, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE TEXT ON SCREEN 

nchar = 64 

kolor = 10 

icol = 15 

irow = 22 

CALL qptxt (nchar, title2, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE TEXT ON SCREEN 

nchar = 7 

kolor - 10 

icol = 15 

irow = 21 

CALL qp tx t( nc ha r,title7, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE TEXT ON SCREEN 
nchar = 40 

title3 = 1 Magnitude Versus Theta’ 
kolor = 10 
icol = 30 
irow = 21 

CALL qp txt(nc ha r,title3, kolor, icol, i row) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SCALE IT Y ,f VALUES 
do 30 i = l,numpts 

aarray(i) = rma r ray( i ) / ( 1 Q*Xsca 1 exy ) 

30 continue 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
GRAPH DATA POINTS 
itype = 1 
npt = numpts 

CALL qtabl (itype, npt, tarray,aarray) 
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WRITE SCALE FACTOR ON SCREEN 

nchar = 1 

ti t 1 e4 = 'S' 

kolor = 14 

i co 1 - 1 

irow “ 13 

CALL- qptxt ( nchar, ti tle4,ko lor, i col, i row) 
title4 = 'c' 
irow = 17 

CALL qptxt(nchar,title4,kolor,icol,i row) 
title4 " 'a' 
irow = 16 

CALL qptxt ( nchar, title4, kolor, icol, irow) 
ti tl e4 “ 
irow - 15 

CALL qptxt ( nchar , ti tle4, kolor , icol, i row) 
title4 = *e' 
irow r 14 

CALL qptxt (nchar, title4, kolor, icol, irow) 
ti tl e4 * ’X 1 
irow ' 12 

CALL qptxt ( nchar , ti tle4, kolor, icol, i row) 
nchar = 2 
ti tl e4 = f 1 0 f 
irow = 10 

CALL qptxt ( nchar, title4, kolor, icol, i row) 
nchar = 1 
ti t 1 e4 = 'X f 
irow - 9 

CALL qptxt (nchar, title4, kolor, icol, irow) 
ti tl e4 - '*' 
irow = 8 

CALL qptxt (nchar, title4, kolor, icol, irow) 

CALL qcmov(2,8) 
write(*,ll) scalefac 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX******** 

HARD COPY OF GRAPH ? 
nchar = 20 

titleS = 'Hardcopy -—> P or p* 
kolor = 14 
icol = 1 
irow = 1 

CALL qptxt ( nchar , title3, kolor, icol, irow) 

CALL qcmov (25,1) 
read( x, 1 00 ) graf prnt 
nchar - 40 

titleS = ' 1 

kolor = 10 
icol = 1 
irow = 1 

CALL qptxt ( nchar, title3,kolor, icol, irow) 

if(grafprnt . eq. 1 P ' .or. grafprnt , eq. 'p') CALL prtsc 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

11 f orma t ( 13) 

100 forma t ( A ) 
return 
end 

SUBROUTINE RECTPVT( title2, title7 , title8 , numpts, aspect) 

XX* XXXXXXXXX***** XX XXXXXXXX XX XXXXXXXXX XX XXXX XXXXXXXXX XXX* XXXXXXXXX 
XXXXXXXXXXXXXXXXXXXXXXXXXXXX ********** XXX X ************* XX xxxxxxxxx 

PHASE ANGLE VERSUS THETA 
RECTANGULAR COORDINATE GRAPH 

XXXXXXXXX XXXXXXXXX XXX XX XXX XXX ****** XXXXXXXXX XX XX XXX XXX XXX XX XXX XXX* 
** XX XXX XXX XXXXXXXXX********** ******* XXXXXXXXX* XXX XX XXX XX xxxxxxxxx* 
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x ********* ******** ****** ******* *********************** ****** ****** 
VARIABLE DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER A , B, C, I, IREG, IVAL, SCALEFAC 

INTEGERX2 ICOL, IOPT, IROW, I 5YMBL , JC0L1, JC0L2, JR0W1 

INTEGERS J RGW2 , KLRSYM, KOLOR, LABEL, LKOLOR, MINOR 

INTEGERS 2 MODE, NCHAR, NDEC , SELECT1 

CHARACTERS SELECT, XTENSION 

CHARACTERS GRAFPRNT 

CHARACTERS DUMMY, FNAME1 , TITLE7 

CHARACTER*12 FNAME2 

CHARACT ER*20 CHAR 

CHARACT ER*6 A TITLE1, TITLE2, TITLE3, TITLE4, TITLE5 
CHARACTERS^ TITLE8, T ITL E6 ( 5 , 1 9 , 1 ) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCANGI 
REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS, NUMPTS 
REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN, YORG, YOVERX, YST 

REAL AARRAYC37), BARRAY ( 37 ) , INCANGL ( 5 ) , PHI (19) 

REAL MAG( 5,19,37), PHASEC 5 , 1 9 , 37 ) , THETA ( 5 , 1 9 , 37 ) 

COMMON/ONE/ RMARRAY ( 37 ) , SCALEFAC, TARRAY ( 37 ) 

COMMON/TWO/ PARRAY ( 37 ) 

COMMON/THREE/ XARRAY ( 37 ) , YARRAY ( 37 ) , SCALEXY, XYMAX 
****************************************************************** 
SET DISPLAY MODE 
mode - 16 
CALL qsmode(mode) 

****************************************************************** 
SET AXIS COLOR 
CALL qpreg(7,6) 

****************************************************************** 

DEFINE GRAPH REGION 

jcoll = 080 

jco!2 = 530 

jrowl = 30 

jrow2 = 320 

xmin = -0.5 

xmax = 180 

ymin = -180 

ymax = 180 

xorg = 0.0 

yorg = 0.0 

iopt =* 1 

yoverx = 1.3 

CALL qplot(jcoll,jcol2,j rowl , j row2, xmi n, xmax, ymin, ymax, xorg, 
lyorg, iopt, yoverx, aspect ) 

****************************************************************** 

SET UP LINE TYPE 

ndots = 0 

Ikolor = 3 

isymbl = -2 

klrsym = A 

CALL qsetupCndots, Ikolor, isymbl , kl rsym ) 
****************************************************************** 
DRAW X AXIS 
xst =0.0 
xfin = 180 
xmajor = A 5 
minor = 0 
label = 1 
ndec = 1 

CALL qxaxisC xst , xf i n , xma i or , minor , label, ndec) 
****************************************************************** 
DRAW Y AXIS 
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yst = -130 
yfin = 180 
ymajor = 90 
minor - 1 
label = -1 

CALL qyaxisCyst, yfin, yma j or, minor, label, ndec) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 
nchar = 20 
char = ' Theta ' 
kolor * 14 

CALL qptxtbC nchar, char, kolor) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 
nchar = 20 
char s f Phase 1 
kolor = 14 

CALL qptxtcCnchar , char, kolor) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE TEXT ON SCREEN 

nchar = 64 

kolor = 10 

icol = 15 

irow = 24 

CALL qptxt(nchar,title8,koIor,icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE TEXT ON SCREEN 
nchar - 64 

title! - 'Scattered Field Complex Amplitudes' 
kolor * 10 
icol - 15 
irow “ 23 

CALL qptxtCnc ha r,titlel, kolor, icol, i row) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE TEXT ON SCREEN 

nchar = 64 

kolor - 10 

icol = 15 

irow = 22 

CALL qptxtC nchar, title2, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WRITE TEXT ON SCREEN 

nchar - 7 

kolor - 10 

icol = 15 

irow = 21 

CALL qptxtC nchar, title7, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITE TEXT ON SCREEN 
nchar = 40 

title3 = ' Phase Versus Theta 1 
kolor = 10 
icol = 30 
irow = 21 

CALL qptxt(nchar,title3,kolor,icol,irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
GRAPH DATA POINTS 
itype = 1 
npt = numpts 

CALL qtablC itype, npt , tar ray, par ray ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
HARD COPY OF GRAPH ? 
nchar = 20 

title3 = 'Hardcopy > P or p' 

kolor = 14 
icol = 40 
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irow = 1 

CALL qptxt(nchar,title3,kolor,icol,irow) 

CALL qcmov (65,1) 
read(X,100) grafprnt 
nchar = AO 

t i t 1 e3 = 1 1 

ko 1 or - 10 
i c o 1 = 4 0 
irow = 1 

CALL qptxtCnchar, tit le3,kolor,icol, irow) 
ifCgrafprnt .eq. ’P* .or. grafprnt . eq . 1 p 1 ) CALL prtsc 

c XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

100 format(A) 
return 
end 

SUBROUTINE POL ARMVTC ti tl e2, ti tl e7 , tit 1 eS , numpts, aspect ) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
MAGNITUDE VERSUS THETA 
POLAR COORDINATE GRAPH 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
VARIABLE DECLARATIONS 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXkXXXXXXXXXXXXXXXXXX 
INTEGER A, B, C, I, IREG, IVAL, SCALEFAC 
INTEGERS ICOL, I0PT, IROW, ISYMBL , JC0L1, JC0L2, JR0W1 
INTEGERS JR0W2, KLRSYM, K0L0R, LABEL, LK0L0R, MINOR 
INTEGERS MODE, NCHAR, NDEC , SELECT! 

CHARACT ERS SELECT, XTENSION 
CHARACTERS GRAFPRNT 
CHARACTERS DUMMY, FNAME1 , TITLE7 
CHARACTERS2 FNAME2 
CHARACTERSO CHAR 

CHARACTERX64 TITLE1 , TITLE2, TITLE3, TITLE4, TITLES 
CHARACTERX64 TITLES, TITL E6 ( 5 , 1 9 , 1 ) 

REAL ASPECT, CALCVAR, DUMMYVAR, GAMAX, GMAX, INCANGI 
REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS, NUMPTS 
REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 
REAL YAMAX, YFIN, YMAJOR, YMAX, YMIN , YORG, YOVERX, YST 
REAL AARRAY ( 37 ) , BARRAYC 37 ) , INCANGL ( 5) , PHIC19) 

REAL MAGC5, 19,37), PHASEC 5 , 1 9 , 37 ) , THETAC 5 , 1 9 , 37 ) 

COMMON/ONE/ RMARRAY ( 37 ) , SCALEFAC, TARRAY ( 37 ) 

COMMON/TWO/ PARRAY ( 37 ) 

COMMON/THREE/ XARRAY ( 37 ) , YARRAY ( 37 ) , SCALEXY, XYMAX 
xxx*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET DISPLAY MODE 
mode = 16 
CALL qsmode(mode) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SET AXIS COLOR 
CALL qpreg(7,6) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SCALE "X" AND "Y" VALUES 
do 10 i = 1, numpts 

aarray(i) = xar ray ( i )/ ( 1 OXXsca 1 exy ) 
barray(i) = yar ray ( i )/ ( 1 0 **sca 1 exy ) 

10 continue 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



DEFINE 


GRAPH REGION 


j CO 1 1 = 


030 


jcol2 = 


400 


jrowl = 


10 
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jrow2 = 280 
xmin = -xymax/2 

xmax - xymax \ 

ymin = -xymax 

ymax = xymax 

x o r g - 0.0 

yorg - 0,0 

i o p t = 1 

yoverx - 0,55 

CALL qplot(jcoll,jcol2,jrowl,j row2, xmin, xmax ; ymin , ymax, xorg, 
lyorg,iopt, yoverx, aspect) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C GET UP LINE TYPE 

ndots = 0 
lkolor - 3 
isymbl = -2 
klrsym = 4 

CALL qsetupC ndots, lkolor , isymbl, klrsym) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C DRAW X AXIS 

xst = 0.0 
xfin s xymax 
xmajor = 2 

if (xymax . le. 5) xmajor = 1 
minor ~ 1 
label ^ 1 
ndec = 1 

CALL qxaxisCxst , xf in , xmajor, minor, label, ndec ) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C DRAW Y AXIS 

yst ~ -xymax 
yfin = xymax 
ymajor = xmajor 
label - 0 

CALL qyaxisCyst , yfin, ymajor , minor, label , ndec ) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar = 10 
char = 1 Theta = 90 1 
kolor = 14 

CALL qptxtbC nchar, char, kolor) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar - 10 
char = r Theta = 0 f 
kolor = 14 
icol = 10 
irow = 19 

CALL qptxtC nchar, char, kolor, icol, irow) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar = 10 
char = 'Theta^lSQ 1 
kolor = 14 
icol = 10 
irow = 1 

CALL qptxtC nchar, char, kolor, icol, irow) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar = 64 
kolor = 10 
icol = 15 
iroi* = 24 

CALL qptxtC nchar, title8, kolor, icol, irow) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 
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nchar = 64 

titlel = 'Scattered Field Complex Amplitudes' 
kolor = 10 
icol = 15 
irow - 23 

CALL qptxtC nchar, titlel, kolor, icol, irow) 

XXX*XXXXXXXXXXXXXX**XXXXXXXXXXXXXX*XXXXXXXXXXXXX*XX**XXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 

nchar = 64 

kolor - 10 

icol = 15 

irow = 22 

CALL qptxtC nchar, title2, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar = 11 

title3 = 'Polar Plot/' 
kolor = 10 
icol - 15 
irow - 21 

CALL qp txtCnc ha r,title3, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar = 7 
kolor = 10 
icol = 27 
irow = 21 

CALL qptxtC nchar, title7, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar - 40 

titleS = ' Magnitude Versus Theta' 
kolor = 10 
icol = 35 
irow = 21 

CALL qptxt(nchar,title3,kolor,icol,irow) 
****************************************************************** 
WRITE SCALE FACTOR ON SCREEN 
nchar = 1 
title4 = 'S' 
kolor = 14 
icol = 10 
irow - 17 

CALL qptxt(nchar,title4, kolor, icol, irow) 
title4 = 'c' 
irow = 16 

CALL qp txtCnc ha r,title4, kolor, icol, irow) 
title4 = 'a' 
irow = 15 

CALL qptxtC nc ha r,title4, kolor, icol, irow) 
title4 = ' 1 ' 
irow = 14 

CALL qptxtC nc ha r,title4, kolor, icol, irow) 
ti 1 1 e4 = ' e ' 
irow = 13 

CALL qptxtC nc ha r,title4, kolor, icol, irow) 
title4 = 'X' 
irow = 11 

CALL qptxtC nchar, title4, kolor, icol ,irow) 
nchar = 2 
ti t 1 e4 = '10' 
irow = 9 

CALL qptxt Cnc ha r , ti tl e4 , kolor, icol, irow) 
nchar = 1 
ti tl e4 = '*' 
irow = 8 
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CALL qptxt C nchar, titled, kolor, icol, irow) 
title* = 
irow = 7 

CALL qptxtCnchar, titled, kolor, icol, irow) 

CALL qcmov (5,7) 
write(*,12) scalefac 

See discussion of QCMOV for problem requiring this fix 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX****** 

GRAPH DATA POINTS 
itype = 1 
npt = numpts 

CALL qtabl ( i type, npt , aarray, bar ray ) 
****************************************************************** 
HARD COPY OF GRAPH ? 
nchar = 20 

title3 = ’Hardcopy > P or p* 

kolor = 14 
icol = 50 
irow = 1 

CALL qptxt(nchar,title3,kolor,icol,irow) 

CALL qcmov ( 7 5,1) 
readC^, 100) graf prnt 
nchar = 26 

title3 = 1 1 

kolor = 10 
icol - 50 
irow = 1 

CALL qptxtC nchar, title3, kolor, icol, irow) 
ifCgrafprnt .eq, f P f -or. grafprnt .eq, ’p 1 ) CALL prtsc 
****************************************************************** 
12 f ormat( 9x, 13) 

100 format(A) 
return 
end 

SUBROUTINE POL ARROTC titleZ, txtle7, title8, numpts, aspect) 
****************************************************************** 
****************************************************************** 
MAGNITUDE VERSUS THETA 
POLAR COORDINATE GRAPH 
AXIS ROTATED MINUS NINETY DEGREES 

***************** ************************************************* 
****************************************************************** 

****************************************************************** 
VARIABLE DECLARATIONS 

****************************************************************** 

INTEGER A, B, C, I, IREG, IVAL, SCALEFAC 

INTEGERX2 ICOL, I0PT, IROW, ISYMBL , JC0L1, JC0L2, JR0W1 

INTEGERX2 JR0W2 , KLRSYM, KOLOR, LABEL, LKOLOR, MINOR 

INTEGERX2 MODE, NCHAR, NDEC, SELECT1 

CHARACTERX2 SELECT, XTENSION 

CHARACTERX3 GRAFPRNT 

CHARACTER'S DUMMY, FNAME1 , TITLE7 

CHARACTERX12 FNAME2 

CHARACT ERX20 CHAR 

CHARACTERX6A TITLE1 , TITLE2, TITLE3, TITLE4, TITLE5 
CHARACTERX6A TITLES, TITLE6C5, 19,1) 

REAL ASPECT, CALCVAR, DUMMYVAR , GAMAX, GMAX, INCANGI 
REAL INCANGLD, PHASEI, PHII, PI, MAGI, NDOTS , NUMPTS 
REAL NUMSETS, RAD, THETAI 

REAL XAMAX, XFIN, XMAJOR, XMAX, XMIN, XORG, XST 

REAL YAMAX, YFIN , YMAJOR, YMAX, YMIN , YORG, YOVERX, YST 

REAL AARRAY ( 37 ) , BARRAY C 37 ) , INCANGL ( 5 ) , PHIC19) 

REAL MAGC5, 19,37), PHASEC 5 , 1 9 , 37 ) , THETA C 5 , 1 9 , 37 ) 

COMMON/ONE/ RMARRAY C 37 ) , SCALEFAC, T ARRAY C 37 ) 
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COMMON/TWO/ PARRAY C 37 ) 

COMMON/THREE/ XARRAY ( 37 ) , YARRAY ( 37 ) , SCALEXY, XYMAX 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C SET DISPLAY MODE 

node - 16 
CALL qsmode(mode) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C SET AXIS COLOR 

CALL qpregC 7 , 6 ) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C DEFINE GRAPH REGION 

jcoll = 080 
j col 2 = 530 
jrowl = 10 
jrow2 = 320 
xmin = -xymax 
xmax = xymax 
ymin = -xymax/ 4 
ymax = xymax 
xorg - 0,0 
yorg = 0,0 
iopt = 1 
yoverx = 0.5 

CALL qplot(jcoll,jcol2, jrowl, j row2, xmi n , xmax, ymin , ymax, xorg, 

1 yorg, iopt , yoverx, aspect ) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C SET UP LINE TYPE 

ndots = 0 
lkolor - 3 
isymbl = -2 
klrsym = A 

CALL qsetupCn dots, lkolor, isymbl , kl rsym) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C DRAW X AXIS 

xst - -xymax 
xfin = xymax 
xmajor = 2 

ifCxymax .le. 5) xmajor = 1 
minor = 1 
label = 0 
ndec = 1 

CALL qxaxisC xst , xfin, xmajor, minor, label, ndec) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C DRAW Y AXIS 

yst = 0.0 
yfin = xymax 
ymajor = xmajor 
label = 1 

CALL qyaxisCyst , yf in , ymajor, mi nor, label, ndec) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar = 10 
char = *Theta=180 f 
kolor = 14 
icol = 65 
irow = 2 

CALL qptxt( nchar, char, kolor, icol, irow) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar = 10 
char = , Theta=90 l 
kolor = 14 

CALL qptxtcC nchar, char, kolor) 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C WRITE TEXT ON SCREEN 

nchar = 10 
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char - 'Theta - O' 
kolor = 14 
icol = 6 
irow = 2 

CALL qptxtCnchar, char, kolor, icol, irow) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE TEXT ON SCREEN 

nchar = 6 A 

kolor - 10 

icol = 15 

irow = ?A 

CALL qptxt( nchar, titled, kolor, icol, irow) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx********** 
WRITE TEXT ON SCREEN 
nchar = 64 

titlel = 'Scattered Field Complex Amplitudes' 
kolor - 10 
icol - 15 
irow = 23 

CALL qptxt( nchar, titlel, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar = 64 
kolor - 10 
icol = 15 
irow = 22 

CALL qptxtC nchar, title2, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar 5 11 

titled = 'Polar Plot, ' 
kolor - 10 
icol = 15 
irow - 21 

CALL qptxtCnchar, title3, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar s 7 
kolor - 10 
icol = 27 
irow = 21 

CALL qptxtCnchar, title7, kolor, icol, irow) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar s 40 

titled = * Magnitude Versus Theta 1 
kolor = 10 
icol = 35 
irow = 21 

CALL qptxt Cnc ha r,title3, kolor, icol, i row ) 
****************************************************************** 
WRITE TEXT ON SCREEN 
nchar - 64 

title5 = 'Axis Rotated Minus Ninety Degrees' 
kolor = 10 
icol - 15 
irow = 20 

CALL qptxtCnchar, title5, kolor, icol, irow) 
****************************************************************** 
SCALE "X" AND ''Y» VALUES 
do 20 i = l,numpts 

aarrayC i ) = -yarrayC i )/C 10X*scalexy ) 
bar ray C i ) = xarrayC i )/ C 1 0**scal exy ) 

20 continue 

****************************************************************** 
GRAPH DATA POINTS 
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itype = 1 
npt = numpts 

CALL qtabl(itype,npt,aarray,barray) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITE SCALE FACTOR ON SCREEN 
nchar = 1 
titled = 'S' 
kolor = 14 
icol = 5 
irow = 16 

CALL qptxt (nchar, title*, kolor, icol, irow) 
titled = f c» 
irow = 15 

CALL qptxt ( nchar, title*, kolor, icol, irow) 
titled = *a f 
irow = 14 

CALL qp t xt (nc ha r,titleA, kolor, icol, irow) 
titled = f l f 
irow = 13 

CALL qptxt C nchar, titled, kolor, icol, irow) 
titled = 1 e 1 
irow - 12 

CALL qptxt(nchar,title*,kolor,icol,irow) 
titled = »X* 
i row = 10 

CALL qptxt ( nchar, title*, kolor, icol, i row) 
nchar = 2 
title* = '10 f 
irow = 8 

CALL qptxt ( nchar, title*, kolor, icol, irow) 
nchar = 1 
title* = 
irow = 7 

CALL qptxt ( nchar, title*, kolor, icol, irow) 
title* = 
irow “ 6 

CALL qptxt (ncha r, title*, kolor, icol, i row) 

Call QCMOV ( * , 6 ) 
write(*,13) scalefac 

See discussion of QCMOV for problem requiring this fix 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

HARD COPY OF GRAPH ? 
nchar = 20 

title3 = 'Hardcopy > P or p 1 

kolor = 1* 
icol = 30 
irow = 1 

CALL qptxt (nc ha r,title3, kolor, icol, i row) 

CALL qcmov (55,1) 
readC^, 100) grafprnt 
nchar = *0 

title3 = 1 1 

kolor = 10 
icol = 30 
irow = 1 

CALL qp txtCnc ha r,title3, kolor, icol, irow) 

ifCgrafprnt . eq. 'P' .or. grafprnt ,eq. 1 p 1 ) CALL prtsc 
C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
13 formats Ax, 13 ) 

100 formatCA) 
return 
end 
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APPENDIX J 



SPHERE SOURCE CODE 



PROGRAM SPHERE 

GENERATING THE RADII a Er, Ur AT ROWS I a 1+1 
FOR A LAYERED DIELECTRIC SPHERE HAVING 

Z-AXIS OFFSET/LAMBDA = D. Written By M. Morgan 2/14/88 

Modified By E. Connolly 4/3/88 



REAL RAD(5),R0(5) 

COMPLEX ER/ UR 

CHARACTERS 4 FNAME, HDR1 , HDR2 
PI =3 , 1415927 
WRI TEC * , 110) 

READC*, 100) FNAME 

OPENC 1, FI LE= FNAME, STATUS=’ UNKNOWN ’,ACCE$$=’ SEQUENTIAL * , 
1 F0RM= f FORMATTED’ ) 

WRITEC*, 111) 

READC*, 100) FNAME 

OPENC 2/ FI L E= FNAME, STATUS =’ UNKNOWN f ,ACCES$=’ SEQUENTIAL f , 
1 F0RM= 'FORMATTED 1 ) 

WRITEC*, 112) 

READC*, 100) HDR1 
WRITEC*, 113) 

READC *,100) HDR2 
WRITEC1/105) HDR1 
WRITEC1/105) HDR2 
WRITEC*, 114) 

READC *, *) NL 
WRITEC 1,101) NL 
DO 11 L=1,NL 

WRITEC * / * ) f Layer No: ’,L 
WRITEC*,*) 

WRITEC*,*) ’Enter Radi us/Lambda : 1 

READC*,*) ROCL) 

WRITEC*,*) ’Enter Real Er of This Layer: ’ 

READC*,*) ERR 

WRITEC*,*) ’Enter Imag Er of This Layer: ’ 

READC*,*) ERI 
ER=CMPLXC ERR, ERI ) 

UR=CMPLXC1.0,-l.E-6) 

WRI i EC 1 , 102) ER, UR 
11 CONTINUE 

WRITEC*,*) ’Enter Z~Axis Of fset/ Lambda : ’ 

READC*,*) D 

WRITEC*,*) ’Enter IBIG: ’ 

READC*,*) IBIG 
WR1TEC2, 104) I BIG, NL + 1 
RMIN=2.*PI*CR0C1) -D) 

RMAX=2.*PI*CR0CNL)+D) 

WRITEC 2, 102) RMIN,RMAX 
DT = PI/C IBIG-1 . 0) 

DO 33 1 = 1, IBIG 
T=CI-1 . 0)*DT 
DS2=C D*SINCT ) )**2 
C=C0SCT) 

DO 22 L = 1 , NL 

22 RADCL )=2. 0*PI*CD*C+SQRTCR0CL )**2-DS2)) 

WRITEC 2, 103) I , C RADCL ) , L =1 , NL ) 

100 FORMAT C A ) 
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101 FORMAT ( 15 ) 

102 FORMAT ( A ( El A 6)) 

103 F0RMAT(I5,5(E14.6)) 

104 F0RMATC2C 15) ) 

105 FORMAT (' 1 , A) 

110 FORMAT(//7X, 'TWO FILES ARE REQUIRED TO HOLD PROGRAM OUPUT f ,//7X, 
l'DATA. PLEASE ENTER THE NAME OF THE FILE TO BE USED TO 1 ,//7X, 

2 * HO L D THE PROGRAM MATERIAL PARAMETER DATA. PLEASE ',//7X, 

3' INCLUDE THE EXTENSION. I. E. FI L ENAME . DAT 1 , //7X ) 

111 F0RMATC//7X, 'PLEASE ENTER THE NAME OF THE FILE TO BE USED f ,//7X, 
1 ' TO HOLD THE PROGRAM OUTPUT DATA. PLEASE INCLUDE THE f ,//7X, 

2 * EXT ENSI ON . I. E. FI L ENAME DAT 1 , //7X ) 

112 F0RMAT(//7X, ’THE HEADERS ALLOW THE USER TO IDENTIFY THIS',//7X, 

1 ' SET OF DATA FROM ALL OTHER SETS. 1 , 

2//7X, 1 PLEASE ENTER HEADER #1 (64 CHARACTERS MAX)’,//7X) 

113 F0RMATC//7X, ’PLEASE ENTER HEADER #2 (64 CHARACTERS MAX) ! ,//7X) 

114 F0RMAT(//7X, 'ENTER THE NUMBER OF LAYERS ( . L E . 5)',//7X) 

33 CONTINUE 

CLOSE(l) 

CL OSE( 2) 

STOP 

END 
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APPENDIX K 



CONE SOURCE CODE 



PROGRAM CONE 

******* xx*xx*xxx*xxxxx*xx*xxx*xxxxxxxxx*xxxxxxxxx*xxxxxxxxxxxxxxxx 



GENERATING THE RADII & Er, Ur AT ROWS I 8 1+1 
FOR A LAYERED DIELECTRIC CONE 
Written By E « Connolly A/19/88 

****************************************************************** 
VARIABLE DECLARATIONS 

****************************************************************** 
REAL RADC1000,5), PI, DELTHETA, DEGTORAD 
REAL THETA, H(5), X(5), Y(5) 

REAL ALPHAK5), ALPHAC5), B, M, RL, GAMMA 
INTEGER I BIG, I 
COMPLEX ER, UR 

CHARACTERS A FNAME, HDR1 , HDR2 

****************************************************************** 
INITIAL VALUES 

****************************************************************** 

RMIN = 10000 

RMAX = -10000 

PI -5 . 1 A15927 

DEGTORAD = PI/180. 



****************************************************************** 
INPUT NECESSARY DATA 

****************************************************************** 
WRITEC * ,110) 

READC *,100) FNAME 

OPENC 1 , FI LE= FNAME, STATUS = 1 UNKNOWN ’ , ACCESS = 'SEQUENTIAL 1 , 

1 F0RM= f FORMATTED' ) 

WRITEC*, 111) 

READC *,100) FNAME 

OPENC 2 , FI LE = FNAME , STATUS = 1 UNKNOWN * , ACC ESS = F SEQUENTIAL 1 , 

1 FORM= ' FORMATTED ' ) 

WRITEC *,112) 

READC 100 ) HDR1 
WRITEC*, 113) 

READC*, 100) HDR2 
WRITEC 1 , 105) HDR1 
WRITEC 1 , 105) HDR2 
WRITEC*, 11A) 

READC*,*) NL 
WRITEC 1 , 101) NL 
DO 11 L=1,NL 
WRITEC*,*) ' 

WRITEC*,*) 

WRITEC*, 10) 

READC*,*) ALPHA1CL) 

ALPHACL) = CALPHA1CD/2. ) 

WRITEC*,*) 

WRI r EC * , 1 5 ) 

READC*,*) HCL) 

WRITEC*,*) 1 
READC*,*) ERR 
WRITEC*,*) 1 



Layer No: 



* DEGTORAD 



Enter Real Er of This Layer: 
Enter Imag Er of This Layer: 



1 9 A 



oo oooo ooooooo o o o oooo 



READC*,*) ERI 
ER=CMPLX(ERR,ERI) 

UR=CMPLX( 1 . 0,-1 . E-6) 

HR I TEC 1 , 102) ER , UR 
11 CONTINUE 

WRITEC*, X) T Enter IBIG: 1 

READ!*,*) IBIG 

xxxxxxxxxxxxxxxxxxxxxxxxxx*x*xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx*xxxx 

CALCULATIONS 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXXXXXXXXXXXXX 

DELTHETA = PI/CIBIG-l) 

DC) 22 1 = 1, IBIG 

THETA = (1-1) * DELTHETA 

DO 33 L = 1 , NL 

Y(L) = HCU/2 

X(L) = ( Y( L ) ) *TAN( ALPHA ( L ) ) 

CHECK FOR HORIZONTAL LINE 
GAMMA = ATANCXC L )/YC L ) ) 

IFC (PI-GAMMA-THETA) ,LT. .0001) THEN 
YFOUND = -Y( L ) 

XFOUND = YFOUND X TANC THETA ) 

RADC I , L ) = 2*PI*SQRT(XF0UND**2 + YF0UND**2) 

ELSE 

SLOPE OF LINE CONTAINING CORNERS OF CONE 
M = ( -H(L)/X(L) ) 

Y-INTERCEPT OF LINE CONTAINING CORNERS OF CONE 
B = Y(L) 

RECALL EQUATION OF ANY LINE IS Y = MX + B 

THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS INTERSECTS 
THE LINE DRAWN FROM THE ORIGIN THROUGH THE DESIRED DATA POINT 
AT THE DESIRED DATA POINT XFOUND AND YFOUND. THE EQUATION FOR 
THE LINE CONTAINING THE I-TH AND K-TH DATA POINTS MUST THEN 
ALSO WORK AT THE DESIRED DATA POINT. WRITING XFOUND AND 
YFOUND IN TERMS OF THEIR ANGLE THETA(N) , AND SUBSTITUTING 
YFOUND = B/ (1-(M* TANC THETA ) ) ) 

XFOUND = YFOUND X TANC THETA ) 

RADC I , L ) = 2*PI*SQRT(XF0UND**2 + YF0UND**2) 

ENDIF 

IFC RADC I , L ) .GT. RMAX) RMAX = RADC I , L ) 

IFC RADC I , L ) .LT. RMIN) RMIN = RADC I , L ) 

33 CONTINUE 
22 CONTINUE 

XXX XX XXX XX XX XXXXXXXXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXXXXXXXXXXXXX 

OUTPUT 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

WRITEC2, 104) IBIG, NL+1 
WRITEC 2, 102) RMIN, RMAX 
DO 44 I = 1 , 1 BIG 

WRITEC 2, 103) I, (RADC I, L), L=1,NL) 

44 CONTINUE 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

10 F0RMATC/3X, 1 Please enter the apex angle, ALPHA1, in degrees, of 1 
l/8X,*the offset cone. Note, apex angle is the angle from one side 1 
2/8X,'of the cone to the other side of the cone. *) 

15 F0RMATC/8X, * Please enter the height, H, of the V 
18X , 1 cone. ») 

100 FORMAT C A ) 

101 FORMAT ( 15 ) 

102 FORMAT C 4 C E14 . 6 ) ) 

103 FORMAT (15, 5( El 4. 6)) 

104 FORMAT C 2 C 15 ) ) 

105 FORMAT C 1 f , A) 

110 F0RMATC//7X, 'TWO FILES ARE REQUIRED TO HOLD PROGRAM 0UPUT»,//7X, 
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1 1 DAT A . PLEASE ENTER THE NAME OF THE FILE TO BE USED TO ',//7X, 

2 1 HOL D THE PROGRAM MATERIAL PARAMETER DATA. PLEASE f ,//7X, 

3 f INCLUDE THE EXTENSION. I. E. FILENAME.DAT 1 ,//7X) 

111 FORMATC//7X, ’PLEASE ENTER THE NAME OF THE FILE TO BE USED f ,//7X, 
1 1 TO HOLD THE PROGRAM OUTPUT DATA. PLEASE INCLUDE THE',//7X, 

2 1 EXTENSION . I. E. FI L ENAME . DAT ' , // 7X ) 

112 FORMATC//7X, 1 THE HEADERS ALLOW THE USER TO IDENTIFY THIS 1 ,//7X, 

1 1 SET OF DATA FROM ALL OTHER SETS. 1 , 

2//7X, 1 PL EASE ENTER HEADER #1 (64 CHARACTERS MAX) f ,//7X) 

113 FORMATC//7X, ’PLEASE ENTER HEADER #2 (64 CHARACTERS MAX)',//7X) 

114 F0RMAT(//7X, ’ENTER THE NUMBER OF LAYERS (.LE. 5)’,//7X) 

CLOSEtl ) 

CL0SEC2) 

STOP 

END 
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OOOO DODO OOOOOOOOOOO 



APPENDIX L 



CYLINDER SOURCE CODE 



PROGRAM CYLINDER 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

GENERATING THE RADII & Er, Ur AT ROWS I 8 1+1 
FOR A LAYERED DIELECTRIC CYLINDER 
Written By E, Connolly 4/19/88 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

VARIABLE DECLARATIONS 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL RADCX000,5), PI, DELTHETA, DEGTORAD 
REAL THETA, H(5), RC5), YC5) 

REAL RL, GAMMA 
INTEGER I BIG, I 
COMPLEX ER, UR 

CHARACTERS 4 FNAME, HDR1 , HDR2 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INITIAL VALUES 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

RMIN = 10000 

RMAX = -10000 

PI=3 . 1415927 

DEGTORAD = PI/180. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INPUT NECESSARY DATA 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
WRITECX, 110) 

READ( x ,100) FNAME 

OPENC 1 , FI L E = FNAME, STATUS= 1 UNKNOWN 1 , ACCESS= f SEQUENTIAL ■ , 

1 FORM= f FORMATTED 1 ) 

WRITEC x, 1 1 1 ) 

READC X, 100) FNAME 

OP EN( 2, FI L E= FNAME, STATUS= f UNKNOWN 1 ,ACCESS= f SEQUENTIAL 1 , 

1 F0RM= 1 FORMATTED ' ) 

WRITEC X, 112) 

READC X, 10 0) HDR1 
WRITECX, 113) 

READC x , 10 0) HDR2 
WRITEC 1,105) HDR1 
WRITEC 1,105) HDR2 
WRITECX, 114) 

READC x, x ) NL 
WRITEC 1,101) NL 
DO 11 L = 1 , NL 

WRITECX, X) 'Layer No: * ,L 
WRITEC x , x ) 

WRITECX, 15) 

READC x , x ) HCL) 

WRITECX, 16) 

READC x , x ) RCL) 

WRITECX, X) 'Enter Real Er of This Layer; 1 
READC X, X) ERR 

WRITECX, x) 'Enter Imag Er of This Layer: ' 

READC x , x ) ERI 
ER = CMPLXC ERR, ERI ) 
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oo oooo o o o noon 



UR=CMPLXC 1 . 0 , -1 E-6) 

WRITEC 1,102) ER, UR 
11 CONTINUE 

WRITEC*,*) ' Enter IBIG: f 

READC *, *) IBIG 

****************************************************************** 

CALCULATIONS 

****************************************************************** 
DELTHETA = PI/CIBIG-1) 

DO 22 1*1, IBIG 

THETA = C 1-13 * DELTHETA 

DO 33 L*1,NL 

YCL) * HCL)/2 

CHECK FOR HORIZONTAL LINE 

GAMMA = ATAN(R(L)/Y(L) ) 

I FC ( PI-GAMMA- THETA) .LT. .0001) THEN 
YFOUND = -YCL) 

XFOUND = YFOUND x TANC THETA ) 

RADCI,L) = 2*PI*$QRTCXF0UND**2 + YFQUND**2) 

CHECK FOR HORIZONTAL LINE 

ELSE IFC C GAMMA-THETA) .GT. 0) THEN 

YFOUND = YCL) 

XFOUND = YFOUND x TANC THETA ) 

RADCI,l) = 2XPI*SQRTCXF0UNDX*2 + YF0UNDXX2) 

ELSE 

CHECK FOR VERTICAL LINE 
XFOUND = RCL) 

YFOUND * XFOUND / TANC THETA ) 

RADCI,L) = 2*PI*SQRTCXF0UND**2 + YF0UND**2) 

ENDIF 

IFCRADCI,L) . GT . RMAX) RMAX = RAD(1,L) 

IFCRADCI,L) . LT . RMIN) RMXN = RAD(I,l) 

33 CONTINUE 
22 CONTINUE 

****************************************************************** 

OUTPUT 

****************************************************************** 
WRITEC 2,104) IBIG,NL + 1 
WRITEC 2,102) RMIN, RMAX 
DO 44 I = 1 , IBIG 

WRITEC 2, 1 03) I,CRADCI,L),L=1,NL) 

44 CONTINUE 

****************************************************************** 

15 FGRMATC/8X, ' Please enter the height, H, of the V 

18X, ' cylinder Cfrom base to base). ' ) 

16 FORMATC/8X, 'Please enter the radius R, of the */ 

18X, ' cylinder. ' ) 

100 FORMAT C A ) 

101 FORMAT C 15) 

102 FORMAT ( ACE1 A .6 ) ) 

103 FORMAT CIS, 5C El A. 6)) 

1 0 A FORMAT C 2C 15) ) 

105 FORMAT C ' f ,A) 

110 F0RMATC//7X, 'TWO FILES ARE REQUIRED TO HOLD PROGRAM 0UPUT',//7X, 

1 ' DAT A . PLEASE ENTER THE NAME OF THE FILE TO BE USED TO ',//7X, 
2'HOLD THE PROGRAM MATERIAL PARAMETER DATA. PLEASE ',//7X, 

3 « INCL UDE THE EXTENSION. I. E. FILENAME . DAT ' ,//7X) 

111 FORMATC//7X, 'PLEASE ENTER THE NAME OF THE FILE TO BE USED ',//7X, 

1 'TO HOLD THE PROGRAM OUTPUT DATA. PLEASE INCLUDE THE',//7X, 

2 ' EXTENSION . I. E. FI L ENAME . DAT ' , //7X) 

112 FORMATC//7X, 'THE HEADERS ALLOW THE USER TO IDENTIFY THIS',//7X, 
I'SET OF DATA FROM ALL OTHER SETS.', 

2//7X, 'PLEASE ENTER HEADER #1 C6 A CHARACTERS MAX)',//7X) 
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113 F0RMATC//7X, 'PLEASE ENTER HEADER f 2 (64 CHARACTERS MAX)',//7X) 

114 FORMAT (//7X, 'ENTER THE NUMBER OF LAYERS C . L E . 5)»,//7X) 

CLOSE(l) 

CL0SEC2) 

STOP 

END 
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APPENDIX M 



SOFTWARE SOURCES 



Ryan-McFarland FORTRAN 

609 Deep Valley Drive 

Rolling Hills Estates, CA 90274 

(213) 541 - 4828 

Microcoxnpatibles, Inc. GRAFMATIC 
301 Prelude Drive 
Silver Spring, MD 20901 
(301) 593 - 5151 

Jewell Technologies, Inc. GRAFPLUS 
Seattle, WA 
(206) 937 - 1081 

West Coast Consultants CURVE DIGITIZER 
4202 Genesee Avenue, Suite 309 
San Diego, CA 92117 
(619) 565 - 1266 

Microsoft FORTRAN 

16011 NE 36th Way 
Box 97017 

Redmond, WA 98073 - 9717 

Prof. M. A. Morgan 
Code 62Mw 

Naval Postgraduate School 
Monterey, CA 93943 

LT. E. M. Connolly 

c/o J. J. Connolly 
208 Kent Drive 
Greenville, NC 27858 
(919) 752 - 2085 

OR 

c/o C. W. Snell 
228 Windsor Road 
Greenville, NC 27858 
(919) 756 - 4918 
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